JP2004147863A - X-ray bone density measuring apparatus - Google Patents

X-ray bone density measuring apparatus Download PDF

Info

Publication number
JP2004147863A
JP2004147863A JP2002316490A JP2002316490A JP2004147863A JP 2004147863 A JP2004147863 A JP 2004147863A JP 2002316490 A JP2002316490 A JP 2002316490A JP 2002316490 A JP2002316490 A JP 2002316490A JP 2004147863 A JP2004147863 A JP 2004147863A
Authority
JP
Japan
Prior art keywords
bone density
body fat
fat percentage
distribution
bone
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
Application number
JP2002316490A
Other languages
Japanese (ja)
Other versions
JP3974842B2 (en
Inventor
Hiroyuki Toho
弘之 東方
Kazumasa Ise
一昌 伊勢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Aloka System Engineering Co Ltd
Original Assignee
Aloka Co Ltd
Aloka System Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Aloka Co Ltd, Aloka System Engineering Co Ltd filed Critical Aloka Co Ltd
Priority to JP2002316490A priority Critical patent/JP3974842B2/en
Publication of JP2004147863A publication Critical patent/JP2004147863A/en
Application granted granted Critical
Publication of JP3974842B2 publication Critical patent/JP3974842B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To obtain an image which represents the relationship between a bone density and a body fat percentage in an easily understandable way by using an X-ray bone density measuring apparatus. <P>SOLUTION: A patient is scanned with high and low energy X-rays to obtain the attenuation of the X-rays at each pixel for each energy (S10 and S12). From the ratio of attenuations of each energy (attenuation ratio) (S14), whether each pixel represents a bone part or a soft tissue (S16) is discriminated. About a bone part pixel, the bone density is calculated by a DXA method (S18). About a soft tissue pixel, the body fat percentage is calculated by using a function for converting the attenuation ratio to the body fat percentage (S26). A bone part image expressing the bone density of each pixel of the bone part in a first color (S20) and a soft tissue image expressing each pixel of the soft tissue with a second color (S26) are combined according to the classification obtained at step S16. Thus, a synthetic image which displays bone density distribution and body fat ratio distribution as a list is prepared (S28). <P>COPYRIGHT: (C)2004,JPO

Description

【0001】
【発明の属する技術分野】
本発明は、X線を用いて被検体の骨密度分布を測定する装置に関する。
【0002】
【従来の技術】
X線を被検体に透過させて骨密度を演算するX線骨密度測定装置は、従来より骨の疾病診断やその予防のために多くの医療機関で利用されている。周知のように骨はX線を減弱するカルシウムなどのミネラル(骨塩)で構成されており、X線透過の際の減衰量(減弱率)を求めることによって、単位面積当たりの骨密度が演算される。骨は筋肉や脂肪などの軟部組織に覆われており、その部分でもX線の減弱が生じる。そのような軟部組織による影響を排除するために、従来の骨密度測定装置では、2種類のエネルギーのX線が利用されている。これは各エネルギーで演算された減衰量に基づき連立方程式を解くことによって骨部分及び軟部組織分それぞれの減衰量を特定するものである。この手法は二重X線吸収測定法(DXA法)と呼ばれる。このDXA法による骨密度計算の原理を図5を用いて説明する。図5には、(b)に示すように骨部12の周囲を軟部組織14が覆っている被検体10に対し、高低2種類のエネルギーのX線を透過させた時の減衰量の分布を(a)に示す。図から分かるように、骨部も軟部組織も低エネルギーX線の方が高エネルギーX線よりも減衰量が高いが、両エネルギー間での減衰量の比は骨部と軟部組織とで明らかに異なる。そこで、DXA法では、軟部組織についての高エネルギーX線の減衰量に対する低エネルギーX線の減衰量の比αを求め、低エネルギーX線の減衰量から、高エネルギーX線の減衰量にαを掛けたものを減算することで、骨部のみの減衰量を求め、これを骨密度(この場合は単位面積当たりの骨塩量)に換算している。特許文献1及び2には、このような手法を用いた骨密度測定装置の一例が示されている。
【0003】
この方式の骨密度測定装置では、例えばベッドに寝た被検体の上方から下方へとX線を照射して検出することで、被検体を2次元投影した場合の各点の骨密度を求めることができ、この2次元分布を表示することができる。
【0004】
また、従来のDXA法の骨密度測定装置には、被検体の外形形状等を把握しやすくするなどの目的で、骨部だけでなく軟部組織を表示できるものもある。すなわち、一般的な骨密度算出では、上述のごとく軟部組織の値が0になるような比αを用いているが、このαの値をそれより小さい値とすれば、軟部組織も正の値となるので、このような小さいαを用いて表示画像を形成すれば、軟部組織をも表した表示が得られる。
【0005】
【特許文献1】
特許第2735507号明細書
【特許文献2】
特許第3256667号明細書
【0006】
【発明が解決しようとする課題】
近年、ラットやマウス等の動物を用いた投薬効果確認等のために、X線骨密度測定装置を用いることも行われている。例えば、やせ薬による筋肉や脂肪の量の変化が骨に与える影響を調べるために、X線骨密度測定装置により骨密度分布を測定する例が知られている。この用途では、体脂肪率の分布と骨密度分布との関係が分かると非常に便利である。しかしながら、従来の骨密度測定装置にはそのような両分布を関連づけて表示する機能を持つものは存在しない。
【0007】
ここで、上記従来技術に示した軟部組織表示は、骨密度分布だけでなく、骨部の周囲の軟部組織もある程度表示できるが、この表示に示される軟部組織の値は、軟部組織によるX線減衰量を骨密度のスケールに換算したものであって、体脂肪率を表してはいない。例えば、体脂肪率は一般に体側部(脇腹など)の方が高い傾向があるが、体側部は体の厚みが小さいためX線の減弱が小さくなるので、上記の軟部組織表示では体側部の値が被検体外部と同等の小さい値となり、体脂肪率の状況を表さない。
【0008】
本発明はこのような問題に鑑みなされたものであり、骨密度と体脂肪率との関係を分かりやすく示すことができるX線骨密度測定装置を提供することを目的とする。
【0009】
【課題を解決するための手段】
上記目的を達成するため、本発明に係る装置は、被検体を透過したX線を検出し、この検出結果に基づきX線の減衰量の2次元分布を求める測定手段と、前記2次元分布から被検体の骨部と軟部組織とを弁別する弁別手段と、前記2次元分布から前記骨部の骨密度分布を求める骨密度算出手段と、前記2次元分布から前記軟部組織の体脂肪率分布を求める体脂肪率算出手段と、前記骨部の骨密度分布と前記軟部組織の体脂肪率分布とを一画像に合成して出力する合成手段とを含む。
【0010】
この構成では、減衰量の2次元分布から、骨部については骨密度分布を、軟部組織については体脂肪率分布をそれぞれ求め、両者を合成して1枚の画像を生成する。したがって、ユーザはその画像により骨密度分布と体脂肪率分布を一覧することができるので、骨密度と体脂肪率の関係の分析が容易になる。
【0011】
好適な態様では、前記測定手段は、異なる2つのエネルギーのX線についてそれぞれ前記減衰量の2次元分布を求め、前記骨密度算出手段は、前記各エネルギーのX線についての前記2次元分布から二重X線吸収測定法に従って骨密度分布を求め、前記体脂肪率分布算出手段は、各エネルギーのX線の減衰量同士の比と体脂肪率との対応関係の情報を有し、前記2次元分布から各点での各エネルギーのX線の減衰量同士の比を求め、その比に対応する体脂肪率を前記対応関係の情報から求めることで前記体脂肪率分布を求める。
【0012】
また別の好適な態様では、前記合成手段は、前記骨密度分布と前記体脂肪率分布とに、互いに異なる表示形態範囲を割り当て、骨密度と体脂肪率とをそれぞれ対応する表示形態範囲内の表示形態で表現する。この態様によれば、骨密度と体脂肪率とが互いに区別しやすい画像を生成できる。
【0013】
更に別の好適な態様では、前記体脂肪率分布算出手段は、前記2次元分布から各点での各エネルギーのX線の減衰量同士の比を求めるに当たり、前記異なる2つのエネルギーのX線のうち、エネルギーが低い方のX線に対する前記2次元分布の減衰量の値を補正する。この態様によれば、軟部組織の特性にあった適切な減衰量補正を施すことができる。
【0014】
【発明の実施の形態】
以下、本発明の実施の形態(以下実施形態という)について、図面に基づいて説明する。
【0015】
本実施形態のX線骨密度測定装置は、被検体の骨密度分布と体脂肪率分布を求め、これらを1枚の画像に合成して表示する機能を備える。この装置は、ハードウエア構成としては従来のDXA法による骨密度測定装置とほぼ同等の構成でよく、体脂肪率分布の計算や合成画像の生成は典型的にはソフトウエア的に実現される(もちろんこれら計算や画像生成の処理の一部又は全部をハードウエア回路化することも不可能ではない)。
【0016】
そこでまず、本発明に係る骨密度・体脂肪率の同時表示方式が適用されるX線骨密度測定装置のハードウエア構成について、図1を参照して説明する。
【0017】
この骨密度測定装置は、例えばラット等の被検体10の胴体の骨密度分布を測定する測定装置であり、固定配置されるベッド16上には被検体10が載置される。実際には、小動物をベッド16に固定するための機構が存在するが、図では繁雑さを避けるため図示を省略している。被検体10は生体であり、骨12は、筋肉や脂肪などからなる軟部組織14によって覆われている。本実施形態において、ベッド16の下側にはX線発生装置18が設けられ、ベッド16の上方にはX線検出装置22が配置されている。X線発生装置18は必要に応じて高エネルギーX線又は低エネルギーX線を発生する装置であり、X線検出装置は被検体10を透過したX線100を検出する装置である。これらのX線発生装置18及びX線検出装置22は走査装置24によって両者一体的に水平方向に駆動される。例えば、ジグザグ走査により被検体の胴体全域を走査し、透過X線を検出する。
【0018】
制御部26は、本装置の動作制御を行うものであり、また後述する各種演算を実行する。骨密度や体脂肪率の計算、合成表示画像の生成などの処理も、この制御部26の演算処理によって実現される。記憶部30には、制御部26の制御・演算処理のための各種のプログラムやデータが保存される。記憶部30は例えばEEPROMなどで構成される。表示器32には制御部26における演算結果などが表示される。
【0019】
図1に示したように、被検体10を測定する場合、まず、被検体10がベッド16上に載置され、その後、走査装置24によってX線発生装置18及びX線検出装置22を水平方向にスキャンさせながら、例えば数ミリ秒程度の所定周期で交互に高エネルギーX線及び低エネルギーX線の照射を行う。これをX線検出装置22で検出することにより、水平面内での各点ごとに、被検体10を透過した高エネルギーX線の計数率(Ih)と低エネルギーX線の計数率(Il)が各々求められる。各点での減衰量は、実測した計数率Ih,Ilで、それぞれ基準計数率Ioh,Iolを除したものの自然対数で表される。すなわち、高エネルギーX線の減衰量Rhはln(Ioh/Ih),低エネルギーX線の減衰量Rlはln(Iol/Il)で求められる。なお、基準計数率Ioh,Iolは、ベッド16上の被検体がない位置に、被検体実測時と同じ条件でX線発生装置18から高エネルギーX線及び低エネルギーX線をそれぞれ照射した時に、X線検出装置22の検出結果から得られる計数率を用いることができる。このようにして求めた高エネルギーX線及び低エネルギーX線の各々の減衰量が、この装置による測定の生データとなる。以下、この生データをもとに、各点の骨密度や体脂肪率の計算を行う。
【0020】
ここで、骨密度の計算は周知のDXA法のアルリゴズムに従って行えばよいので説明を省略し、図2を参照して体脂肪率の計算方式を説明する。
【0021】
体脂肪率の計算のためには、X線透過に関して軟部組織を模擬するファントムを用いてキャリブレーションを行い、減衰量比と体脂肪率との関係を示す減衰量比−体脂肪率関数を求める。このキャリブレーションには、2種類の異なる体脂肪率に対応するファントムを用いる。例えば、体脂肪率0%(すなわち筋肉100%)に対応するファントム(ファントムFat0と呼ぶ)と、体脂肪率100%(筋肉0%)に対応するファントム(ファントムFat100と呼ぶ)の2種類を用いる。そして、体脂肪率0%のファントムFat0をベッド16に載置して高エネルギーX線及び低エネルギーX線を照射し、それぞれの減衰量Rh,Rlを計算し、両者の比Rl/Rhを計算する。この比を減衰量比と呼ぶ。ファントムFat0について求めた減衰量比は、図2ではRR0と示している。図2は、縦軸に体脂肪率、横軸に減衰量比をとってそれら両者の関係を示したものであり、体脂肪率0%に対応する点はAで示されている。同様に、ファントムFat100を測定し、体脂肪率100%に対応する減衰量比RR100を求める。これに対応する点は、図2には点Bとして示している。そして、体脂肪率0%と100%の2点を除く範囲については、これら2点をパラメータとする関数により近似する。この関数のことを換算関数fと呼ぶ。図2の例は、それら2点の間を直線の換算関数fで近似した場合の例である。もちろん、実験等によりより適切な換算関数を見出して利用してもよい。いずれにしてもキャリブレーション処理は、利用する換算関数を特定するのに必要な数の、体脂肪率の異なるファントムを用いて行えばよい。
【0022】
このようにして求められた換算関数fは記憶部30に記憶され、利用される。体脂肪率の算出では、被検体100を実測して減衰量比RRxを求め、換算関数fにこのRRxを代入し、体脂肪率の値Fを求めればよい。なお、換算関数は、上述のように関数形を規定するパラメータの組として記憶してもよいが、各減衰量比に対応する体脂肪率の値を登録したテーブルなどの形で実装してももちろんよい。
【0023】
次に、図3を参照して、本発明に係る装置で行われる処理の手順を説明する。
【0024】
まず、上述のようにX線発生装置18及びX線検出装置22を走査して、被検体100の各点の高エネルギーX線及び低エネルギーX線の減衰量Rh,Rlを求め(S10)、この測定結果に基づきRh,Rlの2次元分布を表すマップをそれぞれ作成する(S12)。このマップは、画像の各画素ごとに減衰量(Rh又はRl)の値を持っている。これらマップの組が測定結果の生データとなる。以降は、これらマップの組をもとに、骨密度測定の処理と体脂肪率測定の処理を実行する。
【0025】
骨密度測定でも体脂肪率測定でも、減衰量の生データから、各画素ごとに減衰量比RR=Rl/Rhを計算する(S14及びS20)。ただし、減衰量比の計算は、単に生データにおける高低各エネルギーの減衰量同士の比をとるのではなく、骨密度測定の場合(S14)と体脂肪率測定(S20)の場合とで、それぞれ適切な補正演算を加える。この補正演算は、一つには体厚の違いによる測定結果(骨密度又は体脂肪率)のずれを補正する(体厚補正とも呼ばれる)ためのものである。ここで、体厚とはX線が被検体を透過する長さのことであり、理想的には体厚が異なっても、同じ骨密度(又は体脂肪率)のものを測定した場合は、同じ減衰量比RRが求められるべき(すなわち両者の関係は直線比例関係)であるが、現実にはそうはならず、両者の関係は非直線的になる。そこで、できるだけ理想に近くなるようにRh,Rlの生データに対して補正を加えた上で、両者の比をとるようにしている。ただし、DXA法の骨密度測定装置の場合、装置構成上被検体各部の体厚を測定することはできない(体厚測定のために専用の測定機構を設けることも考えられるが、コスト的に見合わない)ので、単純に生データの減衰量を、体厚補正された減衰量に換算する構成を採っている。すなわち、減衰量Rh及びRlは、体厚に関して比例関係ではないものの単調に変化することから、生の減衰量の値自身が体厚とほぼ1対1の対応関係にあるとの仮定できるために、このような換算方式が可能である。この換算のためには、同じ材質で厚みの異なるファントムを複数用意し、それらを高低各エネルギーのX線で測定し、その時の減衰量比Rl/Rhが体厚によらずできるだけ均一となるように、Rh又はRlの補正値を求める。
【0026】
この体厚補正は、従来のDXA法の骨密度測定装置でも行われているが、従来装置では、高エネルギーX線についての減衰量Rhに対して補正を加えることで、良好な体厚補正が実現されている。ただし、これは、骨密度測定に対して良好な体厚補正であり、本発明者の実験によれば、同じ補正を体脂肪率測定に適用しても良好な結果は得られなかった。そこで、実験により体脂肪率に適した体厚補正法を調べたところ、低エネルギーX線についての減衰量Rlに対して補正を加えるとが好適であることが分かった。この点については、厳密な理論的解明はまだであるが、大まかなメカニズムは次のようなものと考えられる。すなわち、骨によるX線の減衰は軟部組織による減衰よりもはるかに大きいため、骨の部分(すなわちX線の透過経路に骨が含まれる部分)では高エネルギーのX線でも十分に減衰するため、骨密度には高エネルギーX線の減衰量が支配的となる。このため、高エネルギーX線の減衰量を補正することで良好な体厚補正が可能となるものと考えられる。これに対し、軟部組織はX線の減衰が小さいので、軟部組織の部分(すなわちX線の透過経路に骨が含まれない部分)では体厚の違いによる減衰量の差は小さいものとなるが、そのなかでも高エネルギーX線の方が透過性が高いため、低エネルギーX線よりも体厚の違いによる差が現れにくい。また、軟部組織の脂肪と筋肉ではX線の減衰量の差が小さく、特に透過力の大きい高エネルギーのX線では体脂肪率が変わっても減衰量の変化は極めて小さい。これに対し、透過力がそれより小さい低エネルギーX線の減衰量は、脂肪と筋肉の減衰の差を高エネルギーX線よりも大きく反映する。したがって、軟部組織の部分を対象とする体脂肪率計算では、脂肪と筋肉の比率により敏感な低エネルギーX線の減衰量を補正することで、好適な体厚補正が可能になるものと考えられる。
【0027】
なお、体厚補正のための具体的な換算関数(あるいは換算テーブル)を求めるためには、まず骨密度用と体脂肪率用にそれぞれ上述のごとく厚みの異なる複数のファントムを用意し、それぞれ高低各エネルギーのX線で測定して、各体厚に対応する減衰量Rl,Rhを求める。そして、その測定結果において、減衰量比Rl/Rhがファントムの厚みによらずできるだけ均一となるよう、骨密度測定については高エネルギーX線の各厚みに対応する減衰量Rhの補正値Rh’を、体脂肪率測定については低エネルギーX線の各厚みに対応する減衰量Rlの補正値Rl’を、それぞれ求める。そして、骨密度測定については、各厚みに対するRhとRh’の対応関係を近似する関数を換算関数として決定し、体脂肪率測定については、各厚みに対するRlとRl’の対応関係を近似する関数を換算関数として決定する。このように決定した体厚補正用の換算関数は、記憶部30に予め記憶されている。
【0028】
したがってS14では、生データの各画素のRhに対して骨密度用の換算関数により体厚補正を施し、S22では生データの各画素のRlに対して体脂肪率用の換算関数により体厚補正を施した上で、それぞれ各画素ごとにその補正結果を用いて減衰量比Rl/Rhを計算する。
【0029】
S14で骨部用の補正を施した減衰量比のマップができあがると、次にこのマップの各画素を、当該画素の減衰量比の値に従って、骨部(X線透過経路に骨を含む部分)、軟部組織(X線透過経路に骨を含まないが軟部組織は含む部分)、及び空気層(X線透過経路が空気層のみの部分)の3種類のいずれかに分類する。骨部と軟部組織と空気層では、減衰量比Rl/Rhが大きく異なるので、しきい値弁別によりそれら3種類を弁別できる。この時の弁別しきい値としては、骨部と軟部組織とを分ける減衰量比と、軟部組織と空気層とを分ける減衰量比の値を予め実験等で求めておけばよい。この弁別によれば、各画素は、骨部、軟部組織、空気層のいずれか1つに必ず分類される。
【0030】
このS16により、S12で求めた減衰量のマップとS14及びS22で求めた減衰量比のマップにおいて、どの画素が骨部に該当し、どの画素が軟部組織に該当するかを示す分類マップが得られる。
【0031】
S18では、この分類マップを参照して骨部の画素のみを抽出し、それら骨部の各画素ごとに、S14で求めた減衰量比マップの中の当該画素の値を取り出し、この値に基づき周知のDXA法の手法に従って当該画素の骨密度を計算する。
【0032】
このようにして、骨部の各画素の骨密度が求められると、次に各画素ごとに、その画素の骨密度を所定の第1の色の階調値に変換することで、骨密度分布の画像を形成する(S20)。骨密度と階調値との関係を示す変換関数(又はテーブル)は、予め記憶部30に登録されている。
【0033】
以上で、骨密度分布の画像ができあがる。一方、体脂肪率測定の処理では、S22の処理の後、S16で求めた分類マップを参照して軟部組織の各画素を抽出し、それら軟部組織の各画素ごとに、S22で求めた減衰量比マップの中の当該画素の値を取り出し、この値に基づき、上述の体脂肪率計算方法に従って当該画素の体脂肪率を計算する(S24)。
【0034】
このようにして、軟部組織の各画素の体脂肪率が求められると、次に各画素ごとに、その画素の体脂肪率を、骨密度表示に割り当てた第1の色とは異なる所定の第2の色の階調値に変換することで、体脂肪率分布の画像を形成する(S26)。ここで、体脂肪率と階調値との関係を示す変換関数(又はテーブル)は、予め記憶部30に登録されている。
【0035】
そして、S28では、S20で生成した骨密度分布の画像と、S26で生成した体脂肪率分布の画像とを合成し、1枚の出力画像を生成する。この合成では、S16で作成した分類テーブルを参照し、骨部に該当する画素には骨密度分布における対応画素の値を採用し、軟部組織に該当する画素には体脂肪率分布における対応画素の値を採用する。
【0036】
なお、骨密度に割り当てる第1の色と、体脂肪率に割り当てる第2の色は、例えばそれぞれ表示出力機構の原色から選択すればよい。例えば、RGB表示のディスプレイ装置の場合、第1の色にG(緑)、第2の色にR(赤)を選択するなどである。また、骨密度は、骨との対応が直感的に分かりやすい白黒のグレースケールで表現し、体脂肪率は有彩色(例えば赤など)の階調で表現する等の方式も考えられる。もちろんこれらは一例であり、第1の色と第2の色には、階調を変えても両者の違いが区別できるような色の組合せであれば、どのようなものであっても用いることができる。
【0037】
そして、このようにして生成された合成画像を、例えば表示器32に表示出力したり、プリンタ(図示省略)から印刷出力したりする(S30)。
【0038】
図4は、本実施形態の装置から出力される、骨密度分布及び体脂肪率分布の合成画像の一例を模式的に示している。図では色が表現できないため、第1の色は斜線ハッチングで、第2の色はクロスハッチングでそれぞれ示し、骨密度や体脂肪率の大小はハッチングの密度で表現している。この図に示すように、骨密度分布と体脂肪率分布の合成画像200の近傍には、各色の階調と、骨密度又は体脂肪率の値(又は大小)との対応関係を表すスケール210a、210bが示される。なお、この図では、繁雑さを避けるため階調数を少なくしているが、実際の出力画像では階調を更に細かくすることができる。
【0039】
このように、本実施形態の装置によれば、被検体の骨部に該当する部分には骨密度を、軟部組織に該当する部分には体脂肪率を、それぞれ示した合成画像200を作成し、出力することができる。ユーザは、この合成画像200により骨密度の分布と体脂肪率の分布とを一覧することができるので、体脂肪率と骨密度との関係の把握が容易になる。また、本実施形態では、体脂肪率に適した演算方式及び体厚補正方式を用いているので、従来技術で説明した軟部組織表示とは異なり、信頼性の高い体脂肪率の情報を提供することができる。
【0040】
また、本実施形態では、骨密度と体脂肪率を、互いに異なる第1の色と第2の色の階調でそれぞれ表現したので、合成画像200上で、各画素が骨密度(骨部)を表しているか、体脂肪率(軟部組織)を表しているかが区別しやすい。
【0041】
なお、骨密度と体脂肪率を、互いに異なる第1の色と第2の色の階調で表現するという画像表現方式はあくまで一例である。原理的には、骨密度と体脂肪率とには、互いに重ならない表示形態範囲を割り当てればよい。この表示形態範囲には、上述した例のような1つの色の階調範囲も含まれる。この場合、例えば「第1の色」が1つの表示形態範囲であり、その第1の色の範囲内でのある階調の色が、その表示形態範囲内の1つの表示形態となる。また、寒色や暖色といった色相環上での範囲も表示形態範囲の一種である。すなわち、例えば骨密度には寒色範囲の色相のグラデーションを割り当て、体脂肪率には暖色範囲の色相のグラデーションを割り当てるなどの方式が可能である。また、図4に示したハッチングの区別なども表示形態範囲の一種であり、この場合1つ1つのハッチングの種類(左下がり斜線ハッチングや、クロスハッチング)が1つの表示形態範囲に該当し、あるハッチング種類におけるある密度のハッチングが、その表示形態範囲内の1つの表示形態に該当する。いずれの場合も、骨密度と体脂肪率にそれぞれ割り当てる表示形態範囲は、それら骨密度や体脂肪率が上限から下限までどのような値をとっても、合成画像200上で両者が区別しやすいものであることが望ましい。
【0042】
図1に示す実施形態ではX線発生装置18及びX線検出装置22が水平方向に駆動されているが、上述の説明から明らかなように、本実施形態における骨密度分布と体脂肪率分布の同時表示は、このような検出系の構成によらないものである。例えば、X線発生装置18及びX線検出装置22を固定配置し、ベッド16を水平方向に駆動してもよい。この場合、走査装置24によってベッド16が駆動され、これによって被検体10自体がスキャンされることになる。また、検出系としては、X線発生装置18から1点にフォーカスするペンシルビームを出力し、これを1個の検出器からなるX線検出装置22で検出してもよいし、X線発生装置18から扇型に広がるファンビームを出力し、これを複数の検出器が一列に並んだアレイを有するX線検出装置22で検出してもよい。また、X線発生装置18から円錐状に広がるコーンビームを用いる構成でもよい。また、上述の例では、高エネルギーのX線と低エネルギーのX線とを短い間隔で交互に照射したが、この代わりに広帯域のX線を照射し、X線検出装置22で検出信号の波高弁別を行って、高エネルギーX線に対応する検出信号と低エネルギーX線に対応する検出信号とを弁別して個別計数する構成でもよい。
【0043】
【発明の効果】
以上説明したように、本発明によれば、被検体骨部の骨密度分布と、被検体軟部組織の体脂肪率分布とを、1枚の画像に合成して表示することができるので、ユーザによる体脂肪率と骨密度との関係の把握が容易になる。
【図面の簡単な説明】
【図1】実施形態のX線骨密度測定装置の概略構成を示すブロック図である。
【図2】実施形態の体脂肪率測定の原理を説明するための図である。
【図3】実施形態の装置の処理の流れを説明するための図である。
【図4】実施形態の装置の出力画像の一例を模式的に示す図である。
【図5】DXA法による骨密度測定の原理を説明するための図である。
【符号の説明】
10 被検体、12 骨部、14 軟部組織、16 ベッド、18 X線発生装置、22 X線検出装置、24 走査装置、26 制御部、30 記憶部、32 表示器。
[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to an apparatus for measuring a bone density distribution of a subject using X-rays.
[0002]
[Prior art]
2. Description of the Related Art An X-ray bone density measuring apparatus that calculates bone density by transmitting X-rays through a subject has been conventionally used in many medical institutions for diagnosing and preventing bone diseases. As is well known, bones are composed of minerals (bone minerals) such as calcium that attenuate X-rays, and by calculating the amount of attenuation (attenuation rate) during X-ray transmission, the bone density per unit area is calculated. Is done. The bones are covered with soft tissues such as muscles and fats, and X-rays are attenuated in those areas as well. In order to eliminate such an influence due to soft tissue, a conventional bone density measuring apparatus uses X-rays of two types of energy. In this method, the simultaneous equations are solved based on the amounts of attenuation calculated for the respective energies, thereby specifying the amounts of attenuation of the bone portion and the soft tissue. This method is called a double X-ray absorption measurement method (DXA method). The principle of the bone density calculation by the DXA method will be described with reference to FIG. FIG. 5 shows the distribution of the attenuation when the X-rays of two kinds of high and low energies are transmitted through the subject 10 in which the soft tissue 14 covers the bone 12 as shown in FIG. 5B. (A). As can be seen from the figure, both the bone and the soft tissue attenuate more with low-energy X-rays than with high-energy X-rays, but the ratio of attenuation between the two energies is clearly different between bone and soft tissue. different. Therefore, in the DXA method, the ratio α of the attenuation of low-energy X-rays to the attenuation of high-energy X-rays in soft tissue is determined, and α is determined from the attenuation of low-energy X-rays to the attenuation of high-energy X-rays. By subtracting the product of the multiplication, the amount of attenuation of only the bone is obtained, and this is converted into bone density (in this case, the amount of bone mineral per unit area). Patent Documents 1 and 2 show examples of a bone density measuring device using such a method.
[0003]
In this type of bone density measuring device, for example, X-rays are irradiated from above to below a subject lying on a bed to detect the X-rays, thereby obtaining the bone density at each point when the subject is two-dimensionally projected. And the two-dimensional distribution can be displayed.
[0004]
Some conventional bone density measuring devices using the DXA method can display not only bones but also soft tissues for the purpose of easily grasping the external shape and the like of the subject. That is, in the general bone density calculation, the ratio α is used so that the value of the soft tissue becomes 0 as described above. However, if the value of α is set to a smaller value, the soft tissue also has a positive value. Therefore, if a display image is formed using such a small α, a display also showing a soft tissue can be obtained.
[0005]
[Patent Document 1]
Patent No. 2735507
[Patent Document 2]
Patent No. 3256667
[0006]
[Problems to be solved by the invention]
In recent years, an X-ray bone density measuring device has been used for confirming a drug effect using an animal such as a rat or a mouse. For example, there has been known an example in which a bone density distribution is measured by an X-ray bone density measuring device in order to examine the effect of changes in the amount of muscle and fat caused by a thinner drug on bone. In this application, it is very convenient to know the relationship between the distribution of the body fat percentage and the bone density distribution. However, there is no conventional bone density measurement device having a function of displaying such both distributions in association with each other.
[0007]
Here, the soft tissue display shown in the above-mentioned prior art can display not only the bone density distribution but also the soft tissue surrounding the bone to some extent. The attenuation was converted to a bone density scale, and did not represent the body fat percentage. For example, the body fat percentage generally tends to be higher at the side of the body (such as the flank), but since the thickness of the side of the body is small, the attenuation of X-rays is small. Is a small value equivalent to that of the outside of the subject, and does not represent the state of the body fat percentage.
[0008]
The present invention has been made in view of such a problem, and an object of the present invention is to provide an X-ray bone density measurement device that can clearly show the relationship between bone density and body fat percentage.
[0009]
[Means for Solving the Problems]
In order to achieve the above object, an apparatus according to the present invention detects X-rays transmitted through a subject, and a measuring unit that obtains a two-dimensional distribution of X-ray attenuation based on the detection result, Discriminating means for discriminating between a bone part and soft tissue of a subject; a bone density calculating means for obtaining a bone density distribution of the bone part from the two-dimensional distribution; and a body fat percentage distribution of the soft tissue from the two-dimensional distribution. It includes a body fat percentage calculating means to be obtained and a synthesizing means for synthesizing the bone density distribution of the bone part and the body fat percentage distribution of the soft tissue into one image and outputting the image.
[0010]
In this configuration, from the two-dimensional distribution of attenuation, a bone density distribution is obtained for a bone, and a body fat percentage distribution is obtained for a soft tissue, and the two are combined to generate one image. Therefore, the user can view the bone density distribution and the body fat percentage distribution on the basis of the image, which facilitates the analysis of the relationship between the bone density and the body fat percentage.
[0011]
In a preferred aspect, the measurement means obtains a two-dimensional distribution of the attenuation amount for each of two different X-ray energies, and the bone density calculation means obtains a two-dimensional distribution from the two-dimensional distribution of each X-ray of each energy. A bone density distribution is obtained according to a heavy X-ray absorption measurement method, and the body fat percentage distribution calculating means has information on a correspondence relationship between a ratio between attenuation amounts of X-rays of respective energies and a body fat percentage, and the two-dimensional information. The distribution of the body fat percentage is obtained by obtaining the ratio between the attenuation amounts of the X-rays of each energy at each point from the distribution and obtaining the body fat percentage corresponding to the ratio from the information of the correspondence.
[0012]
In another preferred aspect, the synthesizing unit assigns different display form ranges to the bone density distribution and the body fat percentage distribution, and sets the bone density and the body fat percentage in the corresponding display form ranges respectively. Express in display form. According to this aspect, it is possible to generate an image in which the bone density and the body fat percentage are easily distinguishable from each other.
[0013]
In still another preferred aspect, the body fat percentage distribution calculating means calculates a ratio between attenuation amounts of X-rays of respective energies at respective points from the two-dimensional distribution. The value of the attenuation of the two-dimensional distribution with respect to the lower energy X-ray is corrected. According to this aspect, it is possible to perform appropriate attenuation correction appropriate for the characteristics of the soft tissue.
[0014]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention (hereinafter, referred to as embodiments) will be described with reference to the drawings.
[0015]
The X-ray bone density measuring apparatus according to the present embodiment has a function of obtaining a bone density distribution and a body fat percentage distribution of a subject, synthesizing these into one image, and displaying them. This apparatus may have a hardware configuration substantially the same as a conventional bone density measuring apparatus using the DXA method, and the calculation of the body fat percentage distribution and the generation of a composite image are typically realized by software ( Of course, it is not impossible to implement a part or all of the calculation and image generation processing as a hardware circuit.)
[0016]
Therefore, first, a hardware configuration of an X-ray bone density measuring apparatus to which the simultaneous display method of the bone density and the body fat percentage according to the present invention is applied will be described with reference to FIG.
[0017]
This bone density measuring device is a measuring device for measuring a bone density distribution of a body of a subject 10 such as a rat, for example, and the subject 10 is mounted on a fixedly arranged bed 16. Actually, a mechanism for fixing the small animal to the bed 16 exists, but is not shown in the figure to avoid complication. The subject 10 is a living body, and the bone 12 is covered with a soft tissue 14 composed of muscle, fat, and the like. In the present embodiment, an X-ray generator 18 is provided below the bed 16, and an X-ray detector 22 is arranged above the bed 16. The X-ray generation device 18 is a device that generates high-energy X-rays or low-energy X-rays as necessary, and the X-ray detection device is a device that detects the X-rays 100 transmitted through the subject 10. The X-ray generation device 18 and the X-ray detection device 22 are driven in a horizontal direction integrally by a scanning device 24. For example, the whole body of the subject is scanned by zigzag scanning, and transmitted X-rays are detected.
[0018]
The control unit 26 controls the operation of the present apparatus, and executes various calculations described later. Processing such as calculation of bone density and body fat percentage and generation of a composite display image are also realized by the arithmetic processing of the control unit 26. The storage unit 30 stores various programs and data for control and arithmetic processing of the control unit 26. The storage unit 30 is composed of, for example, an EEPROM. The display 32 displays the result of calculation in the control unit 26 and the like.
[0019]
As shown in FIG. 1, when measuring the subject 10, the subject 10 is first placed on the bed 16, and then the X-ray generation device 18 and the X-ray detection device 22 are horizontally moved by the scanning device 24. , High-energy X-rays and low-energy X-rays are alternately irradiated at a predetermined cycle of, for example, about several milliseconds. By detecting this with the X-ray detector 22, the count rate (Ih) of the high-energy X-rays and the count rate (Il) of the low-energy X-rays transmitted through the subject 10 are determined for each point in the horizontal plane. Each is required. The amount of attenuation at each point is represented by the natural logarithm of the actually measured count rates Ih and Il divided by the reference count rates Ioh and Iol, respectively. That is, the attenuation Rh of the high energy X-ray is obtained by ln (Ioh / Ih), and the attenuation Rl of the low energy X-ray is obtained by ln (Iol / Il). Note that the reference count rates Ioh and Iol are obtained by irradiating high-energy X-rays and low-energy X-rays from the X-ray generator 18 to a position on the bed 16 where there is no subject under the same conditions as when the subject is actually measured. The count rate obtained from the detection result of the X-ray detection device 22 can be used. The attenuation of each of the high-energy X-rays and the low-energy X-rays obtained in this manner becomes raw data measured by this apparatus. Hereinafter, based on the raw data, calculation of bone density and body fat percentage at each point is performed.
[0020]
Here, the calculation of the bone density may be performed in accordance with the well-known algorithm of the DXA method, and a description thereof will be omitted. The calculation method of the body fat percentage will be described with reference to FIG.
[0021]
In order to calculate the body fat percentage, calibration is performed using a phantom that simulates soft tissue with respect to X-ray transmission, and an attenuation ratio-body fat percentage function indicating the relationship between the attenuation ratio and the body fat percentage is obtained. . For this calibration, phantoms corresponding to two different body fat percentages are used. For example, two types of phantoms (called phantom Fat0) corresponding to a body fat percentage of 0% (ie, muscle 100%) and phantoms (called a phantom Fat100) corresponding to a body fat percentage of 100% (muscle 0%) are used. . Then, a phantom Fat0 having a body fat percentage of 0% is placed on the bed 16 and irradiated with high-energy X-rays and low-energy X-rays, the respective attenuation amounts Rh and Rl are calculated, and the ratio Rl / Rh between the two is calculated. I do. This ratio is called an attenuation ratio. The attenuation ratio obtained for the phantom Fat0 is indicated as RR0 in FIG. FIG. 2 shows the relationship between the body fat percentage on the vertical axis and the attenuation ratio on the horizontal axis, and the point corresponding to the body fat percentage of 0% is indicated by A. Similarly, the phantom Fat100 is measured, and the attenuation ratio RR100 corresponding to the body fat percentage of 100% is obtained. The corresponding point is shown as point B in FIG. The range excluding the two points of the body fat percentage of 0% and 100% is approximated by a function using these two points as parameters. This function is called a conversion function f. The example of FIG. 2 is an example of a case where the two points are approximated by a linear conversion function f. Of course, a more appropriate conversion function may be found and used by experiments or the like. In any case, the calibration process may be performed using phantoms having different body fat percentages, which are necessary to specify the conversion function to be used.
[0022]
The conversion function f thus obtained is stored in the storage unit 30 and used. In the calculation of the body fat percentage, the attenuation ratio RRx is obtained by actually measuring the subject 100, and this RRx is substituted for the conversion function f to obtain the value F of the body fat percentage. The conversion function may be stored as a set of parameters defining the function form as described above, but may be implemented in a form such as a table in which values of the body fat percentage corresponding to each attenuation ratio are registered. Of course it is good.
[0023]
Next, with reference to FIG. 3, a procedure of processing performed by the apparatus according to the present invention will be described.
[0024]
First, the X-ray generator 18 and the X-ray detector 22 are scanned as described above to obtain attenuation amounts Rh and Rl of high-energy X-rays and low-energy X-rays at each point of the subject 100 (S10). Based on the measurement results, maps representing the two-dimensional distributions of Rh and Rl are created (S12). This map has the value of the amount of attenuation (Rh or Rl) for each pixel of the image. The set of these maps becomes the raw data of the measurement results. Thereafter, the processing of measuring the bone density and the processing of measuring the body fat percentage are executed based on the set of these maps.
[0025]
In both the bone density measurement and the body fat percentage measurement, the attenuation ratio RR = R1 / Rh is calculated for each pixel from the raw attenuation data (S14 and S20). However, the calculation of the attenuation amount ratio is not simply taking the ratio between the attenuation amounts of the high and low energies in the raw data, but in the case of the bone density measurement (S14) and in the case of the body fat percentage measurement (S20), respectively. Add an appropriate correction operation. This correction operation is performed to correct a deviation of a measurement result (bone density or body fat percentage) due to a difference in body thickness (also referred to as body thickness correction). Here, the body thickness is the length of the X-rays penetrating the subject, and ideally, even when the body thickness is different, when the same bone density (or body fat percentage) is measured, The same attenuation amount ratio RR should be obtained (that is, the relationship between the two is linearly proportional), but in reality this is not the case, and the relationship between the two becomes nonlinear. Therefore, after correcting the raw data of Rh and Rl so as to be as close to ideal as possible, the ratio between the two is determined. However, in the case of the bone density measuring device of the DXA method, the body thickness of each part of the subject cannot be measured due to the configuration of the device (a dedicated measuring mechanism may be provided for measuring the body thickness, but in view of cost, Therefore, a configuration is employed in which the attenuation of the raw data is simply converted to the attenuation corrected for the body thickness. That is, since the attenuations Rh and Rl are not proportional to the body thickness, but change monotonically, it can be assumed that the raw attenuation itself has a one-to-one correspondence with the body thickness. Such a conversion method is possible. For this conversion, a plurality of phantoms of the same material having different thicknesses are prepared, and they are measured with X-rays of high and low energies, so that the attenuation ratio Rl / Rh at that time is as uniform as possible irrespective of the body thickness. Next, a correction value of Rh or Rl is obtained.
[0026]
This body thickness correction is also performed by the conventional bone density measuring device of the DXA method. However, in the conventional device, by correcting the attenuation Rh for high-energy X-rays, good body thickness correction can be performed. Has been realized. However, this is a good correction of the body thickness for the measurement of bone density, and according to the experiments of the inventor, good results were not obtained even if the same correction was applied to the measurement of the body fat percentage. Therefore, when a body thickness correction method suitable for the body fat percentage was examined by an experiment, it was found that it is preferable to correct the attenuation Rl for low energy X-rays. Although a rigorous theoretical elucidation has not yet been made on this point, the rough mechanism is considered to be as follows. That is, since the attenuation of X-rays by bone is much greater than that of soft tissue, the high energy X-rays are sufficiently attenuated in the bone part (ie, the part where the bone is included in the X-ray transmission path). High energy X-ray attenuation is dominant in bone density. For this reason, it is considered that good body thickness correction can be performed by correcting the amount of attenuation of high-energy X-rays. On the other hand, since soft tissue has a small attenuation of X-rays, the difference in attenuation due to the difference in body thickness is small in the soft tissue portion (ie, the portion where bone is not included in the X-ray transmission path). Among them, since high energy X-rays have higher transparency, differences due to differences in body thickness are less likely to appear than low energy X-rays. Further, the difference between the attenuation of X-rays is small in fat and muscle of soft tissue, and the change in attenuation is extremely small even if the body fat percentage changes, especially for high-energy X-rays having high permeability. On the other hand, the attenuation of low-energy X-rays having a smaller penetrating power reflects the difference between attenuation of fat and muscle more greatly than that of high-energy X-rays. Therefore, in the body fat percentage calculation for the soft tissue portion, it is considered that a preferable body thickness correction becomes possible by correcting the attenuation amount of the low energy X-ray which is more sensitive to the ratio of fat and muscle. .
[0027]
In order to obtain a specific conversion function (or conversion table) for correcting the body thickness, first, a plurality of phantoms having different thicknesses as described above are prepared for the bone density and the body fat percentage. Measurement is performed using X-rays of each energy, and the attenuation amounts Rl and Rh corresponding to each body thickness are obtained. Then, in the measurement results, the correction value Rh ′ of the attenuation Rh corresponding to each thickness of the high-energy X-ray is measured for the bone density measurement so that the attenuation ratio Rl / Rh is as uniform as possible irrespective of the thickness of the phantom. For the measurement of the body fat percentage, a correction value Rl 'of the attenuation Rl corresponding to each thickness of the low energy X-ray is obtained. For bone density measurement, a function approximating the correspondence between Rh and Rh 'for each thickness is determined as a conversion function, and for body fat percentage measurement, a function approximating the correspondence between Rl and Rl' for each thickness. Is determined as a conversion function. The conversion function for body thickness correction determined in this way is stored in the storage unit 30 in advance.
[0028]
Therefore, in S14, the body thickness correction is performed on the Rh of each pixel of the raw data using the conversion function for bone density, and in S22, the body thickness correction is performed on the Rl of each pixel of the raw data using the conversion function for the body fat percentage. , And the attenuation ratio Rl / Rh is calculated for each pixel using the correction result.
[0029]
When the map of the attenuation ratio corrected for the bone portion is created in S14, each pixel of this map is then converted to the bone portion (the portion including the bone in the X-ray transmission path) according to the value of the attenuation ratio of the pixel. ), Soft tissue (portion not including bone in the X-ray transmission path but including soft tissue), and air layer (portion where the X-ray transmission path includes only the air layer). Since the attenuation ratio Rl / Rh greatly differs between the bone, the soft tissue, and the air layer, the three types can be discriminated by threshold discrimination. As the discrimination threshold value at this time, the value of the attenuation ratio that separates the bone and the soft tissue and the value of the attenuation ratio that separates the soft tissue and the air layer may be determined in advance by experiments or the like. According to this discrimination, each pixel is always classified into one of a bone, a soft tissue, and an air layer.
[0030]
By this S16, in the attenuation amount map obtained in S12 and the attenuation amount ratio maps obtained in S14 and S22, a classification map showing which pixels correspond to bones and which pixels correspond to soft tissues is obtained. Can be
[0031]
In S18, only the pixels of the bone are extracted with reference to this classification map, and for each pixel of the bone, the value of the pixel in the attenuation ratio map obtained in S14 is extracted, and based on this value, The bone density of the pixel is calculated according to the well-known DXA method.
[0032]
When the bone density of each pixel of the bone portion is obtained in this way, the bone density of each pixel is converted into a gradation value of a predetermined first color for each pixel. Is formed (S20). A conversion function (or table) indicating the relationship between the bone density and the gradation value is registered in the storage unit 30 in advance.
[0033]
Thus, an image of the bone density distribution is completed. On the other hand, in the process of measuring the body fat percentage, after the process of S22, each pixel of the soft tissue is extracted with reference to the classification map obtained in S16, and the attenuation amount obtained in S22 is extracted for each pixel of the soft tissue. The value of the pixel in the ratio map is extracted, and based on the value, the body fat percentage of the pixel is calculated according to the above-described body fat percentage calculation method (S24).
[0034]
When the body fat percentage of each pixel of the soft tissue is determined in this manner, the body fat percentage of the pixel is then determined for each pixel by a predetermined number different from the first color assigned to the bone density display. By converting the image into a gradation value of the second color, an image having a body fat percentage distribution is formed (S26). Here, the conversion function (or table) indicating the relationship between the body fat percentage and the gradation value is registered in the storage unit 30 in advance.
[0035]
Then, in S28, the image of the bone density distribution generated in S20 and the image of the body fat percentage distribution generated in S26 are combined to generate one output image. In this synthesis, the classification table created in S16 is referred to, the value of the corresponding pixel in the bone density distribution is adopted for the pixel corresponding to the bone, and the value of the corresponding pixel in the body fat percentage distribution is used for the pixel corresponding to the soft tissue. Adopt the value.
[0036]
Note that the first color assigned to the bone density and the second color assigned to the body fat percentage may be selected from, for example, the primary colors of the display output mechanism. For example, in the case of a display device of RGB display, G (green) is selected as the first color and R (red) is selected as the second color. In addition, a method of expressing the bone density in a black and white gray scale in which the correspondence with the bone is intuitive and easy to understand, and expressing the body fat percentage in a chromatic (for example, red) gradation can be considered. Of course, these are only examples, and any combination of colors can be used for the first color and the second color as long as they can be distinguished from each other even if the gradation is changed. Can be.
[0037]
Then, the synthesized image generated in this manner is displayed on, for example, the display 32 or printed out from a printer (not shown) (S30).
[0038]
FIG. 4 schematically shows an example of a composite image of the bone density distribution and the body fat percentage distribution, which is output from the device of the present embodiment. Since colors cannot be represented in the figure, the first color is indicated by hatching and the second color is indicated by cross hatching, and the magnitude of bone density and body fat percentage is indicated by hatching density. As shown in this figure, a scale 210a representing the correspondence between the gradation of each color and the value (or magnitude) of the bone density or the body fat percentage is provided near the composite image 200 of the bone density distribution and the body fat percentage distribution. , 210b are shown. Although the number of gradations is reduced in this figure to avoid complexity, the gradation can be further reduced in an actual output image.
[0039]
As described above, according to the apparatus of the present embodiment, a composite image 200 is created that shows the bone density in the portion corresponding to the bone of the subject and the body fat percentage in the portion corresponding to the soft tissue. , Can be output. Since the user can view the distribution of the bone density and the distribution of the body fat percentage from the composite image 200, the user can easily grasp the relationship between the body fat percentage and the bone density. Further, in the present embodiment, since the calculation method and the body thickness correction method suitable for the body fat percentage are used, unlike the soft tissue display described in the related art, highly reliable information on the body fat percentage is provided. be able to.
[0040]
Further, in the present embodiment, since the bone density and the body fat percentage are respectively expressed by gradations of the first color and the second color which are different from each other, on the composite image 200, each pixel has a bone density (bone portion). Or the body fat percentage (soft tissue).
[0041]
It should be noted that the image expression method of expressing the bone density and the body fat percentage by the gradations of the first color and the second color different from each other is merely an example. In principle, a display form range that does not overlap each other may be assigned to the bone density and the body fat percentage. This display form range includes the gradation range of one color as in the example described above. In this case, for example, “first color” is one display form range, and a color of a certain gradation within the first color range becomes one display form within the display form range. A range on the hue circle such as a cool color or a warm color is also a kind of display form range. That is, for example, a method of assigning a gradation of hues in a cool color range to bone density and assigning a gradation of hues in a warm color range to body fat percentage is possible. Also, the distinction of hatching shown in FIG. 4 is a kind of display form range, and in this case, each type of hatching (falling hatched hatch or cross hatching) corresponds to one display form range. Hatching of a certain density in the hatching type corresponds to one display form within the display form range. In any case, the display form range assigned to each of the bone density and the body fat percentage is such that the bone density and the body fat percentage can be easily distinguished on the composite image 200 regardless of the value from the upper limit to the lower limit. Desirably.
[0042]
In the embodiment shown in FIG. 1, the X-ray generator 18 and the X-ray detector 22 are driven in the horizontal direction. The simultaneous display does not depend on such a configuration of the detection system. For example, the X-ray generator 18 and the X-ray detector 22 may be fixedly arranged, and the bed 16 may be driven in the horizontal direction. In this case, the bed 16 is driven by the scanning device 24, thereby scanning the subject 10 itself. As a detection system, a pencil beam for focusing on one point may be output from the X-ray generator 18 and detected by the X-ray detector 22 including one detector. A fan beam spreading in a fan shape may be output from 18, and the fan beam may be detected by the X-ray detector 22 having an array in which a plurality of detectors are arranged in a line. Further, a configuration using a cone beam that spreads conically from the X-ray generator 18 may be used. In the above-described example, high-energy X-rays and low-energy X-rays are alternately irradiated at short intervals. Instead, a wide-band X-ray is irradiated. The discrimination may be performed, and the detection signal corresponding to the high energy X-ray and the detection signal corresponding to the low energy X-ray may be discriminated and individually counted.
[0043]
【The invention's effect】
As described above, according to the present invention, the bone density distribution of the subject's bone and the body fat percentage distribution of the subject's soft tissue can be combined into a single image and displayed. It becomes easy to grasp the relationship between the body fat percentage and the bone density.
[Brief description of the drawings]
FIG. 1 is a block diagram illustrating a schematic configuration of an X-ray bone density measuring device according to an embodiment.
FIG. 2 is a diagram for explaining the principle of measuring body fat percentage according to the embodiment.
FIG. 3 is a diagram for explaining a processing flow of the apparatus of the embodiment.
FIG. 4 is a diagram schematically illustrating an example of an output image of the apparatus according to the embodiment.
FIG. 5 is a diagram for explaining the principle of bone density measurement by the DXA method.
[Explanation of symbols]
DESCRIPTION OF REFERENCE NUMERALS 10 subjects, 12 bones, 14 soft tissues, 16 beds, 18 X-ray generators, 22 X-ray detectors, 24 scanners, 26 controllers, 30 storage units, 32 displays.

Claims (4)

被検体を透過したX線を検出し、この検出結果に基づきX線の減衰量の2次元分布を求める測定手段と、
前記2次元分布から被検体の骨部と軟部組織とを弁別する弁別手段と、
前記2次元分布から前記骨部の骨密度分布を求める骨密度算出手段と、
前記2次元分布から前記軟部組織の体脂肪率分布を求める体脂肪率算出手段と、
前記骨部の骨密度分布と前記軟部組織の体脂肪率分布とを一画像に合成して出力する合成手段と、
を含むX線骨密度測定装置。
Measuring means for detecting X-rays transmitted through the subject and obtaining a two-dimensional distribution of X-ray attenuation based on the detection result;
Discriminating means for discriminating a bone portion and a soft tissue of the subject from the two-dimensional distribution,
Bone density calculation means for obtaining a bone density distribution of the bone from the two-dimensional distribution,
Body fat percentage calculating means for obtaining a body fat percentage distribution of the soft tissue from the two-dimensional distribution,
A synthesizing unit that synthesizes and outputs a bone density distribution of the bone part and a body fat percentage distribution of the soft tissue into one image,
An X-ray bone density measurement device comprising:
前記測定手段は、異なる2つのエネルギーのX線についてそれぞれ前記減衰量の2次元分布を求め、
前記骨密度算出手段は、前記各エネルギーのX線についての前記2次元分布から二重X線吸収測定法に従って骨密度分布を求め、
前記体脂肪率分布算出手段は、各エネルギーのX線の減衰量同士の比と体脂肪率との対応関係の情報を有し、前記2次元分布から各点での各エネルギーのX線の減衰量同士の比を求め、その比に対応する体脂肪率を前記対応関係の情報から求めることで前記体脂肪率分布を求める、
ことを特徴とする請求項1記載のX線骨密度測定装置。
The measuring means obtains a two-dimensional distribution of the attenuation amount for each of two different energy X-rays,
The bone density calculation means obtains a bone density distribution from the two-dimensional distribution for the X-rays of each energy according to a double X-ray absorption measurement method,
The body fat percentage distribution calculating means has information on a correspondence relationship between a ratio between attenuation amounts of X-rays of each energy and a body fat percentage, and attenuates X-rays of each energy at each point from the two-dimensional distribution. Determine the ratio between the amounts, determine the body fat percentage distribution by determining the body fat percentage corresponding to the ratio from the information of the correspondence relationship,
The X-ray bone density measuring device according to claim 1, wherein:
前記合成手段は、前記骨密度分布と前記体脂肪率分布とに、互いに異なる表示形態範囲を割り当て、骨密度と体脂肪率とをそれぞれ対応する表示形態範囲内の表示形態で表現することを特徴とする請求項1又は2記載のX線骨密度測定装置。The synthesizing unit assigns different display form ranges to the bone density distribution and the body fat percentage distribution, and expresses the bone density and the body fat percentage in display forms within the corresponding display form ranges. The X-ray bone density measuring device according to claim 1 or 2, wherein 前記体脂肪率分布算出手段は、前記2次元分布から各点での各エネルギーのX線の減衰量同士の比を求めるに当たり、前記異なる2つのエネルギーのX線のうち、エネルギーが低い方のX線に対する前記2次元分布の減衰量の値を補正することを特徴とする請求項1〜3までのいずれか1項に記載のX線骨密度測定装置。The body fat percentage distribution calculating means calculates the ratio between the amounts of X-ray attenuation of each energy at each point from the two-dimensional distribution. The X-ray bone density measuring apparatus according to any one of claims 1 to 3, wherein a value of an attenuation amount of the two-dimensional distribution with respect to a line is corrected.
JP2002316490A 2002-10-30 2002-10-30 X-ray bone density measuring device Expired - Fee Related JP3974842B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002316490A JP3974842B2 (en) 2002-10-30 2002-10-30 X-ray bone density measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002316490A JP3974842B2 (en) 2002-10-30 2002-10-30 X-ray bone density measuring device

Publications (2)

Publication Number Publication Date
JP2004147863A true JP2004147863A (en) 2004-05-27
JP3974842B2 JP3974842B2 (en) 2007-09-12

Family

ID=32460184

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002316490A Expired - Fee Related JP3974842B2 (en) 2002-10-30 2002-10-30 X-ray bone density measuring device

Country Status (1)

Country Link
JP (1) JP3974842B2 (en)

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006312030A (en) * 2005-04-07 2006-11-16 Toshiba Corp X-ray ct apparatus
JP2007236766A (en) * 2006-03-10 2007-09-20 Shimadzu Corp X-ray photographing apparatus
JP2007333471A (en) * 2006-06-13 2007-12-27 Shimadzu Corp Nuclear medicine diagnosis device
JP2008119332A (en) * 2006-11-15 2008-05-29 Ge Medical Systems Global Technology Co Llc X-ray tomography apparatus
WO2009008102A1 (en) * 2007-07-11 2009-01-15 Hitachi Medical Corporation X ray ct system
JP2009178493A (en) * 2008-02-01 2009-08-13 Ge Medical Systems Global Technology Co Llc X-ray ct apparatus
JP2009213534A (en) * 2008-03-07 2009-09-24 Aloka Co Ltd X-ray ct apparatus and program
JP2011078528A (en) * 2009-10-06 2011-04-21 Aloka Co Ltd X-ray image pickup apparatus
JP2012040160A (en) * 2010-08-19 2012-03-01 Hitachi Aloka Medical Ltd X-ray ct apparatus
JP2012055451A (en) * 2010-09-08 2012-03-22 Hitachi Aloka Medical Ltd Bone salt quantity measuring instrument
WO2012128036A1 (en) * 2011-03-18 2012-09-27 日立アロカメディカル株式会社 Bone density measurement device
WO2012128031A1 (en) * 2011-03-23 2012-09-27 日立アロカメディカル株式会社 Bone diagnosis image display device
JP2012200329A (en) * 2011-03-24 2012-10-22 Hitachi Aloka Medical Ltd Bone mineral content measurement instrument
KR101249862B1 (en) 2011-05-13 2013-04-03 주식회사 비엠텍월드와이드 Bone density measuring apparatus with function of determining measuring position and function of analysing body composition using infrared camera
CN103054601A (en) * 2011-10-21 2013-04-24 三星电子株式会社 X-ray imaging apparatus and method for controlling same
JP2013106648A (en) * 2011-11-17 2013-06-06 Toshiba Corp Medical image diagnostic apparatus
KR101271255B1 (en) * 2009-06-12 2013-06-10 주식회사 엔씨테크론 Jaw Orthopedic Bone Mineral Density Measurement Device
JP2013236962A (en) * 2007-08-15 2013-11-28 Fujifilm Corp Image component separating device, method and program
JP2015019789A (en) * 2013-07-18 2015-02-02 日立アロカメディカル株式会社 X-ray measurement apparatus
WO2016047711A1 (en) * 2014-09-25 2016-03-31 日立アロカメディカル株式会社 Medical x-ray measuring device and method
WO2016190327A1 (en) * 2015-05-28 2016-12-01 株式会社日立製作所 Medical x-ray measurement system
JP2017093879A (en) * 2015-11-26 2017-06-01 株式会社日立製作所 X-ray measurement system and X-ray detection data processing method
JP2018134205A (en) * 2017-02-21 2018-08-30 富士フイルム株式会社 Energy subtraction processing apparatus, method, and program
JP2018153605A (en) * 2017-03-16 2018-10-04 富士フイルム株式会社 Body fat percentage measuring apparatus, method, and program
JP2018192052A (en) * 2017-05-18 2018-12-06 富士フイルム株式会社 Radiography system, radiography method, radiography program, and derivation apparatus
JP2018192051A (en) * 2017-05-18 2018-12-06 富士フイルム株式会社 Radiography system, radiography method, radiography program, and body thickness estimation apparatus
JP2019184450A (en) * 2018-04-12 2019-10-24 コニカミノルタ株式会社 X-ray imaging system
JP2019209028A (en) * 2018-06-08 2019-12-12 富士フイルム株式会社 Bone mineral information acquisition apparatus, method and program
JP2019209027A (en) * 2018-06-08 2019-12-12 富士フイルム株式会社 Bone mineral information acquisition apparatus, method and program
WO2020158659A1 (en) * 2019-02-01 2020-08-06 富士フイルム株式会社 Photographic subject information acquisition device, operating method therefor, and photographic subject information acquisition program
KR20200102028A (en) * 2019-02-20 2020-08-31 하이윈 테크놀로지스 코포레이션 Method for converting scan information of computed tomography scanner into bone parameters
US11234641B2 (en) 2017-03-16 2022-02-01 Fujifilm Corporation Body fat percentage measurement device, method and program
JP2022093509A (en) * 2018-04-12 2022-06-23 コニカミノルタ株式会社 X-ray imaging system

Cited By (59)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006312030A (en) * 2005-04-07 2006-11-16 Toshiba Corp X-ray ct apparatus
JP2007236766A (en) * 2006-03-10 2007-09-20 Shimadzu Corp X-ray photographing apparatus
JP2007333471A (en) * 2006-06-13 2007-12-27 Shimadzu Corp Nuclear medicine diagnosis device
JP4656008B2 (en) * 2006-06-13 2011-03-23 株式会社島津製作所 Nuclear medicine diagnostic equipment
JP2008119332A (en) * 2006-11-15 2008-05-29 Ge Medical Systems Global Technology Co Llc X-ray tomography apparatus
WO2009008102A1 (en) * 2007-07-11 2009-01-15 Hitachi Medical Corporation X ray ct system
CN101686824B (en) * 2007-07-11 2012-11-28 株式会社日立医药 X ray CT system
US7924969B2 (en) 2007-07-11 2011-04-12 Hitachi Medical Corporation X ray CT system
JP2013236962A (en) * 2007-08-15 2013-11-28 Fujifilm Corp Image component separating device, method and program
JP2009178493A (en) * 2008-02-01 2009-08-13 Ge Medical Systems Global Technology Co Llc X-ray ct apparatus
JP2009213534A (en) * 2008-03-07 2009-09-24 Aloka Co Ltd X-ray ct apparatus and program
KR101271255B1 (en) * 2009-06-12 2013-06-10 주식회사 엔씨테크론 Jaw Orthopedic Bone Mineral Density Measurement Device
JP2011078528A (en) * 2009-10-06 2011-04-21 Aloka Co Ltd X-ray image pickup apparatus
JP2012040160A (en) * 2010-08-19 2012-03-01 Hitachi Aloka Medical Ltd X-ray ct apparatus
JP2012055451A (en) * 2010-09-08 2012-03-22 Hitachi Aloka Medical Ltd Bone salt quantity measuring instrument
WO2012128036A1 (en) * 2011-03-18 2012-09-27 日立アロカメディカル株式会社 Bone density measurement device
TWI507175B (en) * 2011-03-18 2015-11-11 Hitachi Aloka Medical Ltd Bone-density measuring device
JP2012192118A (en) * 2011-03-18 2012-10-11 Hitachi Aloka Medical Ltd Bone density measuring apparatus
CN103442638A (en) * 2011-03-18 2013-12-11 日立阿洛卡医疗株式会社 Bone density measurement device
JP2012196401A (en) * 2011-03-23 2012-10-18 Hitachi Aloka Medical Ltd Bone diagnosis image display apparatus
CN103442640A (en) * 2011-03-23 2013-12-11 日立阿洛卡医疗株式会社 Bone diagnosis image display device
WO2012128031A1 (en) * 2011-03-23 2012-09-27 日立アロカメディカル株式会社 Bone diagnosis image display device
JP2012200329A (en) * 2011-03-24 2012-10-22 Hitachi Aloka Medical Ltd Bone mineral content measurement instrument
KR101249862B1 (en) 2011-05-13 2013-04-03 주식회사 비엠텍월드와이드 Bone density measuring apparatus with function of determining measuring position and function of analysing body composition using infrared camera
US8942347B2 (en) 2011-10-21 2015-01-27 Samsung Electronics Co., Ltd. X-ray imaging apparatus and method for controlling the same
CN108992081B (en) * 2011-10-21 2022-08-16 三星电子株式会社 X-ray imaging apparatus and control method thereof
JP2013085967A (en) * 2011-10-21 2013-05-13 Samsung Electronics Co Ltd X-ray imaging system and method for controlling the same
WO2013058567A1 (en) * 2011-10-21 2013-04-25 Samsung Electronics Co., Ltd. X-ray imaging apparatus and method for controlling the same
CN103054601A (en) * 2011-10-21 2013-04-24 三星电子株式会社 X-ray imaging apparatus and method for controlling same
CN108992081A (en) * 2011-10-21 2018-12-14 三星电子株式会社 X-ray imaging equipment and its control method
KR101337339B1 (en) 2011-10-21 2013-12-06 삼성전자주식회사 X-ray imaging apparatus and control method for the same
EP2583625B1 (en) * 2011-10-21 2017-03-22 Samsung Electronics Co., Ltd X-Ray imaging apparatus and method for controlling the same
JP2013106648A (en) * 2011-11-17 2013-06-06 Toshiba Corp Medical image diagnostic apparatus
JP2015019789A (en) * 2013-07-18 2015-02-02 日立アロカメディカル株式会社 X-ray measurement apparatus
WO2016047711A1 (en) * 2014-09-25 2016-03-31 日立アロカメディカル株式会社 Medical x-ray measuring device and method
JP2016064043A (en) * 2014-09-25 2016-04-28 日立アロカメディカル株式会社 Medical x-ray measuring apparatus and method
WO2016190327A1 (en) * 2015-05-28 2016-12-01 株式会社日立製作所 Medical x-ray measurement system
JP2016220850A (en) * 2015-05-28 2016-12-28 株式会社日立製作所 Medical X-ray Measurement System
JP2017093879A (en) * 2015-11-26 2017-06-01 株式会社日立製作所 X-ray measurement system and X-ray detection data processing method
JP2018134205A (en) * 2017-02-21 2018-08-30 富士フイルム株式会社 Energy subtraction processing apparatus, method, and program
JP2018153605A (en) * 2017-03-16 2018-10-04 富士フイルム株式会社 Body fat percentage measuring apparatus, method, and program
US11234641B2 (en) 2017-03-16 2022-02-01 Fujifilm Corporation Body fat percentage measurement device, method and program
CN108937983A (en) * 2017-05-18 2018-12-07 富士胶片株式会社 Radiation image picking-up system, radiation image photography method, program storage medium and guiding device
JP2018192051A (en) * 2017-05-18 2018-12-06 富士フイルム株式会社 Radiography system, radiography method, radiography program, and body thickness estimation apparatus
CN108937983B (en) * 2017-05-18 2023-04-28 富士胶片株式会社 Radiographic imaging system, radiographic imaging method, program storage medium, and deriving device
JP2018192052A (en) * 2017-05-18 2018-12-06 富士フイルム株式会社 Radiography system, radiography method, radiography program, and derivation apparatus
JP7306528B2 (en) 2018-04-12 2023-07-11 コニカミノルタ株式会社 Control device, display method and inspection method
JP2022093509A (en) * 2018-04-12 2022-06-23 コニカミノルタ株式会社 X-ray imaging system
JP7067221B2 (en) 2018-04-12 2022-05-16 コニカミノルタ株式会社 X-ray system
JP2019184450A (en) * 2018-04-12 2019-10-24 コニカミノルタ株式会社 X-ray imaging system
JP7016293B2 (en) 2018-06-08 2022-02-04 富士フイルム株式会社 Bone salt information acquisition device, method and program
US11096645B2 (en) 2018-06-08 2021-08-24 Fujifilm Corporation Bone mineral information acquisition apparatus, bone mineral information acquisition method, and bone mineral information acquisition program
JP7016294B2 (en) 2018-06-08 2022-02-04 富士フイルム株式会社 Bone salt information acquisition device, method and program
JP2019209027A (en) * 2018-06-08 2019-12-12 富士フイルム株式会社 Bone mineral information acquisition apparatus, method and program
JP2019209028A (en) * 2018-06-08 2019-12-12 富士フイルム株式会社 Bone mineral information acquisition apparatus, method and program
WO2020158659A1 (en) * 2019-02-01 2020-08-06 富士フイルム株式会社 Photographic subject information acquisition device, operating method therefor, and photographic subject information acquisition program
JP7284768B2 (en) 2019-02-01 2023-05-31 富士フイルム株式会社 SUBJECT INFORMATION ACQUISITION DEVICE, OPERATION METHOD THEREOF, AND SUBJECT INFORMATION ACQUISITION PROGRAM
KR102167859B1 (en) 2019-02-20 2020-10-21 하이윈 테크놀로지스 코포레이션 Method for converting scan information of computed tomography scanner into bone parameters
KR20200102028A (en) * 2019-02-20 2020-08-31 하이윈 테크놀로지스 코포레이션 Method for converting scan information of computed tomography scanner into bone parameters

Also Published As

Publication number Publication date
JP3974842B2 (en) 2007-09-12

Similar Documents

Publication Publication Date Title
JP3974842B2 (en) X-ray bone density measuring device
US11331068B2 (en) X-ray device, x-ray inspection method, and data processing apparatus
JP4545291B2 (en) Method and apparatus for controlling the dynamic range of digital diagnostic images
EP2548510B1 (en) Dual-energy X-ray imaging system and control method for the same
US20140086383A1 (en) Method and computerized tomography system for determining bone mineral density values
JPH10509075A (en) Bone densitometer with improved indicating characteristics
US20060098866A1 (en) User interface for inspection system with isoluminant regions
CN110770607B (en) Method and device for processing photon counting type X-ray detection data and X-ray device
CN104103055B (en) Automatically optimal output data is obtained
JP2009285356A (en) Image capturing system for medical use, image processing apparatus, image processing method, and program
US10426424B2 (en) System and method for generating and performing imaging protocol simulations
KR20130142654A (en) X-ray image apparatus and x-ray image system
KR101087137B1 (en) Method and support device for jawbone mineral density measurement
JPS62246352A (en) Digital medical diagnostic imaging apparatus and method
Rudolph et al. Influence of geometric distortion and exposure parameters on sensitivity of digital subtraction radiography
EP1253557A1 (en) Apparatus for and method of generating an enhanced contrast information digital image
Kim et al. Comparison between different cone-beam computed tomography devices in the detection of mechanically simulated peri-implant bone defects
Neitzel Management of pediatric radiation dose using Philips digital radiography
KR20140017468A (en) Dual-energy x-ray imaging apparatus and control method for the same
US20040017936A1 (en) Method, system and computer product for calculating mass scores
US5809104A (en) Method and device for measuring the content of bone mineral in the skeleton
JP2004275769A (en) Recording method by determining x-ray exposure value and x-ray diagnostic apparatus
Söderberg Automatic exposure control in CT an investigation between different manufacturers considering radiation dose and image quality
Önem et al. Effect of exposure parameters on the radiopacity of root canal sealers
Schaefer-Prokop et al. Computed radiography/digital radiography: radiologist perspective on controlling dose and study quality

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040727

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20070605

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20070612

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070615

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100622

Year of fee payment: 3

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100622

Year of fee payment: 3

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: 20120622

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130622

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

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

LAPS Cancellation because of no payment of annual fees