JP3279513B2 - Bone density measurement device, bone density measurement method, and bone density measurement phantom - Google Patents
Bone density measurement device, bone density measurement method, and bone density measurement phantomInfo
- Publication number
- JP3279513B2 JP3279513B2 JP32629297A JP32629297A JP3279513B2 JP 3279513 B2 JP3279513 B2 JP 3279513B2 JP 32629297 A JP32629297 A JP 32629297A JP 32629297 A JP32629297 A JP 32629297A JP 3279513 B2 JP3279513 B2 JP 3279513B2
- Authority
- JP
- Japan
- Prior art keywords
- subject
- bone
- projection data
- bone density
- calibration
- 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 - Fee Related
Links
- 230000037182 bone density Effects 0.000 title claims description 162
- 238000001739 density measurement Methods 0.000 title claims description 69
- 238000000034 method Methods 0.000 title claims description 43
- 210000000988 bone and bone Anatomy 0.000 claims description 115
- 238000011088 calibration curve Methods 0.000 claims description 61
- 229910052588 hydroxylapatite Inorganic materials 0.000 claims description 38
- XYJRXVWERLGGKC-UHFFFAOYSA-D pentacalcium;hydroxide;triphosphate Chemical compound [OH-].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O XYJRXVWERLGGKC-UHFFFAOYSA-D 0.000 claims description 38
- 238000005259 measurement Methods 0.000 claims description 30
- 230000001054 cortical effect Effects 0.000 claims description 26
- 230000002093 peripheral effect Effects 0.000 claims description 23
- 238000013480 data collection Methods 0.000 claims description 14
- 239000000126 substance Substances 0.000 claims description 12
- 229910052586 apatite Inorganic materials 0.000 claims 2
- VSIIXMUUUJUKCM-UHFFFAOYSA-D pentacalcium;fluoride;triphosphate Chemical compound [F-].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O VSIIXMUUUJUKCM-UHFFFAOYSA-D 0.000 claims 2
- 125000002887 hydroxy group Chemical group [H]O* 0.000 claims 1
- 229910052500 inorganic mineral Inorganic materials 0.000 claims 1
- 239000011707 mineral Substances 0.000 claims 1
- 241000700159 Rattus Species 0.000 description 13
- 238000012937 correction Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 208000001132 Osteoporosis Diseases 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 238000013425 morphometry Methods 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 210000000689 upper leg Anatomy 0.000 description 3
- 208000020084 Bone disease Diseases 0.000 description 2
- 239000011575 calcium Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000009806 oophorectomy Methods 0.000 description 2
- 210000001672 ovary Anatomy 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 208000030136 Marchiafava-Bignami Disease Diseases 0.000 description 1
- 208000029725 Metabolic bone disease Diseases 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 238000000333 X-ray scattering Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003562 morphometric effect Effects 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 229920005749 polyurethane resin Polymers 0.000 description 1
- 208000001685 postmenopausal osteoporosis Diseases 0.000 description 1
- 229920002803 thermoplastic polyurethane Polymers 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
Description
【０００１】[0001]
【発明の属する技術分野】本発明は、Ｘ線ＣＴ装置を用
いた骨診断技術に関するものであり、特に、高齢化に伴
い生じる骨粗鬆症等の骨疾患の診断に好適に用いられる
骨密度測定装置および骨密度測定方法ならびに骨密度測
定用ファントムに関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a bone diagnosis technique using an X-ray CT apparatus, and more particularly to a bone density measuring apparatus suitably used for diagnosing a bone disease such as osteoporosis caused by aging. The present invention relates to a bone density measurement method and a bone density measurement phantom.
【０００２】[0002]
【従来の技術】従来より、ｐＱＣＴ（peripheral Quant
itative Computed Tomography ）装置を用いて骨の診断
がなされている。このｐＱＣＴ装置は、Ｘ線ＣＴの技術
を応用したものであり、被検体である骨の投影データを
収集し、この投影データに基づいてＣＴ画像を再構成す
るものである。ｐＱＣＴ装置を用いて被検体を測定すれ
ば、被検体の中心部にある海綿骨およびその周囲の皮質
骨それぞれの平均的な骨密度を求めることができる。し
かし、被検体の海綿骨の微細構造（骨梁構造）を観察す
ることはできない。2. Description of the Related Art Conventionally, pQCT (peripheral quantum
Bone diagnosis has been made using an itative computed tomography device. This pQCT apparatus is an application of X-ray CT technology, which collects projection data of a bone as a subject and reconstructs a CT image based on the projection data. If the subject is measured using the pQCT apparatus, the average bone density of the cancellous bone at the center of the subject and the cortical bone around the cannula can be determined. However, the microstructure (trabecular structure) of the cancellous bone of the subject cannot be observed.
【０００３】一方、被検体の海綿骨の微細構造は、非脱
灰薄切標本を用いた骨形態計測法により評価することが
可能である。この評価方法は、高度の専門性と長い評価
時間を要し、また、非破壊検査ではない。On the other hand, the microstructure of cancellous bone of a subject can be evaluated by a bone morphometry using a non-decalcified sliced sample. This evaluation method requires a high degree of expertise and long evaluation time, and is not a nondestructive inspection.
【０００４】ところで、骨粗鬆症や代謝性骨疾患に関す
る骨診断においては、海綿骨の一定部位の骨密度および
微細構造の双方を同時に且つ非破壊的に測定することが
要求されている（例えば、日本骨形態計測学会雑誌，
６，２４３−２５１（１９９６）を参照）。しかしなが
ら、上記のｐＱＣＴおよび非脱灰薄切標本を用いた骨形
態計測法の何れも、このような要求を満たすことはでき
ない。[0004] In the bone diagnosis for osteoporosis and metabolic bone disease, it is required to simultaneously and non-destructively measure both the bone density and the microstructure of a certain portion of cancellous bone (for example, Japanese bone). Journal of the Society of Morphometrics,
6, 243-251 (1996)). However, neither the pQCT nor the bone morphometry using non-decalcified slices can satisfy such a requirement.
【０００５】そこで、以下のような骨密度測定方法が採
られている。すなわち、ポリウレタン樹脂中に既知濃度
のハイドロキシアパタイトを含有させたロッドからなる
骨密度測定用ファントムであって当該ハイドロキシアパ
タイト濃度が互いに異なるものを複数個用意し、この複
数個の骨密度測定用ファントムと被検体とを同時にＸ線
ＣＴ装置により測定して投影データを収集し、この投影
データに基づいてＣＴ画像を再構成し、このＣＴ画像に
おける複数個の骨密度測定用ファントムおよび被検体そ
れぞれのＣＴ値の比較に基づいて被検体の骨密度を測定
する。Therefore, the following method of measuring bone density has been adopted. That is, a plurality of bone density measurement phantoms comprising a rod containing a known concentration of hydroxyapatite in a polyurethane resin and having different hydroxyapatite concentrations are prepared. The subject and the subject are simultaneously measured by an X-ray CT apparatus to collect projection data, and a CT image is reconstructed based on the projection data. The bone density of the subject is measured based on the comparison of the values.
【０００６】[0006]
【発明が解決しようとする課題】しかしながら、上記従
来の骨密度測定用ファントムを用いた骨密度測定方法で
は以下のような問題点がある。第１に、複数個の骨密度
測定用ファントムと被検体とを同時に測定する為に測定
視野を広く採る必要があることから、ＣＴ画像の分解能
が悪く、被検体の海綿骨の微細構造を観測することがで
きない。第２に、ＣＴ画像における骨密度測定用ファン
トム内のロッドのＣＴ値に基づいて被検体の海綿骨の骨
密度を推定するという方法が採られているが、Ｘ線ＣＴ
装置による測定時において被検体の高骨密度部分である
皮質骨でのＸ線散乱等の影響が考慮されていないので、
海綿骨のＣＴ値が実際の値より高く得られ、海綿骨の骨
密度を正確に測定することが困難である。However, the conventional bone density measuring method using the conventional bone density measuring phantom has the following problems. First, since it is necessary to adopt a wide measurement field of view to simultaneously measure a plurality of bone density measurement phantoms and the subject, the resolution of the CT image is poor, and the microstructure of the cancellous bone of the subject is observed. Can not do it. Second, a method of estimating the bone density of cancellous bone of a subject based on the CT value of a rod in a bone density measurement phantom in a CT image has been adopted.
Since the influence of X-ray scattering and the like on cortical bone, which is a high bone density part of the subject, is not taken into account during measurement by the device,
The CT value of cancellous bone is higher than the actual value, and it is difficult to accurately measure the bone density of cancellous bone.
【０００７】本発明は、上記問題点を解消する為になさ
れたものであり、被検体である骨の海綿骨の一定部位の
骨密度および微細構造の双方を同時に、非破壊的に且つ
正確に測定することができる骨密度測定装置および骨密
度測定方法を提供することを目的とする。また、この骨
密度測定装置および骨密度測定方法による骨密度の測定
に際して好適に用いられ得る骨密度測定用ファントムを
提供することをも目的とする。SUMMARY OF THE INVENTION The present invention has been made to solve the above problems, and simultaneously, non-destructively and accurately measures both the bone density and the microstructure of a certain portion of the cancellous bone of a subject bone. An object of the present invention is to provide a bone density measurement device and a bone density measurement method capable of measuring. It is another object of the present invention to provide a bone density measurement phantom that can be suitably used in measuring bone density by the bone density measurement device and the bone density measurement method.
【０００８】[0008]
【課題を解決するための手段】本発明に係る骨密度測定
装置は、(1) 定空間を挟んで配置されたＸ線管と検出器
アレイとを有し、測定空間内に置かれた被検体にＸ線管
によりＸ線を照射し、検出器アレイに到達したＸ線の強
度を測定して、投影データを収集する投影データ収集手
段と、(2) 投影データに基づいて画像を再構成する画像
再構成手段と、(3) 既知であって互いに異なる濃度のハ
イドロキシアパタイトを含有する円柱形状の中心部と中
心部の外周を皮質骨と等価な物質で覆う外周部とを有す
る２以上の所定数の骨密度測定用ファントムそれぞれに
ついて順次に投影データ収集手段により収集された校正
用投影データに基づいて画像再構成手段により再構成さ
れた校正用ＣＴ画像に基づいて、校正用ＣＴ画像の中心
部のＣＴ値と中心部のハイドロキシアパタイト濃度との
間の検量線を作成する検量線作成手段と、(4) 被検体に
ついて投影データ収集手段により収集された被検体投影
データに基づいて画像再構成手段により再構成された被
検体ＣＴ画像の海綿骨のＣＴ値および検量線に基づい
て、被検体の海綿骨の骨密度を獲得する被検体骨密度獲
得手段と、を備えることを特徴とする。According to the present invention, there is provided a bone density measuring apparatus comprising: (1) an X-ray tube and a detector array arranged with a fixed space interposed therebetween; Projection data collecting means for irradiating a sample with X-rays through an X-ray tube, measuring the intensity of the X-rays reaching the detector array, and collecting projection data, and (2) reconstructing an image based on the projection data Image reconstruction means, and (3) at least two or more having a cylindrical central portion containing a known and different concentration of hydroxyapatite and an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone Based on the calibration CT images reconstructed by the image reconstruction means based on the calibration projection data sequentially acquired by the projection data acquisition means for each of the predetermined number of bone density measurement phantoms, the center of the calibration CT image is calculated. Part CT value and central part A calibration curve creating means for creating a calibration curve between the concentration of hydroxyapatite and (4) an object reconstructed by the image reconstructing means based on the subject projection data collected by the projection data collecting means for the subject. And a subject bone density acquiring means for acquiring the bone density of the cancellous bone of the subject based on the CT value and the calibration curve of the cancellous bone of the subject CT image.
【０００９】この骨密度測定装置によれば、既知であっ
て互いに異なる濃度のハイドロキシアパタイトを含有す
る円柱形状の中心部と中心部の外周を皮質骨と等価な物
質で覆う外周部とを有する２以上の所定数の骨密度測定
用ファントムそれぞれについて、順次に投影データ収集
手段により校正用投影データが収集され、この校正用投
影データに基づいて画像再構成手段により校正用ＣＴ画
像が再構成され、そして、検量線作成手段により、この
校正用ＣＴ画像に基づいて、校正用ＣＴ画像の中心部の
ＣＴ値と中心部のハイドロキシアパタイト濃度との間の
検量線が作成される。一方、被検体について投影データ
収集手段により被検体投影データが収集され、この被検
体投影データに基づいて画像再構成手段により被検体Ｃ
Ｔ画像が再構成され、そして、被検体骨密度獲得手段に
より、この被検体ＣＴ画像の海綿骨のＣＴ値および検量
線に基づいて、被検体の海綿骨の骨密度が獲得される。According to this bone density measuring device, there is provided a cylindrical portion having a known and different concentration of hydroxyapatite and having an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone. For each of the predetermined number of bone density measurement phantoms, projection data for calibration is sequentially collected by the projection data collection unit, and a CT image for calibration is reconstructed by the image reconstruction unit based on the projection data for calibration. Then, based on the calibration CT image, a calibration curve is created between the CT value at the center of the calibration CT image and the hydroxyapatite concentration at the center by the calibration curve creating means. On the other hand, projection data of the subject is collected by the projection data collecting means, and the image reconstruction means based on the subject projection data.
The T image is reconstructed, and the subject bone density obtaining means obtains the bone density of the cancellous bone of the subject based on the CT value and the calibration curve of the cancellous bone of the subject CT image.
【００１０】また、本発明に係る骨密度測定装置では、
検量線作成手段は、所定数の骨密度測定用ファントムの
校正用ＣＴ画像それぞれの外周部のＣＴ値で補正した中
心部のＣＴ値と中心部のハイドロキシアパタイト濃度と
の間の検量線を作成するとともに、被検体骨密度獲得手
段は、被検体ＣＴ画像の皮質骨のＣＴ値で補正した海綿
骨のＣＴ値および検量線に基づいて被検体の海綿骨の骨
密度を獲得する、ことを特徴とする。この場合には、検
量線作成手段により、所定数の骨密度測定用ファントム
それぞれの校正用ＣＴ画像の中心部のＣＴ値が補正さ
れ、この補正ＣＴ値とハイドロキシアパタイト濃度との
間の検量線が作成され、そして、被検体骨密度獲得手段
により、被検体ＣＴ画像の海綿骨の補正されたＣＴ値お
よび検量線に基づいて被検体の海綿骨の骨密度が獲得さ
れる。[0010] In the bone density measuring apparatus according to the present invention,
The calibration curve creating means creates a calibration curve between the central CT value and the central hydroxyapatite concentration corrected by the CT value of the outer peripheral portion of each of the calibration CT images of the predetermined number of bone density measurement phantoms. In addition, the subject bone density acquiring means acquires the bone density of the cancellous bone of the subject based on the CT value of the cancellous bone and the calibration curve corrected with the CT value of the cortical bone of the subject CT image. I do. In this case, the calibration curve creating means corrects the CT value of the central part of the calibration CT image of each of the predetermined number of bone density measurement phantoms, and calculates a calibration curve between the corrected CT value and the hydroxyapatite concentration. The created bone density of the subject is acquired by the subject bone density acquiring means based on the corrected CT value and the calibration curve of the cancellous bone of the subject CT image.
【００１１】また、本発明に係る骨密度測定装置では、
投影データ収集手段は、被検体のスライス面に垂直な方
向に一定間隔の各位置それぞれで被検体を測定し、当該
測定結果を加算したものを被検体投影データとする、こ
とを特徴とする。また、投影データ収集手段は、被検体
のスライス面に垂直な方向の各位置それぞれで被検体を
測定し、当該測定結果を補間し積分したものを被検体投
影データとする、ことを特徴とする。これら何れの場合
にも、被検体の移動範囲に亘る平均された測定結果が得
られる。In the bone density measuring device according to the present invention,
The projection data collection means measures the subject at each position at regular intervals in a direction perpendicular to the slice plane of the subject, and adds the measured results to the projection data as subject projection data. Further, the projection data collecting means measures the subject at each position in a direction perpendicular to the slice plane of the subject, and interpolates and integrates the measurement result to obtain subject projection data. . In each case, an averaged measurement result over the moving range of the subject is obtained.
【００１２】本発明に係る骨密度測定方法は、(1) 既知
であって互いに異なる濃度のハイドロキシアパタイトを
含有する円柱形状の中心部と中心部の外周を皮質骨と等
価な物質で覆う外周部とを有する２以上の所定数の骨密
度測定用ファントムそれぞれを順次に、Ｘ線ＣＴ装置の
測定空間内の所定位置に配置し、Ｘ線ＣＴ装置により測
定して校正用投影データを収集する校正用投影データ収
集ステップと、(2) 所定数の骨密度測定用ファントムそ
れぞれについて、校正用投影データに基づいて校正用Ｃ
Ｔ画像を再構成する校正用ＣＴ画像再構成ステップと、
(3) 所定数の骨密度測定用ファントムの校正用ＣＴ画像
に基づいて中心部のＣＴ値と中心部のハイドロキシアパ
タイト濃度との間の検量線を作成する検量線作成ステッ
プと、(4) 被検体である骨を所定位置に配置し、Ｘ線Ｃ
Ｔ装置により測定して被検体投影データを収集する被検
体投影データ収集ステップと、(5) 被検体投影データに
基づいて被検体ＣＴ画像を再構成する被検体ＣＴ画像再
構成ステップと、(6) 被検体ＣＴ画像の海綿骨のＣＴ値
および検量線に基づいて被検体の海綿骨の骨密度を獲得
する被検体骨密度獲得ステップと、を備えることを特徴
とする。The method for measuring bone density according to the present invention comprises the steps of: (1) a center portion of a columnar shape containing a known and different concentration of hydroxyapatite and an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone; A calibration for sequentially arranging two or more predetermined number of bone density measurement phantoms each having a predetermined position in a measurement space of the X-ray CT apparatus and measuring the X-ray CT apparatus to collect calibration projection data (2) for each of the predetermined number of bone density measurement phantoms, based on the calibration projection data,
A calibration CT image reconstruction step for reconstructing a T image;
(3) a calibration curve creating step of creating a calibration curve between the central CT value and the central hydroxyapatite concentration based on the calibration CT images of a predetermined number of bone density measurement phantoms; The bone, which is the specimen, is placed at a predetermined position, and X-ray C
(5) a subject CT image reconstructing step of reconstructing a subject CT image based on the subject projection data, and (6) a subject CT image reconstructing step of reconstructing a subject CT image based on the subject projection data. A) obtaining a bone density of the cancellous bone of the subject based on the CT value of the cancellous bone of the subject CT image and a calibration curve.
【００１３】この骨密度測定方法によれば、既知であっ
て互いに異なる濃度のハイドロキシアパタイトを含有す
る円柱形状の中心部と中心部の外周を皮質骨と等価な物
質で覆う外周部とを有する２以上の所定数の骨密度測定
用ファントムそれぞれについて順次に、校正用投影デー
タ収集ステップにおいて校正用投影データが収集され、
この校正用投影データに基づいて校正用ＣＴ画像再構成
ステップにおいて校正用ＣＴ画像が再構成され、そし
て、検量線作成ステップにおいて、この校正用ＣＴ画像
に基づいて、校正用ＣＴ画像の中心部のＣＴ値と中心部
のハイドロキシアパタイト濃度との間の検量線が作成さ
れる。一方、被検体について、被検体投影データ収集ス
テップにおいて被検体投影データが収集され、この被検
体投影データに基づいて被検体ＣＴ画像再構成ステップ
において被検体ＣＴ画像が再構成され、そして、被検体
骨密度獲得ステップにおいて、この被検体ＣＴ画像の海
綿骨のＣＴ値および検量線に基づいて、被検体の海綿骨
の骨密度が獲得される。According to the method for measuring bone density, the method has a central portion of a columnar shape containing a known and different concentration of hydroxyapatite and an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone. For each of the above predetermined number of bone density measurement phantoms, calibration projection data is collected in the calibration projection data collection step,
Based on the calibration projection data, a calibration CT image is reconstructed in a calibration CT image reconstruction step, and in a calibration curve creating step, a central portion of the calibration CT image is constructed based on the calibration CT image. A calibration curve is created between the CT value and the central hydroxyapatite concentration. On the other hand, for the subject, subject projection data is collected in a subject projection data collection step, and a subject CT image is reconstructed in a subject CT image reconstruction step based on the subject projection data. In the bone density obtaining step, the bone density of the cancellous bone of the subject is obtained based on the CT value of the cancellous bone and the calibration curve of the CT image of the subject.
【００１４】また、本発明に係る骨密度測定方法では、
検量線作成ステップは、所定数の骨密度測定用ファント
ムの校正用ＣＴ画像それぞれの外周部のＣＴ値で補正し
た中心部のＣＴ値と中心部のハイドロキシアパタイト濃
度との間の検量線を作成するとともに、被検体骨密度獲
得ステップは、被検体ＣＴ画像の皮質骨のＣＴ値で補正
した海綿骨のＣＴ値および検量線に基づいて被検体の海
綿骨の骨密度を獲得する、ことを特徴とする。この場合
には、検量線作成ステップにより、所定数の骨密度測定
用ファントムそれぞれの校正用ＣＴ画像の中心部のＣＴ
値が補正され、この補正ＣＴ値とハイドロキシアパタイ
ト濃度との間の検量線が作成され、そして、被検体骨密
度獲得ステップにより、被検体ＣＴ画像の海綿骨の補正
されたＣＴ値および検量線に基づいて被検体の海綿骨の
骨密度が獲得される。In the method for measuring bone density according to the present invention,
The calibration curve creating step creates a calibration curve between the central CT value and the central hydroxyapatite concentration corrected by the CT value of the outer peripheral portion of each of the calibration CT images of the predetermined number of bone density measurement phantoms. In addition, the subject bone density acquiring step acquires the bone density of the cancellous bone of the subject based on the CT value of the cancellous bone and the calibration curve corrected with the CT value of the cortical bone of the subject CT image. I do. In this case, in the calibration curve creating step, the CT of the central part of the calibration CT image of each of the predetermined number of bone density measurement phantoms is determined.
The value is corrected, a calibration curve is created between the corrected CT value and the hydroxyapatite concentration, and the subject bone density acquisition step is performed to obtain the corrected CT value and calibration curve of the cancellous bone of the subject CT image. Based on this, the bone density of the cancellous bone of the subject is obtained.
【００１５】また、本発明に係る骨密度測定方法では、
被検体投影データ収集ステップは、被検体のスライス面
に垂直な方向に一定間隔の各位置それぞれで被検体を測
定し、当該測定結果を加算したものを被検体投影データ
とする、ことを特徴とする。また、被検体投影データ収
集ステップは、被検体のスライス面に垂直な方向の各位
置それぞれで被検体を測定し、当該測定結果を補間し積
分したものを被検体投影データとする、ことを特徴とす
る。これら何れの場合にも、被検体の移動範囲に亘る平
均された測定結果が得られる。Further, in the method for measuring bone density according to the present invention,
The subject projection data collection step measures the subject at each position at regular intervals in a direction perpendicular to the slice plane of the subject, and the sum of the measurement results is used as subject projection data. I do. Further, the subject projection data collection step measures the subject at each position in a direction perpendicular to the slice plane of the subject, and interpolates and integrates the measurement result to obtain subject projection data. And In each case, an averaged measurement result over the moving range of the subject is obtained.
【００１６】本発明に係る骨密度測定用ファントムは、
既知の濃度のハイドロキシアパタイトを含有する円柱形
状の中心部と、中心部の外周を皮質骨と等価な物質で覆
う外周部と、を備えることを特徴とする。この骨密度測
定用ファントムは、上記の本発明に係る骨密度測定用装
置および方法による骨密度の測定に好適に用いられる。The phantom for measuring bone density according to the present invention comprises:
It is characterized by comprising a cylindrical central portion containing a known concentration of hydroxyapatite, and an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone. This bone density measurement phantom is suitably used for the measurement of bone density by the above-described apparatus and method for measuring bone density according to the present invention.
【００１７】[0017]
【発明の実施の形態】以下、添付図面を参照して本発明
の実施の形態を詳細に説明する。尚、図面の説明におい
て同一の要素には同一の符号を付し、重複する説明を省
略する。Embodiments of the present invention will be described below in detail with reference to the accompanying drawings. In the description of the drawings, the same elements will be denoted by the same reference symbols, without redundant description.
【００１８】先ず、本発明に係る骨密度測定用ファント
ムの１実施形態について図１を用いて説明する。図１
は、本実施形態に係る骨密度測定用ファントム１０を説
明する斜視図である。この骨密度測定用ファントム１０
は、中心部１１とその周囲の外周部１２とからなる円柱
形状のものである。First, one embodiment of a bone density measurement phantom according to the present invention will be described with reference to FIG. FIG.
FIG. 1 is a perspective view illustrating a bone density measurement phantom 10 according to the present embodiment. This bone density measurement phantom 10
Has a columnar shape including a central portion 11 and an outer peripheral portion 12 around the central portion.
【００１９】骨密度測定用ファントム１０の中心領域に
ある中心部１１は、直径が例えば４ｍｍφの円柱形状の
ウレタン樹脂にハイドロキシアパタイトを所定の濃度で
均一に含有させたものであり、そのハイドロキシアパタ
イト濃度は既知である。このハイドロキシアパタイト
は、海綿骨の成分比率と類似する成分比率でＣａ（カル
シウム）元素およびＰ（リン）元素を含むものである。The central portion 11 in the central region of the bone density measurement phantom 10 is made of a columnar urethane resin having a diameter of, for example, 4 mmφ, in which hydroxyapatite is uniformly contained at a predetermined concentration. Is known. This hydroxyapatite contains a Ca (calcium) element and a P (phosphorus) element in a component ratio similar to the component ratio of cancellous bone.
【００２０】この中心部１１の外周にある外周部１２
は、被検体の皮質骨と等価な物質で中心部１１を覆った
ものである。この等価な物質として、例えばラットの大
腿骨幹部が用いられる。具体的には、ラット大腿骨幹部
をハンドモータで４分割から６分割程度に縦断し、これ
を中心部１１の外周にワックスで貼り付けたものを外周
部１２とすることができる。An outer peripheral portion 12 on the outer periphery of the central portion 11
Is a material in which the central part 11 is covered with a substance equivalent to the cortical bone of the subject. For example, a rat femoral shaft is used as this equivalent substance. Specifically, the rat femoral shaft can be longitudinally cut into four to six parts by a hand motor, and this can be attached to the outer periphery of the central part 11 with wax to form the outer peripheral part 12.
【００２１】そして、本発明に係る骨密度測定装置およ
び骨密度測定方法による被検体の海綿骨の骨密度測定に
際しては、このような骨密度測定用ファントム１０を複
数個用意する。この複数個の骨密度測定用ファントムそ
れぞれは、中心部１１に含有されるハイドロキシアパタ
イトの濃度が既知であって互いに異なるものであり、そ
のハイドロキシアパタイトの濃度範囲は、被検体の海綿
骨の骨密度範囲をカバーするものとする。なお、必要な
骨密度測定用ファントムの個数は、少なくとも２個であ
るが、好適には３個以上である。When measuring the bone density of cancellous bone of a subject using the bone density measuring apparatus and the bone density measuring method according to the present invention, a plurality of such bone density measuring phantoms 10 are prepared. Each of the plurality of bone density measurement phantoms has a known concentration of hydroxyapatite contained in the central portion 11 and is different from each other. The concentration range of the hydroxyapatite is determined by the bone density of cancellous bone of the subject. Shall cover the range. The required number of bone density measurement phantoms is at least two, but preferably three or more.
【００２２】次に、本発明に係る骨密度測定装置の１実
施形態について、図２を用いて説明する。図２は、本実
施形態に係る骨密度測定装置の構成図である。この骨密
度測定装置は、Ｘ線ＣＴ装置１００、制御演算装置２０
０、画像処理装置３１０およびＸ線管制御装置３２０を
備えて構成される。Next, an embodiment of a bone density measuring apparatus according to the present invention will be described with reference to FIG. FIG. 2 is a configuration diagram of the bone density measurement device according to the present embodiment. This bone density measuring device includes an X-ray CT device 100, a control arithmetic device 20
0, an image processing device 310 and an X-ray tube control device 320.
【００２３】Ｘ線ＣＴ装置１００は、Ｘ線を遮蔽する部
材からなる測定室１０１内に、Ｘ線管１１０、回転ステ
ージ１２１、リニアＸ線検出器１３０およびＸ線イメー
ジインテンシファイア１４０等を備えている。Ｘ線管１
１０は、Ｘ線管制御装置３２０により所定の電圧が供給
されてＸ線を出力するものである。このＸ線管１１０と
して、特に、極めて微小な領域からＸ線を発生するマイ
クロフォーカスＸ線管が、回転ステージ１２１の上に載
置された被検体２０の微細構造を測定する為に好適に用
いられる。The X-ray CT apparatus 100 includes an X-ray tube 110, a rotary stage 121, a linear X-ray detector 130, an X-ray image intensifier 140, and the like in a measurement room 101 made of a member for shielding X-rays. ing. X-ray tube 1
Reference numeral 10 denotes a device to which a predetermined voltage is supplied by the X-ray tube controller 320 to output X-rays. As the X-ray tube 110, in particular, a microfocus X-ray tube that generates X-rays from an extremely small area is preferably used for measuring the fine structure of the subject 20 placed on the rotating stage 121. Can be
【００２４】被検体２０（または、骨密度測定用ファン
トム１０）を載置する回転ステージ１２１はＺ軸移動ス
テージ１２４の上に設けられ、Ｚ軸移動ステージ１２４
はＸ軸移動ステージ１２２の上に設けられ、また、Ｘ軸
移動ステージ１２２はＹ軸移動ステージ１２３の上に設
けられている。ここで、Ｙ軸はＸ線管１１０とリニアＸ
線検出器１３０とを結ぶ方向であり、Ｘ軸およびＺ軸そ
れぞれはＹ軸に直交する方向であり、また、回転ステー
ジ１２１の回転中心はＺ軸である。The rotary stage 121 on which the subject 20 (or the phantom 10 for measuring bone density) is mounted is provided on a Z-axis moving stage 124.
Is provided on an X-axis moving stage 122, and the X-axis moving stage 122 is provided on a Y-axis moving stage 123. Here, the Y axis is the X-ray tube 110 and the linear X
The X-axis and the Z-axis are directions orthogonal to the Y-axis, and the rotation center of the rotary stage 121 is the Z-axis.
【００２５】Ｘ線管１１０から出力され被検体２０を経
たＸ線を検出するリニアＸ線検出器１３０は、Ｘ軸方向
に一次元状に配列された多数のＸ線検出素子からなり、
例えば、チャネル数が５１２で素子間ピッチが５０μｍ
の一次元型フォトダイオードアレイが好適に用いられ
る。このリニアＸ線検出器１３０の受光面の前に設けら
れたコリメータ１３１は、リニアＸ線検出器１３０の受
光面に散乱Ｘ線が入射するのを阻止するものである。こ
のリニアＸ線検出器１３０はＸ軸移動ステージ１３２の
上に設けられ、また、Ｘ軸移動ステージ１３２はＹ軸移
動ステージ１３３の上に設けられている。The linear X-ray detector 130 for detecting X-rays output from the X-ray tube 110 and passing through the subject 20 is composed of a large number of X-ray detecting elements arranged one-dimensionally in the X-axis direction.
For example, the number of channels is 512 and the pitch between elements is 50 μm.
Is preferably used. The collimator 131 provided in front of the light receiving surface of the linear X-ray detector 130 prevents scattered X-rays from entering the light receiving surface of the linear X-ray detector 130. The linear X-ray detector 130 is provided on an X-axis moving stage 132, and the X-axis moving stage 132 is provided on a Y-axis moving stage 133.
【００２６】リニアＸ線検出器１３０の後方にはＸ線イ
メージインテンシファイア１４０が設けられている。こ
のＸ線イメージインテンシファイア１４０は、リニアＸ
線検出器１３０がＸ軸移動ステージ１３２により待避状
態にあるときに、回転ステージ１２１の上に載置された
被検体２０等のＸ線透視像を撮像するものであり、Ｙ軸
移動ステージ１４３の上に設けられている。Behind the linear X-ray detector 130, an X-ray image intensifier 140 is provided. The X-ray image intensifier 140 is a linear X-ray image intensifier.
When the X-ray detector 130 is in the retracted state by the X-axis moving stage 132, the X-ray moving unit 143 captures an X-ray fluoroscopic image of the subject 20 or the like placed on the rotating stage 121. It is provided above.
【００２７】また、Ｙ軸移動ステージ１２３，１３３お
よび１４３それぞれは、測定室１０１の床面に固定され
た架台１５０の上のステージ移動レール１５１の上をＹ
軸方向に移動可能である。これらのＹ軸方向の移動によ
り、リニアＸ線検出器１３０またはＸ線イメージインテ
ンシファイア１４０の撮像倍率は適切に設定され得る。Each of the Y-axis moving stages 123, 133 and 143 moves on a stage moving rail 151 on a gantry 150 fixed to the floor of the measuring chamber 101.
It is movable in the axial direction. By the movement in the Y-axis direction, the imaging magnification of the linear X-ray detector 130 or the X-ray image intensifier 140 can be appropriately set.
【００２８】Ｘ線イメージインテンシファイア１４０に
より撮像された被検体２０等のＸ線透視像は画像処理装
置３１０に送られ、画像処理装置３１０は、このＸ線透
視像に基づいて画像積算、コントラスト強調およびエッ
ジエンハンス等の画像処理を行う。An X-ray fluoroscopic image of the subject 20 or the like captured by the X-ray image intensifier 140 is sent to the image processing device 310, and the image processing device 310 performs image integration and contrast based on the X-ray fluoroscopic image. Image processing such as enhancement and edge enhancement is performed.
【００２９】制御演算装置２００は、Ｘ線ＣＴ装置１０
０内の各移動ステージを駆動制御し、Ｘ線管制御装置３
２０をも制御する。また、制御演算装置２００は、回転
ステージ１２１を回転駆動すると同時にリニアＸ線検出
器１３０から出力された投影データを入力する。回転ス
テージ１２１上に骨密度測定用ファントム１０が置かれ
ているときに得られた投影データを校正用投影データと
し、２以上の骨密度測定用ファントム１０それぞれにつ
いて校正用投影データを収集する。また、回転ステージ
１２１上に被検体２０である骨が置かれているときに得
られた投影データを被検体投影データとする。The control arithmetic unit 200 includes the X-ray CT apparatus 10
Drive control of each moving stage within the X-ray tube control device 3
20 is also controlled. Further, the control operation device 200 inputs the projection data output from the linear X-ray detector 130 at the same time as driving the rotation stage 121 to rotate. Projection data obtained when the bone density measurement phantom 10 is placed on the rotating stage 121 is used as calibration projection data, and calibration projection data is collected for each of the two or more bone density measurement phantoms 10. Further, projection data obtained when the bone as the subject 20 is placed on the rotating stage 121 is referred to as subject projection data.
【００３０】制御演算装置２００内の画像再構成部２１
０は、２以上の骨密度測定用ファントム１０それぞれに
ついて収集された校正用投影データに基づいて校正用Ｃ
Ｔ画像を再構成し、また、被検体２０について収集され
た被検体投影データに基づいて被検体ＣＴ画像を再構成
する。そして、表示部２４０は、これらの校正用ＣＴ画
像や被検体ＣＴ画像を画面上に表示する。The image reconstruction unit 21 in the control arithmetic unit 200
0 is a calibration C based on the calibration projection data collected for each of the two or more bone density measurement phantoms 10.
The T image is reconstructed, and the subject CT image is reconstructed based on the subject projection data collected for the subject 20. Then, the display unit 240 displays the calibration CT image and the subject CT image on the screen.
【００３１】制御演算装置２００内の検量線作成部２２
０は、画像再構成部２１０により再構成された校正用Ｃ
Ｔ画像に基づいて検量線を作成する。この検量線は、骨
密度測定用ファントム１０の中心部１１のＣＴ値とハイ
ドロキシアパタイト濃度との間の関係を表すものであ
る。制御演算装置２００内の被検体骨密度獲得部２３０
は、被検体ＣＴ画像の海綿骨のＣＴ値を上記の検量線に
当てはめて、被検体２０の海綿骨の骨密度を獲得する。
そして、表示部２４０は、被検体２０の海綿骨の微細構
造とともに骨密度の分布を画面上に例えば濃淡で表示す
る。なお、検量線作成および被検体骨密度獲得それぞれ
の方法の詳細については後述する。The calibration curve creating unit 22 in the control arithmetic unit 200
0 is the calibration C reconstructed by the image reconstruction unit 210
A calibration curve is created based on the T image. This calibration curve represents the relationship between the CT value of the central part 11 of the bone density measurement phantom 10 and the hydroxyapatite concentration. Subject bone density acquisition unit 230 in control arithmetic unit 200
Applies the CT value of the cancellous bone of the subject CT image to the above calibration curve to obtain the bone density of the cancellous bone of the subject 20.
Then, the display unit 240 displays the distribution of the bone density together with the fine structure of the cancellous bone of the subject 20 in, for example, shading on the screen. The details of each of the methods for preparing the calibration curve and obtaining the subject bone density will be described later.
【００３２】また、制御演算装置２００は、Ｚ軸移動ス
テージ１２４を駆動することにより、被検体２０のスラ
イス面に垂直な方向すなわちＺ軸方向に一定間隔の各位
置それぞれで収集した投影データを加算したものを上記
被検体投影データとするのも好適である。また、被検体
２０のスライス面に垂直な方向すなわちＺ軸方向の各位
置それぞれで収集した投影データを補間し積分したもの
を上記被検体投影データとするのも好適である。何れの
場合にも、被検体骨密度獲得部２３０により、被検体２
０の測定範囲に亘るＺ軸方向に関して平均された微細構
造と骨密度分布が獲得される。By driving the Z-axis moving stage 124, the control arithmetic unit 200 adds projection data collected at each position at regular intervals in the direction perpendicular to the slice plane of the subject 20, that is, in the Z-axis direction. It is also preferable that the result is used as the object projection data. It is also preferable that the projection data acquired by interpolating and integrating the projection data collected in each direction in the direction perpendicular to the slice plane of the subject 20, that is, in the Z-axis direction, be used as the subject projection data. In any case, the subject bone density acquisition unit 230 causes the subject 2
An averaged microstructure and bone density distribution over the Z-axis direction over a measurement range of zero is obtained.
【００３３】次に、この骨密度測定装置の作用ととも
に、本発明に係る骨密度測定方法の１実施形態につい
て、図３〜図６を用いて説明する。図３は、本実施形態
に係る骨密度測定方法を説明するフローチャート図であ
る。図４は、本実施形態に係る骨密度測定方法に用いら
れる複数個の骨密度測定用ファントムを説明する図表で
ある。図５は、本実施形態に係る骨密度測定方法におけ
る検量線を説明するグラフである。図６は、その検量線
に基づいて被検体の海綿骨の骨密度を獲得する方法を説
明するグラフである。Next, an operation of the bone density measuring apparatus and an embodiment of a bone density measuring method according to the present invention will be described with reference to FIGS. FIG. 3 is a flowchart illustrating the bone density measurement method according to the present embodiment. FIG. 4 is a table illustrating a plurality of bone density measurement phantoms used in the bone density measurement method according to the present embodiment. FIG. 5 is a graph illustrating a calibration curve in the bone density measurement method according to the present embodiment. FIG. 6 is a graph illustrating a method for obtaining the bone density of cancellous bone of a subject based on the calibration curve.
【００３４】この骨密度測定方法においては、先ずステ
ップＳ１（校正用投影データ収集ステップ）において、
上述したような複数個の骨密度測定用ファントム１０そ
れぞれを順次に、Ｘ線ＣＴ装置１００の測定空間内の回
転ステージ１２１上に載置し、回転ステージ１２１を回
転駆動させながらリニアＸ線検出器１３０により測定し
て、校正用投影データを制御演算装置２００により収集
する。なお、複数個の骨密度測定用ファントム１０は、
例えば図４に示すように中心部１１のハイドロキシアパ
タイト濃度がそれぞれ１００，２００，３００，…，１
００ｎ（ｍｇ／ｃｍ3 ）のものである（ただし、ｎは２
以上の整数）。In this bone density measurement method, first, in step S1 (calibration projection data collection step),
Each of the plurality of bone density measurement phantoms 10 described above is sequentially placed on a rotation stage 121 in a measurement space of the X-ray CT apparatus 100, and the linear X-ray detector is driven while rotating the rotation stage 121. The measurement is performed by 130 and the projection data for calibration is collected by the control arithmetic unit 200. The plurality of bone density measurement phantoms 10
For example, as shown in FIG. 4, the hydroxyapatite concentration in the central portion 11 is 100, 200, 300,.
00n (mg / cm 3 ) (where n is 2
Integer greater than or equal to).
【００３５】続くステップＳ２（校正用ＣＴ画像再構成
ステップ）において、ｎ個の骨密度測定用ファントム１
０それぞれについて、ステップＳ１で収集された校正用
投影データに基づいて校正用ＣＴ画像を画像再構成部２
１０により再構成する。続くステップＳ３（検量線作成
ステップ）において、ｎ個の骨密度測定用ファントム１
０の校正用ＣＴ画像それぞれの中心部１１のＣＴ値を外
周部１２のＣＴ値で補正した補正ＣＴ値と中心部１１の
ハイドロキシアパタイト濃度との間の検量線を検量線作
成部２２０により作成する。In the following step S2 (reconstructing CT image for calibration), n phantoms 1 for measuring bone density are used.
0 for each of the image reconstruction units 2 based on the calibration projection data collected in step S1.
10 to reconstruct. In the following step S3 (calibration curve creation step), the n bone density measurement phantoms 1
The calibration curve creation unit 220 creates a calibration curve between the corrected CT value obtained by correcting the CT value of the central part 11 of each of the zero calibration CT images with the CT value of the outer peripheral part 12 and the hydroxyapatite concentration of the central part 11. .
【００３６】この検量線作成の詳細は以下のとおりであ
る。ステップＳ２で再構成されたｎ個の骨密度測定用フ
ァントム１０の校正用ＣＴ画像それぞれにおいて、中心
部１１の平均ＣＴ値および外周部１２の平均ＣＴ値を図
４に示すようにｘk ，ａk とし（k=1,2,3,…,n）、これ
らの値に基づいてThe details of the preparation of the calibration curve are as follows. In each of the calibration CT images of the n bone density measurement phantoms 10 reconstructed in step S2, the average CT value of the central portion 11 and the average CT value of the outer peripheral portion 12 are represented by x k , a as shown in FIG. k (k = 1,2,3,…, n), and based on these values,
【数１】 なる式に従って補正を行い、中心部１１の補正ＣＴ値Ｘ
k を求める。そして、この中心部１１の補正ＣＴ値Ｘk
とハイドロキシアパタイト濃度との間の検量線を、図５
に示すように、（ｎ−１）次の回帰式または（ｎ−１）
次未満の回帰式として求める。(Equation 1) The correction is performed according to the following formula, and the corrected CT value X
Find k . Then, the corrected CT value X k of the central portion 11 is obtained.
The calibration curve between the concentration of hydroxyapatite and the concentration of hydroxyapatite is shown in FIG.
As shown in the following, (n-1) regression equation or (n-1)
Calculate as a regression equation less than
【００３７】なお、この (1)式に示す補正式は、或る媒
質を透過するＸ線の強度Ｉが媒質の吸収係数μ、密度ρ
および透過距離ｔの関数としてIt should be noted that the correction equation shown in the equation (1) indicates that the intensity I of the X-ray transmitted through a certain medium is equal to the absorption coefficient μ of the medium and the density ρ
And as a function of the transmission distance t
【数２】 なる式で表され、また、ＣＴ値が減衰量μρｔに応じた
量であることに基づいて導かれるものである。(Equation 2) And is derived based on the fact that the CT value is an amount corresponding to the amount of attenuation μρt.
【００３８】続くステップＳ４（被検体投影データ収集
ステップ）において、被検体２０である骨を回転ステー
ジ１２１の上に載置し、回転ステージ１２１を回転駆動
させながらリニアＸ線検出器１３０により測定して被検
体投影データを制御演算装置２００により収集する。ス
テップＳ５（被検体ＣＴ画像再構成ステップ）におい
て、この被検体投影データに基づいて被検体ＣＴ画像を
画像再構成部２１０により再構成する。そして、ステッ
プＳ６（被検体骨密度獲得ステップ）において、被検体
ＣＴ画像の海綿骨のＣＴ値を皮質骨のＣＴ値で補正した
補正ＣＴ値および検量線に基づいて被検体２０の海綿骨
の骨密度を、被検体骨密度獲得部２３０により獲得す
る。In the following step S4 (subject projection data collection step), the bone, which is the subject 20, is placed on the rotary stage 121, and is measured by the linear X-ray detector 130 while rotating the rotary stage 121. Then, the control processor 200 collects the subject projection data. In step S5 (subject CT image reconstruction step), the subject CT image is reconstructed by the image reconstruction unit 210 based on the subject projection data. In step S6 (subject bone density acquisition step), the trabecular bone of the subject 20 is corrected based on the corrected CT value obtained by correcting the CT value of the trabecular bone of the subject CT image with the CT value of the cortical bone and the calibration curve. The density is obtained by the subject bone density obtaining unit 230.
【００３９】この被検体２０の海綿骨の骨密度の獲得の
方法の詳細は以下のとおりである。ステップＳ５で再構
成された被検体ＣＴ画像において、被検体２０の海綿骨
のＣＴ値ｑおよび皮質骨の平均ＣＴ値ｒを求め、これら
の値に基づいて、上記 (1)式と等価な補正式であるThe details of the method for obtaining the bone density of the cancellous bone of the subject 20 are as follows. In the subject CT image reconstructed in step S5, the CT value q of the cancellous bone and the average CT value r of the cortical bone of the subject 20 are obtained, and a correction equivalent to the above equation (1) is made based on these values. Is an expression
【数３】 なる式に従って補正を行い、海綿骨の補正ＣＴ値Ｑを求
める。そして、図６に示すように、ステップＳ３で求め
た検量線に海綿骨の補正ＣＴ値Ｑを当てはめて、ハイド
ロキシアパタイト濃度で換算した骨密度ρを得る。(Equation 3) The correction is performed according to the following formula to obtain a corrected CT value Q of the cancellous bone. Then, as shown in FIG. 6, the corrected CT value Q of the cancellous bone is applied to the calibration curve obtained in step S3 to obtain the bone density ρ converted by the hydroxyapatite concentration.
【００４０】ここで、被検体２０の海綿骨の補正ＣＴ値
Ｑは、ステップＳ４における被検体投影データ収集の際
のスライス面上の位置の関数であり、被検体２０の海綿
骨の骨密度ρも、当該スライス面上の位置の関数であ
る。すなわち、以上のステップＳ１〜Ｓ６からなる骨密
度測定方法によれば、被検体２０である骨の海綿骨のス
ライス面上の微細構造を観測することができるだけでな
く、そのスライス面上の骨密度ρの分布を得ることがで
きる。Here, the corrected CT value Q of the cancellous bone of the subject 20 is a function of the position on the slice plane at the time of acquiring the projection data of the subject in step S4. Is also a function of the position on the slice plane. That is, according to the bone density measurement method including the above steps S1 to S6, not only can the microstructure on the slice surface of the cancellous bone of the subject 20 be observed, but also the bone density on the slice surface can be observed. The distribution of ρ can be obtained.
【００４１】なお、上記実施形態の説明においては、ス
テップＳ３において、ｎ個の骨密度測定用ファントム１
０それぞれの中心部１１の平均ＣＴ値ｘk を外周部１２
の平均ＣＴ値ａk で (1)式に従って補正して中心部１１
の補正ＣＴ値Ｘk を求め（k=1,2,3,…,n）、また、ステ
ップＳ６において、被検体２０の海綿骨のＣＴ値ｑを皮
質骨の平均ＣＴ値ｒで (3)式に従って補正し海綿骨の補
正ＣＴ値Ｑを求め、これにより、被検体２０の海綿骨の
骨密度を精度良く得ようとするものであった。しかし、
許容される誤差に依っては、上記補正を行うことなく、
ステップＳ３では複数個の骨密度測定用ファントム１０
それぞれの中心部１１の平均ＣＴ値ｘkそのままを用い
て検量線を作成し、ステップＳ４では被検体２０の海綿
骨のＣＴ値ｑそのままを用いて被検体２０の海綿骨の骨
密度ρを求めてもよい。In the description of the above embodiment, in step S3, the n bone density measurement phantoms 1
0 The average CT value x k of each central portion 11 is
Is corrected according to equation (1) with the average CT value a k of
Seeking compensation CT value X k (k = 1,2,3, ... , n), also in step S6, the CT value q cancellous bone of the subject 20 by the average CT value r of the cortical bone (3) The corrected CT value Q of the cancellous bone was obtained by performing correction according to the formula, and thereby, the bone density of the cancellous bone of the subject 20 was to be obtained with high accuracy. But,
Depending on the permissible error, without performing the above correction,
In step S3, the plurality of bone density measurement phantoms 10
A calibration curve is created using the average CT value x k of each central part 11 as it is, and in step S4, the bone density ρ of the cancellous bone of the subject 20 is obtained using the CT value q of the cancellous bone of the subject 20 as it is. You may.
【００４２】また、上記実施形態の説明においては、被
検体２０の１つのスライス面上の海綿骨の骨密度の測定
について説明したが、被検体２０の２以上のスライス面
について測定して被検体投影データを求めてもよい。In the above embodiment, the measurement of the bone density of cancellous bone on one slice surface of the subject 20 has been described. Projection data may be obtained.
【００４３】例えば、ステップＳ４で、被検体２０のス
ライス面に垂直な方向すなわちＺ軸方向に一定間隔の各
位置それぞれで被検体２０を測定し、当該測定結果を加
算したものを被検体投影データとしてもよい。すなわ
ち、図７に示すように、海綿骨２１とその周囲の皮質骨
２２からなる被検体２０を模式的に円柱形状のものとし
て表すと、被検体２０のスライス面（図中でハッチング
で示す各面）に垂直な方向に一定間隔（この図では、ス
ライス厚ｕに等しい間隔）の各位置それぞれで被検体２
０を測定する。このとき、被検体２０をＺ軸移動ステー
ジ１２４により移動させてもよいし、あるいは、Ｘ線管
１１０およびリニアＸ線検出器１３０を一体としてＺ軸
方向に移動させてもよい。そして、各位置で被検体２０
を測定した結果を加算したものを被検体投影データとす
る。ステップＳ５以降の内容は上述したものと同様であ
る。このようにすることにより、被検体２０の移動方向
に関して当該移動範囲における平均的な海綿骨２１の骨
密度が得られる。For example, in step S4, the subject 20 is measured at each position at regular intervals in the direction perpendicular to the slice plane of the subject 20, that is, in the Z-axis direction, and the result of the measurement is added to the subject projection data. It may be. That is, as shown in FIG. 7, when the subject 20 including the cancellous bone 21 and the cortical bone 22 around the same is schematically represented as a columnar shape, the slice surface of the subject 20 (each hatched portion in the drawing) The subject 2 at each position at a constant interval (in this figure, an interval equal to the slice thickness u) in the direction perpendicular to the
Measure 0. At this time, the subject 20 may be moved by the Z-axis moving stage 124, or the X-ray tube 110 and the linear X-ray detector 130 may be integrally moved in the Z-axis direction. At each position, the subject 20
The sum of the measurement results is used as subject projection data. The contents after step S5 are the same as those described above. In this way, an average bone density of the cancellous bone 21 in the moving range of the subject 20 in the moving direction can be obtained.
【００４４】あるいは、ステップＳ４で、被検体２０の
スライス面に垂直な方向すなわちＺ軸方向の各位置それ
ぞれで被検体２０を測定し、当該測定結果を補間し積分
したものを被検体投影データとしてもよい。すなわち、
図８に示すように、被検体２０のスライス面に垂直な方
向の位置Ｐ0 ，Ｐ1 ，Ｐ2 ，…，Ｐn それぞれにおいて
被検体２０を測定する。これらの位置Ｐ0 ，Ｐ1 ，Ｐ
2 ，…，Ｐn は、一定間隔であってもよいが、これに限
られるものではない。そして、図９に示すように、各位
置Ｐk で被検体２０を測定して得られたＣＴ値ＣＴk を
補間して位置Ｐk- 1 から位置Ｐk までの補間曲線ｆk(p)
を求め（k=0,1,2,…,n）、この補間曲線を位置Ｐ0 から
位置Ｐn までに亘ってAlternatively, in step S4, the subject 20 is measured at each position in the direction perpendicular to the slice plane of the subject 20, ie, in the Z-axis direction, and the measurement results are interpolated and integrated as subject projection data. Is also good. That is,
As shown in FIG. 8, the subject 20 is measured at positions P 0 , P 1 , P 2 ,..., P n in the direction perpendicular to the slice plane of the subject 20. These positions P 0 , P 1 , P
2 ,..., Pn may be at regular intervals, but are not limited thereto. Then, as shown in FIG. 9, an interpolation curve f k (p) from the position P k − 1 to the position P k is obtained by interpolating the CT value CT k obtained by measuring the subject 20 at each position P k.
The calculated (k = 0,1,2, ..., n ), the interpolation curve over from the position P 0 to the position P n
【数４】 の如く積分して、位置Ｐ0 から位置Ｐn までに亘る被検
体２０の海綿骨２１における平均ＣＴ値を求め、これを
被検体投影データとする。ステップＳ５以降の内容は上
述したものと同様である。このようにすることにより、
被検体２０の移動方向に関して当該移動範囲における平
均的な海綿骨２１の骨密度が得られる。(Equation 4) Integrated to as, an average CT value in cancellous bone 21 of the subject 20 over the position P 0 to the position P n, which is the subject projection data. The contents after step S5 are the same as those described above. By doing this,
With respect to the moving direction of the subject 20, an average bone density of the cancellous bone 21 in the moving range is obtained.
【００４５】次に、本発明に係る骨密度測定装置、骨密
度測定方法および骨密度測定用ファントムを用いた実施
例について説明する。なお、この実施例は本発明の範囲
を限定するものではない。Next, an embodiment using a bone density measuring apparatus, a bone density measuring method, and a bone density measuring phantom according to the present invention will be described. Note that this example does not limit the scope of the present invention.
【００４６】この実施例においては、４つの骨密度測定
用ファントム１０を用意し、その中心部１１の径を４ｍ
ｍφとし、それぞれのハイドロキシアパタイト濃度Ｄを
０，１００，２００および３００ｍｇ／ｃｍ3 とした。
ラット大腿骨骨幹部をハンドモータで４分割から６分割
程度に縦断たものをワックスを用いて中心部１１の周囲
に張り付け、これを外周部１２とした。これら４つの骨
密度測定用ファントムそれぞれについて、校正用投影デ
ータを収集して校正用ＣＴ画像を再構成し、中心部１１
の平均ＣＴ値および外周部１２の平均ＣＴ値を測定し
た。上記 (1)式に従って外周部１２の平均ＣＴ値で補正
して得られた中心部１１の平均ＣＴ値は、それぞれ２１
６，３４４，４７０および６００であった。In this embodiment, four bone density measuring phantoms 10 are prepared, and the diameter of the central portion 11 is 4 m.
mφ, and the respective hydroxyapatite concentrations D were 0, 100, 200 and 300 mg / cm 3 .
The femoral shaft of the rat femur was cut into four to six parts by a hand motor and attached around the center part 11 using wax. For each of these four bone density measurement phantoms, calibration projection data is collected, and a calibration CT image is reconstructed.
And the average CT value of the outer peripheral portion 12 were measured. The average CT value of the central portion 11 obtained by correcting the average CT value of the outer peripheral portion 12 according to the above equation (1) is 21
6, 344, 470 and 600.
【００４７】そして、補正された中心部１１の平均ＣＴ
値の検量線は、ハイドロキシアパタイト濃度Ｄの１次回
帰式として、Then, the corrected average CT of the central portion 11 is obtained.
The calibration curve of the value is a linear regression equation of the hydroxyapatite concentration D,
【数５】 が得られ、また、ハイドロキシアパタイト濃度Ｄの３次
回帰式として、(Equation 5) Is obtained, and as a cubic regression equation of the hydroxyapatite concentration D,
【数６】 が得られた。これら (5)式および (6)式から判るよう
に、２次および３次の項の寄与は非常に小さく、例え
ば、ハイドロキシアパタイト濃度Ｄの値が２００ｍｇ／
ｃｍ3 付近でＣＴ値は２％程度の差しか生じないことか
ら、上記 (5)式で表される検量線でも充分な精度を有す
るものとして用いられ得る。(Equation 6) was gotten. As can be seen from equations (5) and (6), the contribution of the second and third order terms is very small. For example, when the value of the hydroxyapatite concentration D is 200 mg /
Since the CT value is only about 2% around cm 3 , even the calibration curve represented by the above equation (5) can be used as having sufficient accuracy.
【００４８】被検体２０としては、１２週齢の雌性ラッ
トを１４匹用いた。第１の雌性ラット７匹については卵
巣摘除手術を行うことにより骨粗鬆症を発症させ、第２
の雌性ラット７匹については卵巣を摘除することなく偽
手術を行うことにより第１の雌性ラットの状態と近似さ
せた。これら１４匹の雌性ラットそれぞれについて、手
術から３週間経過後に大腿骨を摘出し、その大腿骨の遠
位部の被検体投影データを収集し被検体ＣＴ画像を再構
成して観察した。また、この被検体ＣＴ画像の皮質骨お
よび海綿骨それぞれのＣＴ値を測定し、上記 (3)式に従
って皮質骨のＣＴ値で補正して得られた海綿骨のＣＴ値
を上記 (5)式の検量線に当てはめて海綿骨の骨密度を求
めた。As the subject 20, 14 12-week-old female rats were used. Ovariectomy was performed on the first seven female rats to cause osteoporosis, and
Seven female rats were subjected to sham operation without removing the ovaries to approximate the condition of the first female rat. For each of these 14 female rats, the femur was excised three weeks after the operation, subject projection data of the distal part of the femur was collected, and a CT image of the subject was reconstructed and observed. Further, the CT value of cortical bone and trabecular bone of this subject CT image was measured, and the CT value of trabecular bone obtained by correcting the CT value of cortical bone according to the above equation (3) was used to calculate the CT value of the cancellous bone according to the above equation (5). The bone density of the cancellous bone was determined by applying the calibration curve of the above.
【００４９】その結果、卵巣を摘除した第１の雌性ラッ
トは、偽手術を行った第２の雌性ラットと比較して、骨
梁数の減少や骨梁間距離の増加が認められた。また、得
られた海綿骨の骨密度の平均値±標準偏差は、第１の雌
性ラットでは１４５．６±１６．６ｍｇ／ｃｍ3 であ
り、第２の雌性ラットでは２５４．１±５８．０ｍｇ／
ｃｍ3 であった。すなわち、卵巣摘除手術を行うことに
より明らかに海綿骨の骨密度の低下が認められた。この
ように雌性ラットを用いた閉経後骨粗鬆症モデルによ
り、骨の同一部位について内部微細構造を観察するとと
もに骨密度を測定することができることが確認できた。As a result, the number of trabecular bones and the distance between trabecular bones were observed in the first female rat from which the ovaries were removed, as compared with the second female rat in which the sham operation was performed. The average value ± standard deviation of the bone density of the obtained cancellous bone was 145.6 ± 16.6 mg / cm 3 in the first female rat, and 254.1 ± 58.0 mg in the second female rat. /
cm 3 . That is, the ovariectomy operation clearly reduced the bone density of cancellous bone. Thus, the postmenopausal osteoporosis model using female rats confirmed that it was possible to observe the internal microstructure and measure the bone density at the same site of the bone.
【００５０】[0050]
【発明の効果】以上、詳細に説明したとおり、本発明に
係る骨密度測定装置および方法によれば、既知であって
互いに異なる濃度のハイドロキシアパタイトを含有する
円柱形状の中心部と中心部の外周を皮質骨と等価な物質
で覆う外周部とを有する２以上の所定数の骨密度測定用
ファントムそれぞれについて、順次に校正用投影データ
が収集され、この校正用投影データに基づいて校正用Ｃ
Ｔ画像が再構成され、そして、この校正用ＣＴ画像に基
づいて、校正用ＣＴ画像の中心部のＣＴ値と中心部のハ
イドロキシアパタイト濃度との間の検量線が作成され
る。一方、被検体について被検体投影データが収集さ
れ、この被検体投影データに基づいて被検体ＣＴ画像が
再構成され、そして、この被検体ＣＴ画像の海綿骨のＣ
Ｔ値および検量線に基づいて、被検体の海綿骨の骨密度
が獲得される。このような構成としたので、被検体であ
る骨の海綿骨の微細構造および骨密度の双方を同時に、
非破壊的に且つ正確に測定することができる。したがっ
て、骨粗鬆症等の骨疾患の診断に好適に用いられ得る。As described above in detail, according to the apparatus and method for measuring bone density according to the present invention, the center and the outer periphery of the center of a columnar shape containing known and different concentrations of hydroxyapatite are used. The projection data for calibration is sequentially collected for each of two or more predetermined number of phantoms for measuring bone density having an outer peripheral portion which covers the cortical bone with a substance equivalent to cortical bone, and the calibration C data is obtained based on the projection data for calibration.
The T image is reconstructed, and a calibration curve between the CT value at the center of the calibration CT image and the hydroxyapatite concentration at the center is created based on the calibration CT image. On the other hand, subject projection data is collected for the subject, a subject CT image is reconstructed based on the subject projection data, and the cancellous bone C of the subject CT image is reconstructed.
Based on the T value and the calibration curve, the bone density of the cancellous bone of the subject is obtained. With such a configuration, both the microstructure and bone density of the cancellous bone of the subject bone,
It can be measured nondestructively and accurately. Therefore, it can be suitably used for diagnosis of bone diseases such as osteoporosis.
【００５１】また、所定数の骨密度測定用ファントムの
校正用ＣＴ画像それぞれの外周部のＣＴ値で補正した中
心部のＣＴ値と中心部のハイドロキシアパタイト濃度と
の間の検量線を作成するとともに、被検体ＣＴ画像の皮
質骨のＣＴ値で補正した海綿骨のＣＴ値および検量線に
基づいて被検体の海綿骨の骨密度を獲得する場合には、
骨密度の測定精度が高い。In addition, a calibration curve is created between the CT value at the center and the hydroxyapatite concentration at the center corrected with the CT value at the outer periphery of each of the calibration CT images of the predetermined number of phantoms for measuring bone density. When acquiring the bone density of the cancellous bone of the subject based on the CT value of the cancellous bone and the calibration curve corrected by the CT value of the cortical bone of the subject CT image,
High measurement accuracy of bone density.
【００５２】また、被検体のスライス面に垂直な方向に
一定間隔の各位置それぞれで被検体を測定し当該測定結
果を加算したものを被検体投影データとする場合や、被
検体のスライス面に垂直な方向の各位置それぞれで被検
体を測定し当該測定結果を補間し積分したものを被検体
投影データとする場合には、被検体の移動範囲に亘る平
均された測定結果が得られる。In addition, the object is measured at each position at regular intervals in the direction perpendicular to the slice plane of the object, and the result of the measurement is added as object projection data. When the object is measured at each position in the vertical direction, and the result of the measurement is interpolated and integrated to obtain object projection data, an averaged measurement result over the moving range of the object is obtained.
【００５３】また、本発明に係る骨密度測定用ファント
ムは、既知の濃度のハイドロキシアパタイトを含有する
円柱形状の中心部と、中心部の外周を皮質骨と等価な物
質で覆う外周部と、を備える構成としたので、上記の本
発明に係る骨密度測定用装置および方法による骨密度の
測定に好適に用いられる。Further, the phantom for measuring bone density according to the present invention comprises a cylindrical central portion containing a known concentration of hydroxyapatite, and an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone. Since it is provided with the structure, it is suitably used for the measurement of bone density by the above-described apparatus and method for measuring bone density according to the present invention.
【図１】本実施形態に係る骨密度測定用ファントムを説
明する斜視図である。FIG. 1 is a perspective view illustrating a bone density measurement phantom according to an embodiment.
【図２】本実施形態に係る骨密度測定装置の構成図であ
る。FIG. 2 is a configuration diagram of a bone density measurement device according to the present embodiment.
【図３】本実施形態に係る骨密度測定方法を説明するフ
ローチャート図である。FIG. 3 is a flowchart illustrating a bone density measurement method according to the embodiment.
【図４】本実施形態に係る骨密度測定方法に用いられる
複数個の骨密度測定用ファントムを説明する図表であ
る。FIG. 4 is a table illustrating a plurality of bone density measurement phantoms used in the bone density measurement method according to the embodiment.
【図５】本実施形態に係る骨密度測定方法における検量
線を説明するグラフである。FIG. 5 is a graph illustrating a calibration curve in the bone density measurement method according to the embodiment.
【図６】本実施形態に係る骨密度測定方法における検量
線に基づいて被検体の海綿骨の骨密度を獲得する方法を
説明するグラフである。FIG. 6 is a graph illustrating a method for obtaining the bone density of cancellous bone of a subject based on a calibration curve in the bone density measurement method according to the present embodiment.
【図７】本実施形態に係る骨密度測定方法における複数
スライス面の測定方法の第１の例を説明する図である。FIG. 7 is a diagram illustrating a first example of a method for measuring a plurality of slice planes in the bone density measurement method according to the embodiment.
【図８】本実施形態に係る骨密度測定方法における複数
スライス面の測定方法の第２の例を説明する図である。FIG. 8 is a diagram illustrating a second example of a method for measuring a plurality of slice planes in the bone density measurement method according to the embodiment.
【図９】本実施形態に係る骨密度測定方法における複数
スライス面の測定方法の第２の例を説明するグラフであ
る。FIG. 9 is a graph illustrating a second example of a method for measuring a plurality of slice planes in the bone density measurement method according to the embodiment.
１０…骨密度測定用ファントム、１１…中心部、１２…
外周部、２０…被検体、２１…海綿骨、２２…皮質骨、
１００…Ｘ線ＣＴ装置、１０１…測定室、１１０…Ｘ線
管、１２１…回転テーブル、１２２…Ｘ軸移動ステー
ジ、１２３…Ｙ軸移動ステージ、１２４…Ｚ軸移動ステ
ージ、１３０…リニアＸ線検出器、１３１…コリメー
タ、１３２…Ｘ軸移動ステージ、１３３…Ｙ軸移動ステ
ージ、１４０…Ｘ線イメージインテンシファイア、１４
３…Ｙ軸移動ステージ、１５０…架台、１５１…ステー
ジ移動レール、２００…制御演算装置、２１０…画像再
構成部、２２０…検量線作成部、２３０…被検体骨密度
獲得部、２４０…表示部、３１０…画像処理装置、３２
０…Ｘ線管制御装置。10 phantom for bone density measurement, 11 central part, 12…
Outer peripheral part, 20: subject, 21: cancellous bone, 22: cortical bone,
100 X-ray CT apparatus, 101 Measurement room, 110 X-ray tube, 121 Rotary table, 122 X-axis moving stage, 123 Y-axis moving stage, 124 Z-axis moving stage, 130 Linear X-ray detection 131, collimator, 132, X-axis moving stage, 133, Y-axis moving stage, 140, X-ray image intensifier, 14
3 Y-axis moving stage, 150 gantry, 151 stage moving rail, 200 control arithmetic unit, 210 image reconstructing unit, 220 calibration curve preparing unit, 230 subject bone density obtaining unit, 240 display unit .. 310 image processing device 32
0: X-ray tube controller.
───────────────────────────────────────────────────── フロントページの続き (72)発明者 菊池 一夫 大阪府堺市緑町４−152 株式会社日鉄 エレックス内 (72)発明者 平井 秀和 大阪府堺市緑町４−152 株式会社日鉄 エレックス内 (72)発明者 濱崎 貴生 大阪府堺市緑町４−152 株式会社日鉄 エレックス内 (56)参考文献 特開 平７−148143（ＪＰ，Ａ) 特開 平１−299533（ＪＰ，Ａ) 特開 平３−286743（ＪＰ，Ａ) 実開 平１−117315（ＪＰ，Ｕ) (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 6/00 - 6/14 ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Kazuo Kikuchi 4-152 Midoricho, Sakai-shi, Osaka Nippon Steel Elex Co., Ltd. (72) Inventor Hidekazu Hirai 4-152 Midoricho, Sakai-shi, Osaka Nippon Steel Elex Co., Ltd. 72) Inventor Takao Hamazaki 4-152 Midoricho, Sakai-shi, Osaka Nippon Steel Elex Co., Ltd. (56) Reference JP-A-7-148143 (JP, A) JP-A 1-299533 (JP, A) JP-A 3-286743 (JP, A) Japanese Utility Model 1-117315 (JP, U) (58) Field surveyed (Int. Cl. 7 , DB name) A61B 6/00-6/14
Claims (9)
出器アレイとを有し、前記測定空間内に置かれた被検体
に前記Ｘ線管によりＸ線を照射し、前記検出器アレイに
到達したＸ線の強度を測定して、投影データを収集する
投影データ収集手段と、 前記投影データに基づいて画像を再構成する画像再構成
手段と、 既知であって互いに異なる濃度のハイドロキシアパタイ
トを含有する円柱形状の中心部と前記中心部の外周を皮
質骨と等価な物質で覆う外周部とを有する２以上の所定
数の骨密度測定用ファントムそれぞれについて順次に前
記投影データ収集手段により収集された校正用投影デー
タに基づいて前記画像再構成手段により再構成された校
正用ＣＴ画像に基づいて、前記校正用ＣＴ画像の前記中
心部のＣＴ値と前記中心部のハイドロキシアパタイト濃
度との間の検量線を作成する検量線作成手段と、 前記被検体について前記投影データ収集手段により収集
された被検体投影データに基づいて前記画像再構成手段
により再構成された被検体ＣＴ画像の海綿骨のＣＴ値お
よび前記検量線に基づいて、前記被検体の海綿骨の骨密
度を獲得する被検体骨密度獲得手段と、 を備えることを特徴とする骨密度測定装置。1. An X-ray tube and a detector array arranged with a measurement space interposed therebetween, and an X-ray is radiated to a subject placed in the measurement space by the X-ray tube. Projection data collection means for measuring the intensity of X-rays reaching the array and collecting projection data; image reconstruction means for reconstructing an image based on the projection data; The projection data collecting means sequentially performs two or more predetermined numbers of bone density measurement phantoms each having a cylindrical central portion containing apatite and an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone. Based on the calibration CT image reconstructed by the image reconstruction means based on the acquired calibration projection data, the CT value of the central portion of the calibration CT image and the hydroxyl value of the central portion are obtained. Calibration curve creating means for creating a calibration curve between the apatite concentration; and a subject CT reconstructed by the image reconstruction means based on subject projection data collected by the projection data collecting means for the subject. A bone mineral density acquisition unit for acquiring the bone density of the trabecular bone of the subject based on the CT value of the trabecular bone of the image and the calibration curve.
密度測定用ファントムの校正用ＣＴ画像それぞれの前記
外周部のＣＴ値で補正した前記中心部のＣＴ値と前記中
心部のハイドロキシアパタイト濃度との間の検量線を作
成するとともに、 前記被検体骨密度獲得手段は、前記被検体ＣＴ画像の皮
質骨のＣＴ値で補正した海綿骨のＣＴ値および前記検量
線に基づいて前記被検体の海綿骨の骨密度を獲得する、 ことを特徴とする請求項１記載の骨密度測定装置。2. The method according to claim 1, wherein the calibration curve creating means includes: a CT value of the central portion corrected by a CT value of the outer peripheral portion of each of the calibration CT images of the predetermined number of bone density measurement phantoms; and a hydroxyapatite of the central portion. A calibration curve between the CT density and the CT curve of the trabecular bone corrected by the CT value of the cortical bone of the CT image of the subject and the calibration curve. The bone density measuring apparatus according to claim 1, wherein the bone density of cancellous bone is obtained.
のスライス面に垂直な方向に一定間隔の各位置それぞれ
で前記被検体を測定し、当該測定結果を加算したものを
前記被検体投影データとする、ことを特徴とする請求項
１記載の骨密度測定装置。3. The projection data collection means measures the subject at each position at regular intervals in a direction perpendicular to the slice plane of the subject, and adds the measurement result to the subject projection data. The bone density measuring device according to claim 1, wherein:
のスライス面に垂直な方向の各位置それぞれで前記被検
体を測定し、当該測定結果を補間し積分したものを前記
被検体投影データとする、ことを特徴とする請求項１記
載の骨密度測定装置。4. The projection data collection unit measures the subject at each position in a direction perpendicular to the slice plane of the subject, and interpolates and integrates the measurement result to obtain the subject projection data and The bone density measurement device according to claim 1, wherein
ロキシアパタイトを含有する円柱形状の中心部と前記中
心部の外周を皮質骨と等価な物質で覆う外周部とを有す
る２以上の所定数の骨密度測定用ファントムそれぞれを
順次に、Ｘ線ＣＴ装置の測定空間内の所定位置に配置
し、前記Ｘ線ＣＴ装置により測定して校正用投影データ
を収集する校正用投影データ収集ステップと、 前記所定数の骨密度測定用ファントムそれぞれについ
て、前記校正用投影データに基づいて校正用ＣＴ画像を
再構成する校正用ＣＴ画像再構成ステップと、 前記所定数の骨密度測定用ファントムの校正用ＣＴ画像
に基づいて前記中心部のＣＴ値と前記中心部のハイドロ
キシアパタイト濃度との間の検量線を作成する検量線作
成ステップと、 被検体である骨を前記所定位置に配置し、前記Ｘ線ＣＴ
装置により測定して被検体投影データを収集する被検体
投影データ収集ステップと、 前記被検体投影データに基づいて被検体ＣＴ画像を再構
成する被検体ＣＴ画像再構成ステップと、 前記被検体ＣＴ画像の海綿骨のＣＴ値および前記検量線
に基づいて前記被検体の海綿骨の骨密度を獲得する被検
体骨密度獲得ステップと、 を備えることを特徴とする骨密度測定方法。5. A predetermined number of at least two bones each having a cylindrical central part containing a known and different concentration of hydroxyapatite and an outer peripheral part covering the outer periphery of the central part with a substance equivalent to cortical bone. A calibration projection data collection step of sequentially arranging each of the density measurement phantoms at a predetermined position in a measurement space of the X-ray CT apparatus, and collecting calibration projection data by measuring with the X-ray CT apparatus; For each of the number of bone density measurement phantoms, a calibration CT image reconstruction step of reconstructing a calibration CT image based on the calibration projection data, and a calibration CT image of the predetermined number of bone density measurement phantoms. A calibration curve creating step of creating a calibration curve between the CT value of the central part and the hydroxyapatite concentration of the central part based on the bone; X-ray CT
An object projection data collection step of collecting object projection data by measuring with an apparatus; an object CT image reconstruction step of reconstructing an object CT image based on the object projection data; and the object CT image Obtaining a bone density of the cancellous bone of the subject based on the CT value of the cancellous bone of the subject and the calibration curve.
の骨密度測定用ファントムの校正用ＣＴ画像それぞれの
前記外周部のＣＴ値で補正した前記中心部のＣＴ値と前
記中心部のハイドロキシアパタイト濃度との間の検量線
を作成するとともに、 前記被検体骨密度獲得ステップは、前記被検体ＣＴ画像
の皮質骨のＣＴ値で補正した海綿骨のＣＴ値および前記
検量線に基づいて前記被検体の海綿骨の骨密度を獲得す
る、 ことを特徴とする請求項５記載の骨密度測定方法。6. The calibration curve creating step comprises the steps of: correcting the CT value of the central portion corrected by the CT value of the outer peripheral portion of each of the predetermined number of CT images for calibration of the bone density measurement phantom; and the hydroxyapatite of the central portion. And a calibration curve between the concentration and the concentration of the trabecular bone. The step of obtaining the bone density of the subject is performed based on the CT value of the cancellous bone corrected by the CT value of the cortical bone of the subject CT image and the calibration curve. The bone density measurement method according to claim 5, wherein the bone density of the cancellous bone is obtained.
前記被検体のスライス面に垂直な方向に一定間隔の各位
置それぞれで前記被検体を測定し、当該測定結果を加算
したものを前記被検体投影データとする、ことを特徴と
する請求項５記載の骨密度測定方法。7. The object projection data collecting step,
6. The object projection data, wherein the object is measured at respective positions at regular intervals in a direction perpendicular to a slice plane of the object, and a result obtained by adding the measurement results is used as the object projection data. Bone density measurement method.
前記被検体のスライス面に垂直な方向の各位置それぞれ
で前記被検体を測定し、当該測定結果を補間し積分した
ものを前記被検体投影データとする、ことを特徴とする
請求項５記載の骨密度測定方法。8. The subject projection data collection step includes:
6. The object projection data according to claim 5, wherein the object is measured at each position in a direction perpendicular to a slice plane of the object, and a result obtained by interpolating and integrating the measurement result is used as the object projection data. Bone density measurement method.
含有する円柱形状の中心部と、前記中心部の外周を皮質
骨と等価な物質で覆う外周部と、を備えることを特徴と
する骨密度測定用ファントム。9. A method for measuring bone density, comprising: a cylindrical central portion containing a known concentration of hydroxyapatite; and an outer peripheral portion covering the outer periphery of the central portion with a substance equivalent to cortical bone. phantom.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP32629297A JP3279513B2 (en) | 1997-11-27 | 1997-11-27 | Bone density measurement device, bone density measurement method, and bone density measurement phantom |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP32629297A JP3279513B2 (en) | 1997-11-27 | 1997-11-27 | Bone density measurement device, bone density measurement method, and bone density measurement phantom |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11155852A JPH11155852A (en) | 1999-06-15 |
| JP3279513B2 true JP3279513B2 (en) | 2002-04-30 |
Family
ID=18186150
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP32629297A Expired - Fee Related JP3279513B2 (en) | 1997-11-27 | 1997-11-27 | Bone density measurement device, bone density measurement method, and bone density measurement phantom |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3279513B2 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004184122A (en) * | 2002-11-29 | 2004-07-02 | Shimadzu Corp | Center axis calibrating fixture in x-ray ct equipment and calibration method using the fixture |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4748405B2 (en) * | 2001-04-18 | 2011-08-17 | 学校法人日本大学 | Material hardness measuring device |
| JP4831559B2 (en) * | 2004-03-30 | 2011-12-07 | 独立行政法人産業技術総合研究所 | Nanoapatite phantom and its use |
| US8139829B2 (en) | 2005-02-09 | 2012-03-20 | Hitachi Medical Corporation | Diagnostic imaging support system and diagnostic imaging support program |
| JP4573248B2 (en) * | 2006-09-28 | 2010-11-04 | 国立大学法人岐阜大学 | Phantom for ultrasonic bone measuring device and manufacturing method thereof |
| JP4998944B2 (en) * | 2007-06-20 | 2012-08-15 | 国立大学法人 筑波大学 | CT value-water equivalent thickness conversion apparatus and method, storage medium used therefor, and radiotherapy apparatus and method |
| KR101179540B1 (en) * | 2010-05-13 | 2012-09-04 | 국방과학연구소 | Method for analyzing the Internal Density of Material using a X-ray Computerized Tomography |
| JP6454626B2 (en) * | 2015-10-06 | 2019-01-16 | 株式会社リガク | Bone mineral density analyzer, analysis method, and analysis program |
| KR102227938B1 (en) * | 2018-06-20 | 2021-03-16 | 한국생산기술연구원 | Prediction method of osteoporosis |
-
1997
- 1997-11-27 JP JP32629297A patent/JP3279513B2/en not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004184122A (en) * | 2002-11-29 | 2004-07-02 | Shimadzu Corp | Center axis calibrating fixture in x-ray ct equipment and calibration method using the fixture |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH11155852A (en) | 1999-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jaffray et al. | Cone‐beam computed tomography with a flat‐panel imager: initial performance characterization | |
| US7457450B2 (en) | Method for recording and evaluating image data with the aid of a tomography machine | |
| US8290232B2 (en) | System and method for quantitative imaging of chemical composition to decompose more than two materials | |
| JP4152649B2 (en) | Method and apparatus for CT scout image processing | |
| US7885373B2 (en) | System and method for quantitative imaging of chemical composition to decompose multiple materials | |
| US8705822B2 (en) | Method for creating images indicating material decomposition in dual energy, dual source helical computed tomography | |
| JP6162241B2 (en) | X-ray CT system | |
| EP1149558A2 (en) | Apparatus for region of interest multislice CT scan | |
| WO2007074772A1 (en) | X-ray ct device | |
| EP2196147A1 (en) | X-ray computed tomography apparatus, medical image processing apparatus, X-ray computed tomography method, and medical image processing method | |
| US6493646B1 (en) | High order primary decay correction for CT imaging system detectors | |
| JP4429694B2 (en) | X-ray CT system | |
| JP3279513B2 (en) | Bone density measurement device, bone density measurement method, and bone density measurement phantom | |
| EP3738512A1 (en) | Systems and methods for automatic tube potential selection in dual energy imaging | |
| DE60132826T2 (en) | METHOD AND DEVICE FOR CALIBRATING X-RAY EFFECTORS IN A COMPUTER TOMOGRAPHIC PICTURE SYSTEM | |
| JP5022690B2 (en) | Radiography equipment | |
| JP4554185B2 (en) | X-ray CT system | |
| US6980681B1 (en) | Methods and apparatus for helical reconstruction for multislice CT scan | |
| US6904120B2 (en) | Method and apparatus for correcting bone induced spectral artifacts | |
| JP4887132B2 (en) | X-ray CT system | |
| US11517278B2 (en) | System and method for basis material decomposition with general physical constraint for multi-energy computed tomography | |
| JP3789728B2 (en) | Projection data correction method and apparatus, and radiation tomography apparatus | |
| JP2001212133A (en) | Synthetic tomographic image generating method | |
| JP2008017964A (en) | X-ray ct apparatus | |
| JP2000051205A (en) | Method for reconstituting image from measured value obtained by spiral scanning of examining object using ct device and ct device for executing this method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |