JPH07289542A - 骨塩量測定方法および装置 - Google Patents
骨塩量測定方法および装置Info
- Publication number
- JPH07289542A JPH07289542A JP6082040A JP8204094A JPH07289542A JP H07289542 A JPH07289542 A JP H07289542A JP 6082040 A JP6082040 A JP 6082040A JP 8204094 A JP8204094 A JP 8204094A JP H07289542 A JPH07289542 A JP H07289542A
- Authority
- JP
- Japan
- Prior art keywords
- bone
- profile
- data
- projection
- value
- 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
- 238000000034 method Methods 0.000 title claims abstract description 71
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 88
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 42
- 239000011707 mineral Substances 0.000 claims description 42
- 230000005855 radiation Effects 0.000 claims description 38
- 230000005540 biological transmission Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 9
- 230000001678 irradiating effect Effects 0.000 claims description 8
- 238000010521 absorption reaction Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 abstract description 18
- 238000012937 correction Methods 0.000 abstract description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 26
- 229910052782 aluminium Inorganic materials 0.000 description 26
- 238000010586 diagram Methods 0.000 description 16
- 210000000236 metacarpal bone Anatomy 0.000 description 14
- 238000004458 analytical method Methods 0.000 description 11
- 230000037182 bone density Effects 0.000 description 11
- 238000002835 absorbance Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 238000010191 image analysis Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 230000000994 depressogenic effect Effects 0.000 description 2
- 210000002411 hand bone Anatomy 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 206010068975 Bone atrophy Diseases 0.000 description 1
- 206010061728 Bone lesion Diseases 0.000 description 1
- 208000001132 Osteoporosis Diseases 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000008468 bone growth Effects 0.000 description 1
- 230000004097 bone metabolism Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000002082 fibula Anatomy 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 210000002758 humerus Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 208000005368 osteomalacia Diseases 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 210000002320 radius Anatomy 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000002303 tibia Anatomy 0.000 description 1
- 210000000623 ulna Anatomy 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
同等の高速計測を行うことのできる骨塩量測定方法を提
供することを目的とする。 【構成】 第1の工程(170)で、投影プロファイル
を微分して、最大微分値と左極大値の区間および右極大
値と最小微分値の区間を算出する。次に、第2の工程
(180)で、各区間内の直線上に少なくとも3点の位
置をそれぞれ選択し、これらの各点の位置に対応した投
影プロファイルのデータ量をそれぞれ算出する。さら
に、第3の工程(190)で、算出された各データ量に
それぞれ最小二乗法を用いて投影プロファイルを適正化
して、左右のエッジの補正を行う。そして、適正化され
たプロファイルから被検骨の骨塩量を算出する。
Description
できる骨塩量測定方法および骨塩量測定装置に関し、特
に投影プロファイルの左右のエッジの補正方法に特徴の
ある測定方法および装置に関する。
は骨粗鬆症、骨軟化症等の骨病変の種類の判定やその症
状の進行度、治療時の効果の確認等の種々の骨塩量計測
を行う場合がある。
やDIP法などが知られている。MD法は、人間の両手
の中央にアルミスケール(スロープまたは階段状のアル
ミニウム板)を置いてX線撮影を行い、この撮影済みX
線フィルムから計測専用装置を用いて骨塩量を解析する
方法である。解析は、第2中手骨の中央をマイクロデン
シトメータで走査し、得られた濃度をアルミスケール濃
度に換算して、骨塩量に関する種々の指標を得る。一般
にそれらの指標の内、MCI(骨皮質幅の骨幅に対する
割合)とΣGs/D(骨幅あたりの骨陰影度)が最もよ
く用いられる。MD法の詳細は、“X線像による骨萎縮
度判定の試み、井上他、骨代謝第13巻、1980−
3.”等を参照されたい。
人間の両手を撮影したX線フィルムを、DIPシステム
の透過照明装置の上に置き、第2中手骨の画像を測定用
ビデオカメラを通じて高分解能画像処理装置に送り、M
D法と同等の骨量に関する種々の指標を自動解析によっ
て得る方法である。DIP法はMD法を改良した方法で
あり、MD法に比べて測定精度は落ちるものの、測定時
間はMD法の1/10以下に短縮される。DIP法の詳
細は、“新しい骨量測定法−DIP法とは?、林、Bo
ne&Joint Vol.3 No.5.”等を参照
されたい。
解析手法を具体的に説明すると次のようになる。まず第
2中手骨の中央に対してマイクロデンシトメータを走査
させながら、光を照射して得られる透過光の強度を測定
し、その走査された部位に対応した透過光の強度または
吸光度の線図を所定のチャート紙上に記入する。さらに
第2中手骨と共にX線撮影されたアルミスケールに対し
てマイクロデンシトメータを走査させ、得られた透過光
の強度または吸光度の線図についてもチャート紙に記入
する。このようにして得られたチャート紙上における第
2中手骨に関する吸光度とアルミスケールに関する吸光
度の各々の線図を、デジタイザを用いてコンピュータに
入力し、各点での第2中手骨の吸光度をアルミスケール
濃度に換算する。このようにして得られた図を用いて、
対象部位での骨形態を表す種々の指標がコンピュータ内
で計算され、計算結果が出力される。以上のようにMD
法による骨塩量の解析は、チャート紙上に第2中手骨に
関する吸光度などの線図を記入する必要があり非常に手
間がかかり問題であった。
は、チャート紙に線図を記入する代わりに測定用ビデオ
カメラで第2中手骨の画像を撮影して高分解能画像処理
装置に送っている。このため、非常に高速に解析を行う
ことができる。ところが、DIP法は精度面で以下の問
題があった。つまり、DIP法による解析は、測定用ビ
デオカメラを用いているため、X線フィルムエッジが線
源の大きさによりわずかにボケたり、撮像光学系による
ボケの影響で撮影が不鮮明になる恐れがあった。又、X
線フィルムのSN比が悪い場合もあり、データにノイズ
成分によるゆらぎがあったりした。このため、解析で得
られる指標の一つであるMCI(Metacarpal Index)の
値を予め小さな値に計算されるように補正して対処する
必要があった。しかしながら、この補正によってもMD
法なみに精度が上げることはできなかった。
法の問題を解決して、MD法と同等の精度でDIP法と
同等の高速計測を行うことのできる骨塩量測定方法およ
び骨塩量測定装置を提供することを目的とする。
に、本発明の骨塩量測定方法は、被検骨に放射線を照射
することによって得られる透過放射線像を撮影し、撮影
された透過放射線像の長手方向のほぼ中心を通る軸に垂
直な直線上で放射線吸収特性に基づく投影プロファイル
を測定することにより、この投影プロファイルから被検
骨の骨塩量を算出する骨塩量測定方法において、(a)
投影プロファイルを微分して、最大微分値と左極大値の
区間および右極大値と最小微分値の区間を算出する第1
の工程と、(b)両区間内の投影プロファイル上に少な
くとも3点ずつのデータを選択する第2の工程と、
(c)選択された少なくとも3点のデータから最小二乗
法を用いて曲線パターンを両区間で求め、これらの曲線
パターンで投影プロファイルのエッジ領域を適正化し、
適正化された適正プロファイルから被検骨の骨塩量を算
出する第3の工程とを備える。
被検骨に放射線を照射することにより得られる透過放射
線像を撮影する撮影手段と、(b)撮影手段によって撮
影された前記被検骨の透過放射線像を入力し、この透過
放射線像の長手方向のほぼ中心を通る軸に垂直な直線上
で放射線吸収特性に基づく投影プロファイルを測定する
測定手段と、(c)測定手段によって測定された投影プ
ロファイルを微分して、最大微分値と左極大値の区間お
よび右極大値と最小微分値の区間を算出する第1の演算
手段と、(d)第1の演算手段で算出された両区間内の
投影プロファイル上に少なくとも3点ずつのデータを選
択する選択手段と、(e)選択手段で選択された少なく
とも3点のデータから最小二乗法を用いて曲線パターン
を両区間で求め、これらの曲線パターンで投影プロファ
イルのエッジ領域を適正化し、適正化された適正プロフ
ァイルから被検骨の骨塩量を算出する第2の演算手段と
を備える。
過放射線像を撮影し、この透過放射線像の長手方向のほ
ぼ中心を通る軸に垂直な直線上で放射線吸収特性に基づ
く投影プロファイルを測定すると、被検骨は中央部と周
辺部で放射線吸収特性が異なるため、測定された投影プ
ロファイルは中央が陥没した山形の形状となる。この投
影プロファイルは、被検骨に放射線を照射するための放
射線源の大きさによるボケや、被検骨の透過放射線像を
撮影するための撮影撮像光学系によるボケのため、特に
左右のエッジ部分に誤差が生じ易い。
の補正を行っている。まず、第1の工程で、投影プロフ
ァイルを微分して、最大微分値と左極大値の区間および
右極大値と最小微分値の区間を算出している。これらの
区間が高精度の情報を含んでおり誤差が少ないからであ
る。次に、第2の工程で、両区間内の投影プロファイル
上に少なくとも3点ずつのデータを選択する。最低3点
の位置データがあれば、投影プロファイルの適正化が行
えるからである。さらに、第3の工程で、選択された少
なくとも3点のデータから最小二乗法を用いて曲線パタ
ーンを両区間で求め、これらの曲線パターンで投影プロ
ファイルのエッジ領域を適正化して、左右のエッジの補
正を行う。そして、適正化された適正プロファイルから
被検骨の骨塩量を算出する。
の骨塩量測定方法を適用した装置である。つまり、被検
骨の透過放射線像が撮影手段で撮影され、撮影された透
過放射線像の長手方向のほぼ中心を通る軸に垂直な直線
上で放射線吸収特性に基づく投影プロファイルが測定手
段で測定される。投影プロファイルは第1の演算手段に
与えられ、最大微分値と左極大値の区間および右極大値
と最小微分値の区間が算出される。そして、算出された
区間データは選択手段に与えられ、両区間内の投影プロ
ファイル上に少なくとも3点ずつのデータが選択され
る。さらに、選択された各データは第2の演算手段に与
えられ、最小二乗法によって投影プロファイルのエッジ
領域が適正化される。このように適正化された適正プロ
ファイルから被検骨の骨塩量が算出される。
参照して説明する。図1は、本実施例に係る骨塩量測定
装置の外観を示す斜視図である。同図に示すように本実
施例の骨塩量測定装置は、画像計測部10と画像解析部
20とから構成されている。画像計測部10は、X線フ
ィルム11を撮影する高解像度のビデオカメラ12と、
ビデオカメラ12を制御するカメラコントロールユニッ
ト13と、X線フィルム11を上部に載置して裏面から
照射する高周波点灯蛍光灯14と、X線フィルム11を
高周波点灯蛍光灯14の上部に載置させるフィルムロー
ディングステージ15とを備えている。ビデオカメラ1
2、カメラコントロールユニット13および高周波点灯
蛍光灯14は暗箱16の中に配設され、フィルムローデ
ィングステージ15は暗箱16の外から高周波点灯蛍光
灯14上部に導入できるよう配設されている。また、画
像解析部20は、ビデオカメラ12の出力画像を解析す
るコンピュータ21と、解析に必要なデータを入力する
キーボード22、マウス23と、解析結果を表示するモ
ニタ24と、解析結果を出力するプリンタ25とを備え
ている。
を示すブロック図である。暗箱15には、前述した機器
(ビデオカメラ12〜高周波点灯蛍光灯14)以外に、
X線フィルム11の周辺部を覆うマスク17と主電源回
路18が配設されており、主電源回路18はカメラコン
トロールユニット13および高周波点灯蛍光灯14に電
力を供給している。主電源回路18は、さらにコンピュ
ータ21およびモニタ24の各電源部211 ,241 に
も電力を供給している。コンピュータ21には、ビデオ
カメラ12の出力画像を取り込む画像取込装置21
2 と、取り込んだ画像を記憶するメインメモリ21
3 と、解析プログラム等が格納されたハードディスクユ
ニット214 とが内蔵されている。さらに画像解析結果
等を記録するフロッピーディスクユニット215 と、モ
ニタ24に表示する画像データを一時的に記憶するビデ
オRAM216 と、画像解析結果等をプリンタ25から
出力させるためのプリンタI/F(インタフェース)2
17 と、各機器(ハードディスクユニット214 〜プリ
ンタI/F217 )を制御するマイクロプロセッサ21
8 とが内蔵されている。
測定方法について、図3〜図6を参照して説明する。ま
ず、X線フィルム11の上にアルミスケール30と両手
掌を載せ、50KV、100mA、平均0.05秒でX
線撮影を行う(図3(a))。X線フィルム11を自動
現像装置で現像すると、中央にアルミスケール30のX
線像、左右に手骨のX線像が得られる。X線像が焼き付
けられたX線フィルム11はフィルムローディングステ
ージ15に載置され、暗箱16内の高周波点灯蛍光灯1
4上面に導入される。高周波点灯蛍光灯14の照射によ
ってX線フィルム11上に手骨及びアルミスケール30
のX線像が浮かび上がり、これらのX線像は高周波点灯
蛍光灯14上部のビデオカメラ12で撮影される(図3
(b)、ステップ100:図6)。撮影されたX線像の
画像データは画像取込装置212によって取り込まれ、
メインメモリ213 に記憶される。
下、ハードディスクユニット214に格納された解析プ
ログラムが起動され、メインメモリ213 に記憶され
た画像データの解析が行われる。まず、画像データ中の
アルミスケール画像の輝度が読み取られる(ステップ1
10:図6)。アルミスケール30には、10:1のス
ロープが付いており(空間的に10mmの距離が1mm
Alの厚さが増すことに相当する)、この幾何学的な条
件によってどの輝度値がアルミ厚に換算すると何mmの
厚さに相当するかが計算され、計算結果より変換テーブ
ルが作成される。そして、この変換テーブルに従ってメ
インメモリ213 上の画像データのすべての輝度データ
がアルミ厚に換算した輝度データに変換される(図3
(c)、ステップ120:図6)。
はビデオRAM216 上に移され、モニタ24に表示さ
れる。ここで操作者がマウス23を操作して、モニタ2
4に表示された画像データ中の所望の測定部位(通常は
右手第2中手骨)の上下の端を指定すると、指定した2
点の中間に特定の幅を持った(通常は全長の10%の幅
とする)測定ウインドウがモニタ24に表示される(図
4(d)、ステップ130:図6)。
リ213 上の画像データに対して、骨長手方向と直交す
る投影プロファイルが複数本読み出され、これらの投影
プロファイルがメインメモリ213 の画像データと異な
る領域に記憶される(図4(e))。読み出される投影
プロファイルの本数は通常30〜40本程度である。投
影プロファイルをグラフにすると、中央が陥没した山形
の形状となる(図5(f)、ステップ140:図6)。
ファイルが一本ずつ読み出され、それぞれの投影プロフ
ァイルの空間的な左右のエッジ点と輝度方向におけるバ
ックグラウンドのオフセット分が検出される。この検出
値よりΣGs(図中斜線部)が求められ、ΣGsをエッ
ジ間距離Dで除算することにより、ΣGs/Dが算出さ
れる(図5(g))。ΣGs/Dは、骨塩定量に用いる
骨密度パラメータの一つである。さらに、その他の骨密
度パラメータであるMCI,ΔGsmax ,ΔGsmin ,
D,dがそれぞれ算出される(ステップ150:図
6)。これらの骨密度パラメータは、図7に示すよう
に、投影プロファイルのグラフから求められる。これら
の骨密度パラメータの概要を次に示す。
計算され、さらに、測定精度を向上させるため、各々の
パラメータごとに平均化して、測定シートという電子フ
ァイルとしてハードディスクユニット214 もしくはフ
ロッピーディスクユニット215 に記録される。さら
に、必要に応じて算出された骨密度パラメータのデータ
をプリンタ25から出力させることもできる。ステップ
100〜ステップ150までの処理は、操作者がキーボ
ード22あるいはマウス23を用いて、処理終了の入力
を行うまで繰り返される(ステップ160:図6)。
程での、投影プロファイルの左右のエッジ点の検出法に
ある。つまり、本実施例の検出法を用いれば、従来のD
IP法で問題となっていたビデオカメラ12で撮影する
ことによる撮影画像のエッジのボケをほぼ完全に補正す
ることができるのである。このエッジ検出法の原理につ
いて、図8,9を参照して説明する。図8は第2中手骨
の構造モデルの透過放射線像をアルミ厚で換算する処理
を示す図である。また、図9はデータ処理区間の求め方
を示す図である。
の骨のモデル(図8(a))を用いる。この骨のモデル
に対して上部からX線を照射すると、左右対象な輝度分
布を持つ透過放射線像(図8(b))がX線フィルム1
1上に得られる。次に、X線フィルム11の輝度とアル
ミ厚との関係からなる変換テーブル(図8(c))に透
過放射線像を通すことにより、投影プロファイル(図8
(d))が得られる。このようにして得られた投影プロ
ファイルのエッジを求めればよいが、実際のエッジは、
X線源の大きさによるボケや撮像光学系によるボケのた
め、なまった状態になってしまう(図9(a))さらに
充分なSNが得られない場合が多く、ショットノイズが
重畳した状態になっている。本実施例の方式は、投影プ
ロファイルの微分最大値と左極大値の区間(図9(b)
(c)のaの区間)と、右極大値と微分最小値の区間
(図9(b)(c)のbの区間)とが、骨測定に必要な
大局的かつ高精度な情報を含んでいることに着目し、こ
れらの各区間の少なくとも3点のデータを用いてエッジ
を算出している。
10のフローチャートを用いて説明する。まず、図6の
ステップ140で作成した投影プロファイルを微分し
て、微分最大値と左極大値の区間、および右極大値と微
分最小値の区間を算出する(ステップ170)。これら
の区間が高精度の情報を含んでおり誤差が少ないからで
ある。次に、それぞれの区間から投影プロファイル上に
少なくとも3点ずつのデータを選択する(ステップ18
0)。最低3点の位置データがあれば、投影プロファイ
ルの適正化が行えるからである。さらに、選択された少
なくとも3点のデータから最小二乗法を用いて曲線パタ
ーンを両区間で求め、これらの曲線パターンで投影プロ
ファイルのエッジ領域を適正化して、左右のエッジの補
正を行う。そして、適正化された適正プロファイルから
骨密度パラメータを算出する(ステップ190)。この
ように適正プロファイルから骨密度パラメータを算出し
ているので、エッジのボケやノイズによるゆらぎに対す
る補正がほぼ完全に行え、極めて正確な画像解析を行う
ことができる。
ルとして、楕円形状のモデルを用いた例を図11に示
す。楕円形状の構造モデルを図11(a)に、アルミ厚
換算した投影プロファイルを図11(b)にそれぞれ示
す。ここで、楕円形状のモデルを用いたのは、縦軸がア
ルミ厚換算した時点で、空間座標系である横軸とアルミ
換算厚み輝度の縦軸ではスケールが異なるためである。
´の区間だけである。式を、 y=b/a(a2 −x2 )1/2 …… の様に変形し、このy座標よりアルミの厚みをlとした
場合、yの2倍がlの大きさとなる。そこで式のyに
2lを代入すると、 l=2b/a(a2 −x2 )1/2 …… となる。さらに式を変形すれば、 l2 =4b2 /a2 (a2 −x2 ) …… となる。この式のx2 をpに、l2 をqにそれぞれ置
き換えると、 q=−Ap+B …… が得られる。この式にa´,b´の区間データを用い
てフィッティングすれば理想モデルの式が求まるが、実
際にはエッジがなまるため、この区間を用いることがで
きない。従って、先に述べた図8(b)のaとbの区間
で最小2乗法を用いて式を求め、式から式を求め
ればよい。
から直線式あるいは曲線式を近似する最も一般的な手法
である。具体的には、複数点(xi ,yi )から一つの
直線 y=ax+b までのy軸に平行な直線に沿っての距離di の平方の和
直線式を得るものである。曲線式についても同様であ
る。
肉組織のアルミ換算の厚みがバックグラウンド(背景
部)として重畳されている。したがって、左右のエッジ
を算出して骨密度パラメータを求めるためには、投影プ
ロファイルから背景部を取り除かなければならない。こ
の処理の概要を図12のフローチャート、及び図13,
14のグラフに示す。
称に傾いている場合が多い。したがって、微分最大値の
左側及び微分最小値の右側のそれぞれ一定距離以上離れ
た部分の区間で、左右別々に背景部の境界を成す直線を
最小2乗法により求め、これらを直線L1,L2とする
(ステップ200:図12)。さらに、これらの直線と
微分最大値及び微分最小値の垂線の交点の高さをL-OFFS
ET,R-OFFSETとし、微分最大値でのアルミ換算厚みl1
からL-OFFSETを差し引いた値を式のlに代入して最小
2乗法にて式の係数(a,b)を求めることにより、微
分最大値を通る曲線L3が得られる。同様に、微分最小
値でのアルミ換算厚みl2 からR-OFFSETを差し引いた値
を式のlに代入して最小2乗法にて式の係数(a,
b)を求めることにより、微分最小値を通る曲線L4が
得られる(ステップ210:図12)。
曲線L3の交点P1と、直線L2と曲線L4の交点P2
とを求め、P1とP2とを結んで直線L5とし、P1と
P2との間隔をDとする(ステップ220:図12)。
曲線L3,L4及び直線L5と投影プロファイルとで囲
まれた図中の斜線部分をアルミ換算厚の総和ΣGsとす
る。ここで本方式では、微分最大値と微分最小値の区間
では投影プロファイルを積分して区間内の面積を算出
し、P1 と微分最大値の区間では曲線L3と直線L5で
囲まれた面積を、微分最小値とP2 の区間では曲線L4
と直線L5で囲まれた面積をそれぞれ算出し、このよう
に算出した値を加算することにより、総和ΣGsを求め
ることができる。そして、総和ΣGsをエッジ間距離D
で除算することにより、ΣGs/Dが算出される。さら
に、その他の骨密度パラメータであるMCI,ΔGs
max ,ΔGsmin ,D,dがそれぞれ算出される(ステ
ップ230:図12)。
となく、種々の変形が可能である。例えば、X線フィル
ム11を用いる代わりに、被検骨にX線やγ線の放射線
を照射して得られる透過放射線量自体を二次電子増倍管
等の検知手段によって検知された画像データをコンピュ
ータ20で入力してもよい。また、X線撮影装置(図示
せず)にイメージセンサを設け、そのセンサからの画像
信号をA/D変換して得られる画像データをコンピュー
タ20で入力してもよい。
手骨が用いられるが、上腕骨、撓骨、尺骨、大腿骨、脛
骨、腓骨等の長骨であってもよい。
ば、放射線画像のようなノイズを含んだ画像に対しても
測定された投影プロファイルの内、比較的誤差の少ない
最大微分値と左極大値の区間および右極大値と最小微分
値の区間のデータを用いて、投影プロファイルを適正化
し、適正化されたプロファイルを用いて被検骨の骨塩量
を算出している。
の放射線源の大きさによるボケや、被検骨の透過放射線
像を撮影するための撮影撮像光学系によるボケのために
投影プロファイルのエッジ部分に誤差を含んでいる場合
でも、被検骨の骨塩量のパラメータであるΣGs/Dや
MCIをMD法と同等の精度で算出することができる。
放射線像に基づいて測定しており、所定のチャート紙上
に吸光度の線図を記入する必要がないので、DIP法と
同等の高速計測を行うことができる。
視図である。
ク図である。
る。
る。
る。
トである。
ミ厚で換算する処理を示す図である。
ーチャートである。
ミ厚で換算する処理を示す図である。
密度パラメータを求める処理を示すフローチャートであ
る。
を示す図である。
る。
カメラ、13…カメラコントロールユニット、14…高
周波点灯蛍光灯、15…フィルムローディングステー
ジ、16…暗箱、17…マスク、18…主電源回路、2
0…画像解析部、21…コンピュータ、211 ,241
…電源部、212 …画像取込装置、213…メインメモ
リ、214 …ハードディスクユニット、215 …フロッ
ピーディスクユニット、216 …ビデオRAM、217
…プリンタI/F、218 …マイクロプロセッサ、22
…キーボード、23…マウス、24…モニタ、25…プ
リンタ、30…アルミスケール。
Claims (4)
- 【請求項1】 被検骨に放射線を照射することによって
得られる透過放射線像を撮影し、撮影された透過放射線
像の長手方向のほぼ中心を通る軸に垂直な直線上で放射
線吸収特性に基づく投影プロファイルを測定することに
より、この投影プロファイルから前記被検骨の骨塩量を
算出する骨塩量測定方法において、 前記投影プロファイルを微分して、最大微分値と左極大
値の区間および右極大値と最小微分値の区間を算出する
第1の工程と、 前記両区間内の前記投影プロファイル上に少なくとも3
点ずつのデータを選択する第2の工程と、 選択された少なくとも3点のデータから最小二乗法を用
いて曲線パターンを前記両区間で求め、これらの曲線パ
ターンで前記投影プロファイルのエッジ領域を適正化
し、適正化された適正プロファイルから前記被検骨の骨
塩量を算出する第3の工程とを備えることを特徴とする
骨塩量測定方法。 - 【請求項2】 前記第3の工程では、前記最大微分値と
前記最小微分値の区間の左右の外側の前記投影プロファ
イル上にそれぞれ複数点のデータを選択し、これらのデ
ータから最小二乗法を用いて左右の背景の直線パターン
をそれぞれ求め、これらの直線パターンと前記適正プロ
ファイルとの2つの交点を左右のエッジとして、これら
のエッジ間の適正プロファイルから前記被検骨の骨塩量
を算出していることを特徴とする請求項1記載の骨塩量
測定方法。 - 【請求項3】 被検骨に放射線を照射することにより得
られる透過放射線像を撮影する撮影手段と、 前記撮影手段によって撮影された前記被検骨の透過放射
線像を入力し、この透過放射線像の長手方向のほぼ中心
を通る軸に垂直な直線上で放射線吸収特性に基づく投影
プロファイルを測定する測定手段と、 前記測定手段によって測定された前記投影プロファイル
を微分して、最大微分値と左極大値の区間および右極大
値と最小微分値の区間を算出する第1の演算手段と、 前記第1の演算手段で算出された前記両区間内の前記投
影プロファイル上に少なくとも3点ずつのデータを選択
する選択手段と、 前記選択手段で選択された少なくとも3点のデータから
最小二乗法を用いて曲線パターンを前記両区間で求め、
これらの曲線パターンで前記投影プロファイルのエッジ
領域を適正化し、適正化された適正プロファイルから前
記被検骨の骨塩量を算出する第2の演算手段とを備える
ことを特徴とする骨塩量測定装置。 - 【請求項4】 前記第2の演算手段では、前記最大微分
値と前記最小微分値の区間の左右の外側の前記投影プロ
ファイル上にそれぞれ複数点のデータを選択し、これら
のデータから最小二乗法を用いて左右の背景の直線パタ
ーンをそれぞれ求め、これらの直線パターンと前記適正
プロファイルとの2つの交点を左右のエッジとして、こ
れらのエッジ間の適正プロファイルから前記被検骨の骨
塩量を算出していることを特徴とする請求項3記載の骨
塩量測定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP6082040A JP2793502B2 (ja) | 1994-04-20 | 1994-04-20 | 骨塩量測定方法および装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP6082040A JP2793502B2 (ja) | 1994-04-20 | 1994-04-20 | 骨塩量測定方法および装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH07289542A true JPH07289542A (ja) | 1995-11-07 |
JP2793502B2 JP2793502B2 (ja) | 1998-09-03 |
Family
ID=13763413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP6082040A Expired - Lifetime JP2793502B2 (ja) | 1994-04-20 | 1994-04-20 | 骨塩量測定方法および装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2793502B2 (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003515377A (ja) * | 1999-12-01 | 2003-05-07 | ロナルド・イー・マッシー | 歯および整形のデンシトメトリーモデリングのシステムおよび方法 |
JP2009106425A (ja) * | 2007-10-29 | 2009-05-21 | Aloka Co Ltd | 前腕用x線測定装置 |
JP2012143387A (ja) * | 2011-01-12 | 2012-08-02 | Hiroshima Univ | 骨粗鬆症診断支援装置及び骨粗鬆症診断支援プログラム |
JP2013504341A (ja) * | 2009-09-11 | 2013-02-07 | ストラックスコープ ピーティワイ リミテッド | 解析方法及び解析システム |
CN103239248A (zh) * | 2012-02-09 | 2013-08-14 | 富士胶片株式会社 | 分析对象部位的确定方法、图像处理装置及记录介质 |
JP2014054476A (ja) * | 2012-09-14 | 2014-03-27 | Fujifilm Corp | 骨塩量計測装置および方法 |
JP2014054478A (ja) * | 2012-09-14 | 2014-03-27 | Fujifilm Corp | 骨塩量計測装置および方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9427202B2 (en) | 2014-07-03 | 2016-08-30 | Fujifilm Corporation | Arithmetic processor and bone density measuring device |
-
1994
- 1994-04-20 JP JP6082040A patent/JP2793502B2/ja not_active Expired - Lifetime
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003515377A (ja) * | 1999-12-01 | 2003-05-07 | ロナルド・イー・マッシー | 歯および整形のデンシトメトリーモデリングのシステムおよび方法 |
JP2009106425A (ja) * | 2007-10-29 | 2009-05-21 | Aloka Co Ltd | 前腕用x線測定装置 |
JP2016104159A (ja) * | 2009-09-11 | 2016-06-09 | ストラックスコープ ピーティワイ リミテッドStraxcorp Pty Ltd | 解析方法及び解析システム |
JP2013504341A (ja) * | 2009-09-11 | 2013-02-07 | ストラックスコープ ピーティワイ リミテッド | 解析方法及び解析システム |
US9064320B2 (en) | 2009-09-11 | 2015-06-23 | Straxcorp Pty Ltd | Method and system for image analysis of selected tissue structures |
US9378566B2 (en) | 2009-09-11 | 2016-06-28 | Straxcorp Pty Ltd | Method and system for image analysis |
US9566038B2 (en) | 2009-09-11 | 2017-02-14 | Straxcorp Pty Ltd | Method and system for image analysis |
JP2018008094A (ja) * | 2009-09-11 | 2018-01-18 | ストラックスコープ ピーティワイ リミテッドStraxcorp Pty Ltd | 解析方法及び解析システム |
US9877691B2 (en) | 2009-09-11 | 2018-01-30 | Straxcorp Pty Ltd | Method and system for image analysis |
US10182781B2 (en) | 2009-09-11 | 2019-01-22 | Straxcorp Pty Ltd | Method and system for image analysis |
JP2012143387A (ja) * | 2011-01-12 | 2012-08-02 | Hiroshima Univ | 骨粗鬆症診断支援装置及び骨粗鬆症診断支援プログラム |
CN103239248A (zh) * | 2012-02-09 | 2013-08-14 | 富士胶片株式会社 | 分析对象部位的确定方法、图像处理装置及记录介质 |
JP2013162822A (ja) * | 2012-02-09 | 2013-08-22 | Fujifilm Corp | 骨塩定量分析における分析対象部位の特定方法および、この方法を実施するための画像処理装置並びに記録媒体 |
JP2014054476A (ja) * | 2012-09-14 | 2014-03-27 | Fujifilm Corp | 骨塩量計測装置および方法 |
JP2014054478A (ja) * | 2012-09-14 | 2014-03-27 | Fujifilm Corp | 骨塩量計測装置および方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2793502B2 (ja) | 1998-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4653737B2 (ja) | 3次元再構成のための放射線画像法並びにその方法を実施するコンピュータプログラムおよび装置 | |
US7069066B2 (en) | Bone densitometer providing improved longitudinal studies | |
US6038281A (en) | Bone densitometer with improved point characterization | |
Sharp et al. | Computer‐based methods for measuring joint space and estimating erosion volume in the finger and wrist joints of patients with rheumatoid arthritis | |
EP0648467B1 (en) | Osteometry and osteometric apparatus | |
JPH0331061B2 (ja) | ||
JP2002263097A (ja) | X線ct装置 | |
WO1990009761A1 (en) | Method of measuring bone shape, apparatus therefor and system for evaluating bone | |
JP2793502B2 (ja) | 骨塩量測定方法および装置 | |
JP2007151609A (ja) | 骨計測装置及びその方法 | |
JP2006334319A (ja) | X線ct装置とその前処理方法、及びデータ作成装置とその方法、並びに制御プログラム | |
US7046834B2 (en) | Method for measuring bone mineral density by using X-ray image | |
JPH10509074A (ja) | フィルムカセット付き骨濃度計 | |
JP2983421B2 (ja) | 骨計測方法及び装置 | |
KR0164600B1 (ko) | 뼈형태계측장치 및 뼈평가시스템 | |
JP2732979B2 (ja) | 骨計測装置 | |
JP2620396B2 (ja) | 骨計測方法及び装置 | |
JP2627098B2 (ja) | 骨塩定量分析方法および装置 | |
JP2812875B2 (ja) | 骨計測方法 | |
JP2771452B2 (ja) | 骨計測方法及びそのための装置 | |
JPH02220640A (ja) | 骨計測装置及び骨評価システム | |
JP2763449B2 (ja) | 骨計測方法及び装置 | |
JPH09262232A (ja) | 骨計測方法 | |
JP2812873B2 (ja) | 骨計測方法 | |
JPH0595940A (ja) | 骨計測装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313532 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090619 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100619 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110619 Year of fee payment: 13 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120619 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120619 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130619 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130619 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140619 Year of fee payment: 16 |
|
EXPY | Cancellation because of completion of term |