JPH04300525A - 骨塩定量分析方法 - Google Patents
骨塩定量分析方法Info
- Publication number
- JPH04300525A JPH04300525A JP3064289A JP6428991A JPH04300525A JP H04300525 A JPH04300525 A JP H04300525A JP 3064289 A JP3064289 A JP 3064289A JP 6428991 A JP6428991 A JP 6428991A JP H04300525 A JPH04300525 A JP H04300525A
- Authority
- JP
- Japan
- Prior art keywords
- bone
- bone mineral
- average value
- image
- image signal
- 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.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/20—Measuring radiation intensity with scintillation detectors
- G01T1/2012—Measuring radiation intensity with scintillation detectors using stimulable phosphors, e.g. stimulable phosphor sheets
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/505—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of bone
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Conversion Of X-Rays Into Visible Images (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【0001】
【産業上の利用分野】本発明は人体等の骨塩の定量方法
、さらに詳しくはエネルギーサブトラクションの手法を
用いて骨塩定量分析を行なう方法に関するものである。
、さらに詳しくはエネルギーサブトラクションの手法を
用いて骨塩定量分析を行なう方法に関するものである。
【0002】
【従来の技術】骨塩定量、すなわち骨の中のカルシウム
の量を定量的に測定することは骨折予防のために必要で
ある。すなわち、骨中のカルシウムの微量変化を知るこ
とは骨粗しょう症の早期発見を可能にし、骨折予防の効
果がある。
の量を定量的に測定することは骨折予防のために必要で
ある。すなわち、骨中のカルシウムの微量変化を知るこ
とは骨粗しょう症の早期発見を可能にし、骨折予防の効
果がある。
【0003】そこで従来、以下に列挙するように数々の
骨塩定量の方法が提案され、実施されている。
骨塩定量の方法が提案され、実施されている。
【0004】i) MD法(Microdensit
ometry) :◆これは、中指骨をアルミのステッ
プウエッジ(段階状パターン)とともにX線撮影し、濃
度計により濃度を測定し、アルミステップウエッジと対
応させてX線吸収量を換算し、さらに骨幅によりその値
を補正して骨塩を定量にするものであり、装置構成が簡
便であるが、定量の精度に問題がある上、骨粗しょう症
を最もよく表わす椎骨の測定ができないという欠点があ
る。
ometry) :◆これは、中指骨をアルミのステッ
プウエッジ(段階状パターン)とともにX線撮影し、濃
度計により濃度を測定し、アルミステップウエッジと対
応させてX線吸収量を換算し、さらに骨幅によりその値
を補正して骨塩を定量にするものであり、装置構成が簡
便であるが、定量の精度に問題がある上、骨粗しょう症
を最もよく表わす椎骨の測定ができないという欠点があ
る。
【0005】ii) SPA法(Single Pho
ton Absorptiometry):◆これは低
エネルギーγ線を骨に透過させた後15cm程離れたシ
ンチレーション検出器で検出し、γ線のカウント数の変
化によってアナログ計算から骨の単位長さ当りの重量を
求めるものであり、MD法に比べて正確な測定が可能で
あるが、これも椎骨の測定ができないという欠点がある
上、ラジオアイソトープを使用するため特別な管理が必
要であり、線源に半源期があるため線源を交換しなけれ
ばならないという難点がある。
ton Absorptiometry):◆これは低
エネルギーγ線を骨に透過させた後15cm程離れたシ
ンチレーション検出器で検出し、γ線のカウント数の変
化によってアナログ計算から骨の単位長さ当りの重量を
求めるものであり、MD法に比べて正確な測定が可能で
あるが、これも椎骨の測定ができないという欠点がある
上、ラジオアイソトープを使用するため特別な管理が必
要であり、線源に半源期があるため線源を交換しなけれ
ばならないという難点がある。
【0006】iii) DPA法(Dual Pho
ton Absorptiometry):◆これは4
4Ke Vと100 KeVの2種のエネルギーピーク
を有する核種である153 Gl を線源とし、この2
種のエネルギー線の骨の透過量の差異によって骨塩量を
測定するものであり、腰椎,大腿骨頚部の骨塩の測定、
および全身の骨塩量、脂肪量の高精度の測定が可能であ
るという利点があるが、これもラジオアイソトープを使
用することに伴う困難がある。また、放射線の照射が走
査方式であるため腰椎の場合10数分、全身では30〜
40分と検査時間がかかるという問題がある。
ton Absorptiometry):◆これは4
4Ke Vと100 KeVの2種のエネルギーピーク
を有する核種である153 Gl を線源とし、この2
種のエネルギー線の骨の透過量の差異によって骨塩量を
測定するものであり、腰椎,大腿骨頚部の骨塩の測定、
および全身の骨塩量、脂肪量の高精度の測定が可能であ
るという利点があるが、これもラジオアイソトープを使
用することに伴う困難がある。また、放射線の照射が走
査方式であるため腰椎の場合10数分、全身では30〜
40分と検査時間がかかるという問題がある。
【0007】iv) QDR法(Quantitati
ve Digited Radiography):◆ (別名DPX法)◆これはDPA法とほぼ同じであるが
、ラジオアイソトープの代りにパルス状X線をフィルタ
と組み合わせることにより2種類のエネルギーを得てお
り、再現性がよい上、検査時間を短縮(DPAの約1/
3 )する効果がある。簡便度,性能の両面から最も期
待されている方式であるが、検査時間は短縮されたとい
っても、腰椎撮影に約6分を要し、さらに短縮すること
が望まれている。
ve Digited Radiography):◆ (別名DPX法)◆これはDPA法とほぼ同じであるが
、ラジオアイソトープの代りにパルス状X線をフィルタ
と組み合わせることにより2種類のエネルギーを得てお
り、再現性がよい上、検査時間を短縮(DPAの約1/
3 )する効果がある。簡便度,性能の両面から最も期
待されている方式であるが、検査時間は短縮されたとい
っても、腰椎撮影に約6分を要し、さらに短縮すること
が望まれている。
【0008】v) QCT法(Quantitati
ve Computer Tomography):◆
これはX線CTを用いて、CTナンバーにより主に第3
腰椎の骨塩定量を行なうものであり、断面による定量化
が可能であるが、装置が大規模になってしまうという難
点がある。
ve Computer Tomography):◆
これはX線CTを用いて、CTナンバーにより主に第3
腰椎の骨塩定量を行なうものであり、断面による定量化
が可能であるが、装置が大規模になってしまうという難
点がある。
【0009】vi) DQCT法(Dual ener
gy Quantitative Computer
Tomography):◆これはQCT法において2
種のエネルギーを利用してエネルギーサブトラクション
を行なうことにより骨塩定量を行なうもので、骨組織内
の脂肪の影響を除いた定量化が可能であるという利点が
あるが、これも装置が大規模になってしまうという問題
がある。
gy Quantitative Computer
Tomography):◆これはQCT法において2
種のエネルギーを利用してエネルギーサブトラクション
を行なうことにより骨塩定量を行なうもので、骨組織内
の脂肪の影響を除いた定量化が可能であるという利点が
あるが、これも装置が大規模になってしまうという問題
がある。
【0010】以上列挙した通り、従来の骨塩定量の方法
は、簡便なものは精度が低く、高精度のものは装置が大
がかりになり、検査時間も長くなるといった問題があっ
た。
は、簡便なものは精度が低く、高精度のものは装置が大
がかりになり、検査時間も長くなるといった問題があっ
た。
【0011】そこで本出願人によりエネルギーサブトラ
クションを用いた骨塩定量分析方法が提案されている(
特願平2−114,611 号参照)。このエネルギー
サブトラクションを用いた方法とは、2枚以上の蓄積性
蛍光体シート、X線フイルム等の記録シートのそれぞれ
に、軟部組織と骨部組織を含む被写体を透過したそれぞ
れエネルギーが異なる放射線を照射して、これらの記録
シートに前記被写体の放射線画像を記録し、これらの記
録シートに励起光を走査して前記放射線画像を光電的に
読み取ってデジタル画像信号に変換し、各画像の対応す
る画素間でこのデジタル画像信号の減算を行なって放射
線画像の前記骨部組織のみの画像を形成する差信号を得
るエネルギーサブトラクションにおいて、前記被写体の
放射線画像を得る際に骨塩量が段階的に変化した人骨を
模擬した骨塩レファレンスを同時に写し込んでおき、前
記骨部組織のみの画像(骨部画像)上で骨部組織の陰影
の濃度と骨塩レファレンスの濃度とを比較することによ
り骨塩量を定量化する方法である。
クションを用いた骨塩定量分析方法が提案されている(
特願平2−114,611 号参照)。このエネルギー
サブトラクションを用いた方法とは、2枚以上の蓄積性
蛍光体シート、X線フイルム等の記録シートのそれぞれ
に、軟部組織と骨部組織を含む被写体を透過したそれぞ
れエネルギーが異なる放射線を照射して、これらの記録
シートに前記被写体の放射線画像を記録し、これらの記
録シートに励起光を走査して前記放射線画像を光電的に
読み取ってデジタル画像信号に変換し、各画像の対応す
る画素間でこのデジタル画像信号の減算を行なって放射
線画像の前記骨部組織のみの画像を形成する差信号を得
るエネルギーサブトラクションにおいて、前記被写体の
放射線画像を得る際に骨塩量が段階的に変化した人骨を
模擬した骨塩レファレンスを同時に写し込んでおき、前
記骨部組織のみの画像(骨部画像)上で骨部組織の陰影
の濃度と骨塩レファレンスの濃度とを比較することによ
り骨塩量を定量化する方法である。
【0012】
【発明が解決しようとする課題】しかし、上記エネルギ
ーサブトラクションを用いる方法は簡便かつそれ以前の
方法と比べかなり高精度に骨塩量を定量化することがで
きるが、さらに高精度化を図るには放射線撮影の際に該
放射線が被写体の各組織により散乱されること等に起因
し同一の骨塩量の骨部組織にあっても放射線画像上の位
置によりその濃度が多少異なり、したがって骨部組織の
陰影と骨塩レファレンスの陰影の放射線画像上の位置の
相違により測定誤差を含むという問題がある。
ーサブトラクションを用いる方法は簡便かつそれ以前の
方法と比べかなり高精度に骨塩量を定量化することがで
きるが、さらに高精度化を図るには放射線撮影の際に該
放射線が被写体の各組織により散乱されること等に起因
し同一の骨塩量の骨部組織にあっても放射線画像上の位
置によりその濃度が多少異なり、したがって骨部組織の
陰影と骨塩レファレンスの陰影の放射線画像上の位置の
相違により測定誤差を含むという問題がある。
【0013】本発明は上記問題に鑑み、散乱線等の影響
を少なくすることによりさらに高精度化を図ったエネル
ギーサブトラクションを用いた骨塩定量分析方法を提供
することを目的とするものである。
を少なくすることによりさらに高精度化を図ったエネル
ギーサブトラクションを用いた骨塩定量分析方法を提供
することを目的とするものである。
【0014】
【課題を解決するための手段】上記目的を達成するため
の本発明の骨塩定量分析方法は、軟部組織と骨部組織と
からなる被検体と骨塩量を段階的に模擬した骨塩レファ
レンスとを被写体とする、互いに異なるエネルギーの放
射線により形成された複数の放射線画像に基づいて、前
記被検体中の骨部組織の陰影が抽出もしくは強調された
骨部画像を生成し、該骨部画像の前記骨塩レファレンス
の部分の画像信号の平均値と、前記骨部組織および前記
骨塩レファレンスの存在しないバックグラウンド領域の
画像信号であるバックグラウンド信号のうち、前記骨塩
レファレンスの周辺のバックグラウンド信号の平均値と
の差分値を、前記骨塩レファレンスの骨塩量の段階毎に
求め、該差分値と、所望とする前記骨部組識の画像信号
の平均値と該所望とする骨部組織の周辺のバックグラウ
ンド信号の平均値との差分値とを比較することにより前
記所望とする骨部組識の骨塩量を求めることを特徴とす
るものである。
の本発明の骨塩定量分析方法は、軟部組織と骨部組織と
からなる被検体と骨塩量を段階的に模擬した骨塩レファ
レンスとを被写体とする、互いに異なるエネルギーの放
射線により形成された複数の放射線画像に基づいて、前
記被検体中の骨部組織の陰影が抽出もしくは強調された
骨部画像を生成し、該骨部画像の前記骨塩レファレンス
の部分の画像信号の平均値と、前記骨部組織および前記
骨塩レファレンスの存在しないバックグラウンド領域の
画像信号であるバックグラウンド信号のうち、前記骨塩
レファレンスの周辺のバックグラウンド信号の平均値と
の差分値を、前記骨塩レファレンスの骨塩量の段階毎に
求め、該差分値と、所望とする前記骨部組識の画像信号
の平均値と該所望とする骨部組織の周辺のバックグラウ
ンド信号の平均値との差分値とを比較することにより前
記所望とする骨部組識の骨塩量を求めることを特徴とす
るものである。
【0015】
【作用】本発明の骨塩定量分析方法は、骨部画像の骨塩
レファレンスの部分の画像信号の平均値と、骨塩レファ
レンスの周辺のバックグラウンド信号の平均値との差分
値を骨塩レファレンスの骨塩量の段階毎に求め、この差
分値と、所望とする骨部組織の画像信号の平均値とその
周辺のバックグラウンド信号の平均値との差分値とを比
較することにより所望とする骨部組織の骨塩量を求める
ようにしたため、散乱放射線等を含むバックグラウンド
信号の影響を受けることなく高精度で骨塩量を求めるこ
とが可能となる。
レファレンスの部分の画像信号の平均値と、骨塩レファ
レンスの周辺のバックグラウンド信号の平均値との差分
値を骨塩レファレンスの骨塩量の段階毎に求め、この差
分値と、所望とする骨部組織の画像信号の平均値とその
周辺のバックグラウンド信号の平均値との差分値とを比
較することにより所望とする骨部組織の骨塩量を求める
ようにしたため、散乱放射線等を含むバックグラウンド
信号の影響を受けることなく高精度で骨塩量を求めるこ
とが可能となる。
【0016】
【実施例】以下、図面に示す実施例に基づいて本発明を
詳細に説明する。尚、ここでは蓄積性蛍光体シートを用
いた例について説明する。
詳細に説明する。尚、ここでは蓄積性蛍光体シートを用
いた例について説明する。
【0017】図1(a) は2枚の蓄積性蛍光体シート
A、Bに、軟部組織と骨とを有する同一の被写体1を透
過したX線2を、それぞれエネルギーを変えて照射する
状態を示す。すなわち第1の蓄積性蛍光体シートAに被
写体1のΧ線透過像を蓄積記録し、次いで短時間内で蓄
積性蛍光体シートA、Bを素早く取り替えると同時に、
Χ線源3の管電圧を変えて、透過Χ線のエネルギーが異
なる被写体1のΧ線画像を蓄積性蛍光体シートBに蓄積
記録する。このとき蓄積性蛍光体シートAとBとで被写
体1の位置関係は同じとする。
A、Bに、軟部組織と骨とを有する同一の被写体1を透
過したX線2を、それぞれエネルギーを変えて照射する
状態を示す。すなわち第1の蓄積性蛍光体シートAに被
写体1のΧ線透過像を蓄積記録し、次いで短時間内で蓄
積性蛍光体シートA、Bを素早く取り替えると同時に、
Χ線源3の管電圧を変えて、透過Χ線のエネルギーが異
なる被写体1のΧ線画像を蓄積性蛍光体シートBに蓄積
記録する。このとき蓄積性蛍光体シートAとBとで被写
体1の位置関係は同じとする。
【0018】このとき、被写体1とともに、段階的にX
線吸収量が異なるパターンからなるX線吸収量の既知の
骨塩レファレンス(ファントム)5を蓄積性蛍光体シー
トA,B上に置いて、被写体1のX線画像とともにこの
ファントム5のX線画像も蓄積性蛍光体シートA、Bに
蓄積記録する。このファントム5は、図2に示すように
、骨塩の量すなわちCa CO3 の含有量(wt%)
が段階的に異なるセクション5a,5b ……5fを並
べた構造をしており、このCa CO3 の含有量は予
め知られているものである。
線吸収量が異なるパターンからなるX線吸収量の既知の
骨塩レファレンス(ファントム)5を蓄積性蛍光体シー
トA,B上に置いて、被写体1のX線画像とともにこの
ファントム5のX線画像も蓄積性蛍光体シートA、Bに
蓄積記録する。このファントム5は、図2に示すように
、骨塩の量すなわちCa CO3 の含有量(wt%)
が段階的に異なるセクション5a,5b ……5fを並
べた構造をしており、このCa CO3 の含有量は予
め知られているものである。
【0019】また、図1(b) は2枚の蓄積性蛍光体
シートA、Bを重ね、この間に放射線エネルギーを一部
吸収するフィルタFを介在させて被写体1とファントム
5を透過したX線2を照射する状態を示すもので、これ
によりエネルギーの大きさの異なる放射線を蓄積性蛍光
体シートA、Bに同時に照射するもの(いわゆるワンシ
ョットエネルギーサブトラクション)である。ワンショ
ットエネルギーサブトラクションについては特開昭59
−83486号に詳細が開示されている。
シートA、Bを重ね、この間に放射線エネルギーを一部
吸収するフィルタFを介在させて被写体1とファントム
5を透過したX線2を照射する状態を示すもので、これ
によりエネルギーの大きさの異なる放射線を蓄積性蛍光
体シートA、Bに同時に照射するもの(いわゆるワンシ
ョットエネルギーサブトラクション)である。ワンショ
ットエネルギーサブトラクションについては特開昭59
−83486号に詳細が開示されている。
【0020】上記図1(a) もしくは図1(b) の
いずれかの方法を用いて、2つの放射線画像を2枚の蓄
積性蛍光体シートA、Bに蓄積記録する。次にこれら2
枚の蓄積性蛍光体シートA、Bから、図3に示すような
画像読取手段によってX線画像を読み取り、画像を表わ
すデジタル画像信号を得る。先ず、蓄積性蛍光体シート
Aを矢印Yの方向に副走査のために移動させながら、レ
ーザー光源10からのレーザー光11を走査ミラー12
によってX方向に主走査させ、蛍光体シートAから蓄積
X線エネルギーを、蓄積記録されたX線画像にしたがっ
て輝尽発光光13として発散させる。輝尽発光光13は
透明なアクリル板を成形して作られた集光板14の一端
面からこの集光板14の内部に入射し、中を全反射を繰
返しつつフォトマル15に至り、輝尽発光光13の発光
量が画像信号Sとして出力される。この出力された画像
信号Sは増幅器とA/D変換器を含む対数変換器16に
より対数値(log S)のデジタル画像信号log
SA に変換される。このデジタル画像信号log S
A は例えば磁気ディスク等の記憶媒体17に記憶され
る。次に、全く同様にして、もう1枚の蓄積性蛍光体シ
ートBの記録画像が読み出され、そのデジタル画像信号
log SB が同様に記憶媒体17に記憶される。
いずれかの方法を用いて、2つの放射線画像を2枚の蓄
積性蛍光体シートA、Bに蓄積記録する。次にこれら2
枚の蓄積性蛍光体シートA、Bから、図3に示すような
画像読取手段によってX線画像を読み取り、画像を表わ
すデジタル画像信号を得る。先ず、蓄積性蛍光体シート
Aを矢印Yの方向に副走査のために移動させながら、レ
ーザー光源10からのレーザー光11を走査ミラー12
によってX方向に主走査させ、蛍光体シートAから蓄積
X線エネルギーを、蓄積記録されたX線画像にしたがっ
て輝尽発光光13として発散させる。輝尽発光光13は
透明なアクリル板を成形して作られた集光板14の一端
面からこの集光板14の内部に入射し、中を全反射を繰
返しつつフォトマル15に至り、輝尽発光光13の発光
量が画像信号Sとして出力される。この出力された画像
信号Sは増幅器とA/D変換器を含む対数変換器16に
より対数値(log S)のデジタル画像信号log
SA に変換される。このデジタル画像信号log S
A は例えば磁気ディスク等の記憶媒体17に記憶され
る。次に、全く同様にして、もう1枚の蓄積性蛍光体シ
ートBの記録画像が読み出され、そのデジタル画像信号
log SB が同様に記憶媒体17に記憶される。
【0021】図4は上記のようにして得られた2つのデ
ィジタル画像信号log SA 、log SB に基
づくサブトラクション処理を模式的に表わした図である
。画像18は、画像信号log SA が担持する画像
であって撮影に際しX線源3の管電圧を低圧(例えば6
0kV) として撮影により得られた画像である。また
画像19はX線源の管電圧を高圧(120kV)とした
撮影により得られた画像である。これらの2つの画像1
8、19のそれぞれには、人体の腰部の骨の陰影18a
、19a 、軟部組織の陰影18b 、19b、骨塩
レファレンスの陰影18c 、19c が写し込まれて
いる。
ィジタル画像信号log SA 、log SB に基
づくサブトラクション処理を模式的に表わした図である
。画像18は、画像信号log SA が担持する画像
であって撮影に際しX線源3の管電圧を低圧(例えば6
0kV) として撮影により得られた画像である。また
画像19はX線源の管電圧を高圧(120kV)とした
撮影により得られた画像である。これらの2つの画像1
8、19のそれぞれには、人体の腰部の骨の陰影18a
、19a 、軟部組織の陰影18b 、19b、骨塩
レファレンスの陰影18c 、19c が写し込まれて
いる。
【0022】ここで、軟部組織18b 、19b は各
個人によりその厚さ等が大幅に異なりこのままでは骨塩
量の定量化が難しいため、サブトラクション処理、即ち
2つの画像18、19の互いに対応する各画素毎に◆
log S=A logSA
−B logSB +C ……
(1) ◆(但し、A、B、Cは係数を表わす)◆の演
算処理を行なうことにより軟部組織の陰影18b 、1
9b が消去された骨部画像20を得る。
個人によりその厚さ等が大幅に異なりこのままでは骨塩
量の定量化が難しいため、サブトラクション処理、即ち
2つの画像18、19の互いに対応する各画素毎に◆
log S=A logSA
−B logSB +C ……
(1) ◆(但し、A、B、Cは係数を表わす)◆の演
算処理を行なうことにより軟部組織の陰影18b 、1
9b が消去された骨部画像20を得る。
【0023】図5は上記のようにして求められた骨部画
像20より、骨塩レファレンス20c の部分Cの画像
信号の平均値SC と骨塩レファレンス20c の周辺
部Dのバックグラウンド信号の平均値SD との差分値
および骨部20a の部分C′の画像信号の平均値SC
′と骨部20a の周辺部D′のバックグラウンド信
号の平均値SD ′との差分値の求め方を示した図であ
る。
像20より、骨塩レファレンス20c の部分Cの画像
信号の平均値SC と骨塩レファレンス20c の周辺
部Dのバックグラウンド信号の平均値SD との差分値
および骨部20a の部分C′の画像信号の平均値SC
′と骨部20a の周辺部D′のバックグラウンド信
号の平均値SD ′との差分値の求め方を示した図であ
る。
【0024】まず図5(a) に示す骨塩レファレンス
20c のセクション5aの部分Cの画像信号の平均値
SC と、骨塩レファレンス20c のセクション5a
の周辺部Dのバックグラウンド信号の平均値SD とを
求め、この画像信号の平均値SC とバックグラウンド
信号の平均値SD との差分値SC −SD を求める
。この差分値を求めるということは骨塩レファレンス2
0c に重畳されているバックグラウンド信号から、こ
のバックグラウンド信号と略等しい骨塩レファレンス2
0c の周辺部Dのバックグラウンド信号SD を除去
することであるため、差分値SC −SD は散乱放射
線の影響を受けているバックグラウンド信号に依存しな
い精度の高い画像信号となる。
20c のセクション5aの部分Cの画像信号の平均値
SC と、骨塩レファレンス20c のセクション5a
の周辺部Dのバックグラウンド信号の平均値SD とを
求め、この画像信号の平均値SC とバックグラウンド
信号の平均値SD との差分値SC −SD を求める
。この差分値を求めるということは骨塩レファレンス2
0c に重畳されているバックグラウンド信号から、こ
のバックグラウンド信号と略等しい骨塩レファレンス2
0c の周辺部Dのバックグラウンド信号SD を除去
することであるため、差分値SC −SD は散乱放射
線の影響を受けているバックグラウンド信号に依存しな
い精度の高い画像信号となる。
【0025】上記のようにして骨塩レファレンス20c
のセクション5aから5fの画像信号の平均値と、周
辺部のバックグラウンド信号の平均値との差分値を求め
る。
のセクション5aから5fの画像信号の平均値と、周
辺部のバックグラウンド信号の平均値との差分値を求め
る。
【0026】また、同様にして所望とする骨部20a
の部分C′の画像信号の平均値SC ′と骨部20a
の周辺部D′のバックグラウンド信号の平均値SD ′
との差分値SC ′−SD ′を求める。
の部分C′の画像信号の平均値SC ′と骨部20a
の周辺部D′のバックグラウンド信号の平均値SD ′
との差分値SC ′−SD ′を求める。
【0027】図6は上記のようにして求められた骨塩レ
ファレンス20c の差分値に対する骨塩量を表わす直
線を描いた図である。
ファレンス20c の差分値に対する骨塩量を表わす直
線を描いた図である。
【0028】この直線24により、例えば骨部画像20
の所望とする骨部20a の画像信号の平均値とその周
辺部のバックグラウンド信号の平均値との差分値がΔS
1 の場合、その骨部の骨塩量はV1と求められる。
の所望とする骨部20a の画像信号の平均値とその周
辺部のバックグラウンド信号の平均値との差分値がΔS
1 の場合、その骨部の骨塩量はV1と求められる。
【0029】なお、上記実施例は蓄積性蛍光体シートを
用いた例であるが、蓄積性蛍光体シートに代え、X線フ
イルム等を用いてもよいことはもちろんである。
用いた例であるが、蓄積性蛍光体シートに代え、X線フ
イルム等を用いてもよいことはもちろんである。
【0030】また上記実施例においては段階的にX線吸
収量が異なるパターンが一体となっている骨塩レファレ
ンス(ファントム)5を用いたが、図7に示すようにX
線吸収量が異なる段階毎に互いに間隔を開けて配置させ
るようにしてもよい。
収量が異なるパターンが一体となっている骨塩レファレ
ンス(ファントム)5を用いたが、図7に示すようにX
線吸収量が異なる段階毎に互いに間隔を開けて配置させ
るようにしてもよい。
【0031】さらに上記実施例において骨塩レファレン
ス20c および骨部20a の周辺部のバックグラウ
ンド信号の平均値SD およびSD ′を求める際に、
骨塩レファレンス20c および骨部20a の周辺の
一部分より求めるようにしているが、骨塩レファレンス
20c および骨部20a の周辺であればいかなる部
分でもよく、例えば図8に示すように骨塩レファレンス
20c および骨部20a の周囲を囲む部分よりバッ
クグラウンド信号の平均値SD およびSD ′を求め
るようにしてもよい。
ス20c および骨部20a の周辺部のバックグラウ
ンド信号の平均値SD およびSD ′を求める際に、
骨塩レファレンス20c および骨部20a の周辺の
一部分より求めるようにしているが、骨塩レファレンス
20c および骨部20a の周辺であればいかなる部
分でもよく、例えば図8に示すように骨塩レファレンス
20c および骨部20a の周囲を囲む部分よりバッ
クグラウンド信号の平均値SD およびSD ′を求め
るようにしてもよい。
【0032】また上記校正線24はこの実施例では直線
で近似してあるが、曲線であってもよいのはもちろんで
ある。
で近似してあるが、曲線であってもよいのはもちろんで
ある。
【0033】
【発明の効果】以上詳細に説明したように、本発明の骨
塩定量分析方法は、骨部画像の骨塩レファレンスの部分
の画像信号の平均値と、骨塩レファレンスの周辺のバッ
クグラウンド信号の平均値との差分値を骨塩レファレン
スの骨塩量の段階毎に求め、この差分値と、所望とする
骨部組織の画像信号の平均値と、その周辺のバックグラ
ウンド信号の平均値との差分値とを比較して所望とする
骨部組織の骨塩量を求めるようにしたため、散乱放射線
等を含むバックグラウンド信号の影響を除去することが
でき、散乱放射線等に起因する測定精度の低下が改善さ
れ、高精度に骨塩量を求めることができる。
塩定量分析方法は、骨部画像の骨塩レファレンスの部分
の画像信号の平均値と、骨塩レファレンスの周辺のバッ
クグラウンド信号の平均値との差分値を骨塩レファレン
スの骨塩量の段階毎に求め、この差分値と、所望とする
骨部組織の画像信号の平均値と、その周辺のバックグラ
ウンド信号の平均値との差分値とを比較して所望とする
骨部組織の骨塩量を求めるようにしたため、散乱放射線
等を含むバックグラウンド信号の影響を除去することが
でき、散乱放射線等に起因する測定精度の低下が改善さ
れ、高精度に骨塩量を求めることができる。
【図1】本発明の骨塩定量分析方法の一実施例を実行す
るための撮影ステップを示す図
るための撮影ステップを示す図
【図2】上記撮影に使用する骨塩レファレンスを示す図
【図3】本発明の一例を実行するための読取ステップを
示す図
示す図
【図4】本発明の一例を実行する際のサブトラクション
処理を模式的に表わした図
処理を模式的に表わした図
【図5】骨塩レファレンスの画像信号の平均値と周辺の
バックグラウンド信号の平均値との差分値および骨部の
画像信号の平均値と周辺のバックグラウンド信号の平均
値との差分値の求め方を示した図
バックグラウンド信号の平均値との差分値および骨部の
画像信号の平均値と周辺のバックグラウンド信号の平均
値との差分値の求め方を示した図
【図6】差分値に対する骨塩量を表わす直線を描いた図
【図7】本発明の実施例による骨塩レファレンスの別の
配置の方法を示した図
配置の方法を示した図
【図8】骨塩レファレンスの画像信号の平均値と周辺の
バックグラウンド信号の平均値との差分値および骨部の
画像信号の平均値と周辺のバックグラウンド信号の平均
値との差分値の別の求め方を示した図
バックグラウンド信号の平均値との差分値および骨部の
画像信号の平均値と周辺のバックグラウンド信号の平均
値との差分値の別の求め方を示した図
1 被写体
2 X線
3 X線源
4 ファントム(骨塩レファレンス)5
骨塩ブロック 10 レーザ光源 12 走査ミラー 13 輝尽発光光 15 フォトマル 20 骨部画像
骨塩ブロック 10 レーザ光源 12 走査ミラー 13 輝尽発光光 15 フォトマル 20 骨部画像
Claims (1)
- 【請求項1】 軟部組織と骨部組織とからなる被検体
と骨塩量を段階的に模擬した骨塩レファレンスとを被写
体とする、互いに異なるエネルギーの放射線により形成
された複数の放射線画像に基づいて、前記被検体中の骨
部組織の陰影が抽出もしくは強調された骨部画像を生成
し、該骨部画像の前記骨塩レファレンスの部分の画像信
号の平均値と、前記骨部組織および前記骨塩レファレン
スの存在しないバックグラウンド領域の画像信号である
バックグラウンド信号のうち、前記骨塩レファレンスの
周辺のバックグラウンド信号の平均値との差分値を、前
記骨塩レファレンスの骨塩量の段階毎に求め、該差分値
と、所望とする前記骨部組識の画像信号の平均値と該所
望とする骨部組織の周辺のバックグラウンド信号の平均
値との差分値とを比較することにより前記所望とする骨
部組識の骨塩量を求めることを特徴とする骨塩定量分析
方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3064289A JP2641078B2 (ja) | 1991-03-28 | 1991-03-28 | 骨塩定量分析方法 |
| US07/858,634 US5187731A (en) | 1991-03-28 | 1992-03-27 | Method for quantitatively analyzing bone calcium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3064289A JP2641078B2 (ja) | 1991-03-28 | 1991-03-28 | 骨塩定量分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04300525A true JPH04300525A (ja) | 1992-10-23 |
| JP2641078B2 JP2641078B2 (ja) | 1997-08-13 |
Family
ID=13253927
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3064289A Expired - Fee Related JP2641078B2 (ja) | 1991-03-28 | 1991-03-28 | 骨塩定量分析方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5187731A (ja) |
| JP (1) | JP2641078B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014054301A (ja) * | 2012-09-11 | 2014-03-27 | Hitachi Aloka Medical Ltd | 医療用x線測定装置 |
Families Citing this family (39)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2981701B2 (ja) * | 1992-09-08 | 1999-11-22 | 富士写真フイルム株式会社 | 放射線画像形成方法およびそれを用いた骨塩定量分析方法 |
| US5335260A (en) * | 1992-11-25 | 1994-08-02 | Arnold Ben A | Calibration phantom and improved method of quantifying calcium and bone density using same |
| US5420441A (en) * | 1993-11-23 | 1995-05-30 | Eastman Kodak Company | Automated technique for calibrating a storage phosphor reader |
| US5481587A (en) * | 1994-05-09 | 1996-01-02 | Lunar Corporation | Radiographic phantom for vertebral morphometry |
| WO1996007161A1 (en) * | 1994-08-29 | 1996-03-07 | Eskofot A/S | A method of estimation |
| US5506884A (en) * | 1995-08-25 | 1996-04-09 | Phantom Laboratory, Incorporated | Radiation phantom and test methods employing the same |
| US6064716A (en) * | 1997-09-05 | 2000-05-16 | Cyberlogic, Inc. | Plain x-ray bone densitometry apparatus and method |
| US5917877A (en) * | 1997-09-05 | 1999-06-29 | Cyberlogic, Inc. | Plain x-ray bone densitometry apparatus and method |
| US5852647A (en) * | 1997-09-24 | 1998-12-22 | Schick Technologies | Method and apparatus for measuring bone density |
| US6320931B1 (en) | 1998-03-02 | 2001-11-20 | Image Analysis, Inc. | Automated x-ray bone densitometer |
| US6173038B1 (en) | 1999-12-01 | 2001-01-09 | Cyberlogic, Inc. | Plain x-ray bone densitometry apparatus and method |
| US6510197B1 (en) | 2000-01-11 | 2003-01-21 | Alara, Inc. | Method and apparatus for osteoporosis screening |
| US7467892B2 (en) | 2000-08-29 | 2008-12-23 | Imaging Therapeutics, Inc. | Calibration devices and methods of use thereof |
| US20020186818A1 (en) * | 2000-08-29 | 2002-12-12 | Osteonet, Inc. | System and method for building and manipulating a centralized measurement value database |
| US6904123B2 (en) | 2000-08-29 | 2005-06-07 | Imaging Therapeutics, Inc. | Methods and devices for quantitative analysis of x-ray images |
| US7050534B2 (en) * | 2000-08-29 | 2006-05-23 | Imaging Therapeutics, Inc. | Methods and devices for quantitative analysis of x-ray images |
| CN1498385A (zh) * | 2000-08-29 | 2004-05-19 | X图像量析的方法与装置 | |
| US7660453B2 (en) | 2000-10-11 | 2010-02-09 | Imaging Therapeutics, Inc. | Methods and devices for analysis of x-ray images |
| US20070047794A1 (en) * | 2000-10-11 | 2007-03-01 | Philipp Lang | Methods and devices for analysis of x-ray images |
| US8639009B2 (en) * | 2000-10-11 | 2014-01-28 | Imatx, Inc. | Methods and devices for evaluating and treating a bone condition based on x-ray image analysis |
| EP1389947B1 (en) * | 2001-05-25 | 2009-08-26 | Imaging Therapeutics, Inc. | Methods to diagnose treat and prevent bone loss |
| US6990222B2 (en) * | 2001-11-21 | 2006-01-24 | Arnold Ben A | Calibration of tissue densities in computerized tomography |
| US6570955B1 (en) | 2002-01-08 | 2003-05-27 | Cyberlogic, Inc. | Digital x-ray material testing and bone densitometry apparatus and method |
| DE10219751A1 (de) * | 2002-05-02 | 2003-11-13 | Siemens Ag | Röntgendetektor mit einer Konverterschicht |
| EP1546982A1 (en) * | 2002-09-16 | 2005-06-29 | Imaging Therapeutics, Inc. | Imaging markers in musculoskeletal disease |
| US7840247B2 (en) | 2002-09-16 | 2010-11-23 | Imatx, Inc. | Methods of predicting musculoskeletal disease |
| US8965075B2 (en) | 2002-09-16 | 2015-02-24 | Imatx, Inc. | System and method for predicting future fractures |
| EP1605824A2 (en) | 2003-03-25 | 2005-12-21 | Imaging Therapeutics, Inc. | Methods for the compensation of imaging technique in the processing of radiographic images |
| US7488109B2 (en) * | 2003-03-27 | 2009-02-10 | Wright State University | Osteoporosis screening using radiographic absorptiometry of the mandible |
| US8073521B2 (en) | 2003-09-19 | 2011-12-06 | Imatx, Inc. | Method for bone structure prognosis and simulated bone remodeling |
| US8290564B2 (en) * | 2003-09-19 | 2012-10-16 | Imatx, Inc. | Method for bone structure prognosis and simulated bone remodeling |
| GB0420224D0 (en) * | 2004-09-11 | 2004-10-13 | Koninkl Philips Electronics Nv | Coherent scatter imaging |
| GB0420222D0 (en) * | 2004-09-11 | 2004-10-13 | Koninkl Philips Electronics Nv | Coherent scatter imaging |
| CA2580726A1 (en) | 2004-09-16 | 2006-03-30 | Imaging Therapeutics, Inc. | System and method of predicting future fractures |
| US7657000B2 (en) * | 2006-11-27 | 2010-02-02 | Cmt Medical Technologies Ltd. | Method and apparatus for dual energy radiography |
| US7959742B2 (en) * | 2007-07-11 | 2011-06-14 | Whirlpool Corporation | Outer support body for a drawer-type dishwasher |
| DE102009008793A1 (de) * | 2009-02-13 | 2010-11-11 | Inselspital-Stiftung | Verfahren zum Erzeugen einer Schichtbildaufnahme |
| US8939917B2 (en) * | 2009-02-13 | 2015-01-27 | Imatx, Inc. | Methods and devices for quantitative analysis of bone and cartilage |
| CN112914505A (zh) * | 2020-03-09 | 2021-06-08 | 南京理工大学 | 一种基于光学方法的骨质疏松症检测装置及方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5983486A (ja) * | 1982-11-04 | 1984-05-14 | Fuji Photo Film Co Ltd | 放射線画像のエネルギ−・サブトラクシヨン方法およびその方法に用いられる蓄積性螢光体シ−ト、蓄積性螢光体シ−ト積層体並びに蓄積性螢光体シ−トフイルタ積層体 |
| US5049746A (en) * | 1989-10-19 | 1991-09-17 | Fuji Photo Film Co., Ltd. | Method and apparatus for displaying energy subtraction images |
| US5122664A (en) * | 1990-04-27 | 1992-06-16 | Fuji Photo Film Co., Ltd. | Method and apparatus for quantitatively analyzing bone calcium |
-
1991
- 1991-03-28 JP JP3064289A patent/JP2641078B2/ja not_active Expired - Fee Related
-
1992
- 1992-03-27 US US07/858,634 patent/US5187731A/en not_active Expired - Lifetime
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014054301A (ja) * | 2012-09-11 | 2014-03-27 | Hitachi Aloka Medical Ltd | 医療用x線測定装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| US5187731A (en) | 1993-02-16 |
| JP2641078B2 (ja) | 1997-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2641078B2 (ja) | 骨塩定量分析方法 | |
| US6195413B1 (en) | Method and arrangement for detecting X-rays | |
| US5917877A (en) | Plain x-ray bone densitometry apparatus and method | |
| US5122664A (en) | Method and apparatus for quantitatively analyzing bone calcium | |
| Stewart et al. | Single‐exposure dual‐energy computed radiography | |
| Huddleston | Quantitative methods in bone densitometry | |
| JPS6244224A (ja) | 画像処理方法および装置 | |
| JPH0924039A (ja) | 骨塩定量分析方法および装置 | |
| JP2981701B2 (ja) | 放射線画像形成方法およびそれを用いた骨塩定量分析方法 | |
| JPH08266529A (ja) | 骨塩定量分析方法 | |
| JPH0599829A (ja) | 物質の定量測定法 | |
| JP2627098B2 (ja) | 骨塩定量分析方法および装置 | |
| JP2932019B2 (ja) | 骨塩定量分析方法 | |
| Attaelmanan et al. | Assessments of the physical performance of 2 generations of 2 direct digital intraoral sensors | |
| JPH06133223A (ja) | 骨部画像解析方法 | |
| JP2641077B2 (ja) | 骨塩定量分析方法 | |
| JP2694583B2 (ja) | 骨塩定量分析方法 | |
| CN117438050A (zh) | 放射线图像处理装置、方法及计算机可读取记录介质 | |
| JP2694578B2 (ja) | 骨塩定量分析方法 | |
| JP2694579B2 (ja) | 骨塩定量分析方法 | |
| RU2400141C2 (ru) | Способ определения минеральной плотности костной ткани | |
| KR101239133B1 (ko) | 의료용 방사선 촬영을 위한 디지털 팬텀, 디지털 팬텀을 이용한 의료 영상 처리 방법 및 시스템 | |
| JPH08274994A (ja) | 画像信号補正方法 | |
| JPH08266528A (ja) | 骨塩定量分析方法 | |
| JP2004223158A (ja) | X線撮像方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19970318 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |