JP2017159032A - 骨密度計測装置及び方法 - Google Patents
骨密度計測装置及び方法 Download PDFInfo
- Publication number
- JP2017159032A JP2017159032A JP2017026832A JP2017026832A JP2017159032A JP 2017159032 A JP2017159032 A JP 2017159032A JP 2017026832 A JP2017026832 A JP 2017026832A JP 2017026832 A JP2017026832 A JP 2017026832A JP 2017159032 A JP2017159032 A JP 2017159032A
- Authority
- JP
- Japan
- Prior art keywords
- light
- bone density
- unit
- intensity
- skin
- 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
- 230000037182 bone density Effects 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims description 9
- 238000001739 density measurement Methods 0.000 title description 4
- 238000001514 detection method Methods 0.000 claims abstract description 25
- 238000002835 absorbance Methods 0.000 claims description 18
- 230000003287 optical effect Effects 0.000 claims description 16
- 230000001678 irradiating effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 2
- 210000000988 bone and bone Anatomy 0.000 abstract description 23
- 238000005259 measurement Methods 0.000 abstract description 19
- 238000000691 measurement method Methods 0.000 abstract description 3
- 241000283690 Bos taurus Species 0.000 description 9
- 210000001519 tissue Anatomy 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 4
- 229940028435 intralipid Drugs 0.000 description 4
- 238000004380 ashing Methods 0.000 description 3
- 229920001296 polysiloxane Polymers 0.000 description 3
- 210000000689 upper leg Anatomy 0.000 description 3
- VEXZGXHMUGYJMC-UHFFFAOYSA-N Hydrochloric acid Chemical compound Cl VEXZGXHMUGYJMC-UHFFFAOYSA-N 0.000 description 2
- 210000001185 bone marrow Anatomy 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 241000282465 Canis Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 208000001132 Osteoporosis Diseases 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000011088 calibration curve Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 239000006059 cover glass Substances 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000004590 silicone sealant Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
【課題】皮膚の下にある骨の骨密度の計測精度が高く、コンパクトなハンディタイプにもなる骨密度の計測装置及びそれを用いた骨密度の計測方法の提供を目的とする。【解決手段】皮膚の下側に存在する骨密度の計測装置であって、被検体に所定の強度の光を入射する光照射部と、前記被検体の皮膚から反射された散乱光を除去するフィルター手段と、前記散乱光が除去された反射光の集光手段と、当該集光手段にて集光された光の強度を計測する光検出部とを有することを特徴とする。【選択図】 図１
Description
本発明は骨密度を計測する装置及びそれを用いた計測方法に関し、特に光を用いた光学式骨密度計測に係る。
現在、骨粗鬆症診断にはＸ線を用いた装置が使用されているが、このような装置は大型であり、骨密度低下の早期発見等の簡易測定としては運用に難がある。
これに対して、近赤外光は生体内の情報を被破壊的に得る手段として、様々な医療分野で検討されている。
本発明者は、これまでに光学式の骨密度計測置として特許文献１，２を提案している。
本発明はさらに計測精度を向上させたものであり、特に骨の外側に存在する皮膚層の厚みや色等の生体組織の影響を抑えたものである。
これに対して、近赤外光は生体内の情報を被破壊的に得る手段として、様々な医療分野で検討されている。
本発明者は、これまでに光学式の骨密度計測置として特許文献１，２を提案している。
本発明はさらに計測精度を向上させたものであり、特に骨の外側に存在する皮膚層の厚みや色等の生体組織の影響を抑えたものである。
本発明は、皮膚の下にある骨の骨密度の計測精度が高く、コンパクトなハンディタイプにもなる骨密度の計測装置及びそれを用いた骨密度の計測方法の提供を目的とする。
本発明に係る骨密度計測装置は、皮膚の下側に存在する骨密度の計測装置であって、被検体に所定の強度の光を入射する光照射部と、前記被検体の皮膚から反射された散乱光を除去するフィルター手段と、前記散乱光が除去された反射光の集光手段と、当該集光手段にて集光された光の強度を計測する光検出部とを有することを特徴とする。
本発明は、光照射部から被検体に向けて光を照射すると、骨部にまで深く透過した光がはね返ってくる準直進光と、皮膚層にてはね返ってくる散乱光が存在することに着目し、この散乱光をフィルター手段にて除去した点に特徴がある。
またの準直進光は骨内にて減衰してはね返ってくるが、その強度分布は光の入射位置から半径方向に向かって減衰する。
この半径方向の減衰は骨密度が高い程、急激に減少することに着目したものである。
またの準直進光は骨内にて減衰してはね返ってくるが、その強度分布は光の入射位置から半径方向に向かって減衰する。
この半径方向の減衰は骨密度が高い程、急激に減少することに着目したものである。
本発明において、前記光照射部、光検出部、フィルター手段及び集光手段を組み込んだ光学ユニットを有し、前記被検体から前記光学ユニットの距離を移動制御出来る移動制御手段と、前記光照射部と光検出部とから吸光度を演算する吸光度演算手段とを有するようにすることもできる。
このようにすると、コンパクトな構造になり、ハンディタイプの計測装置となる。
このようにすると、コンパクトな構造になり、ハンディタイプの計測装置となる。
本発明において、前記集光手段は一対のレンズであり、前記フィルター手段はスリット板であってもよい。
このようにするとスリット板のスリット部からはね返ってくる準直線光を通過させつつ、スリット部以外の部分で皮膚層からはね返ってくる散乱光を除去することができる。
また、レンズは光検出部に向けて集光させるものであれば、各種レンズの組み合せが可能である。
このようにするとスリット板のスリット部からはね返ってくる準直線光を通過させつつ、スリット部以外の部分で皮膚層からはね返ってくる散乱光を除去することができる。
また、レンズは光検出部に向けて集光させるものであれば、各種レンズの組み合せが可能である。
本発明において皮膚から反射された散乱光を除去できるが、さらに皮膚の厚み情報を得るには、前記光照射部は波長の異なる複数の照射源を有し、前記光検出部は前記複数の照射源間での前記光の強度差の検出手段を有しているのが好ましい。
本発明に係る骨密度計測方法は、請求項１〜３のいずれかの骨密度計測装置を用いた骨密度計測方法であって、前記光照射部と被検体からの距離を変化させ、前記光照射部から照射された入射光の入射強度に対する前記光検出部にて検出された検出強度の変化に基づいて骨密度を計測することを特徴とする。
これにて得られる光検出部にて検出された検出強度の変化は、例えば光が照射された位置から半径方向の減衰の傾きが骨密度の高低により変化するので、予め骨密度が分かっている各種サンプルの計測値の検量線を用いて被検体の骨密度を求めることができる。
これにて得られる光検出部にて検出された検出強度の変化は、例えば光が照射された位置から半径方向の減衰の傾きが骨密度の高低により変化するので、予め骨密度が分かっている各種サンプルの計測値の検量線を用いて被検体の骨密度を求めることができる。
また、請求項４に記載の骨密度計測装置を用いた骨密度計測方法であって、前記光照射部と被検体からの距離を変化させ、前記光照射部から照射された入射光の入射強度に対する前記光検出部にて検出された検出強度から求められた光強度分布と前記光の強度差の検出手段にて求められた皮膚厚情報とに基づいて骨密度を計測することができる。
本発明は、皮膚組織層による影響をスリット板で除去しつつ、スリット板及びレンズを被検体に向けて移動させつつ、吸光度の強度分布を計測できるようにできたので、ハンディタイプのコンパクトな計測装置となる。
また、例えば近赤外光と可視光等、波長の異なる複数の照射源を用いて、その光強度差から皮膚の厚み情報を検出することで、さらに骨密度の計測精度が向上する。
また、例えば近赤外光と可視光等、波長の異なる複数の照射源を用いて、その光強度差から皮膚の厚み情報を検出することで、さらに骨密度の計測精度が向上する。
本発明に係る骨密度の計測装置の構成例及び計測方法例を以下説明するが、本発明はこれに限定されない。
図１に光学ユニットの構成例を示し、図２に計測原理の説明図を示す。
被検体１に向けて光を照射する光照射部１３と、被検体から返ってくる光から散乱光によるノイズを除去するためのフィルター手段として第１スリット板１１ａと第２スリット板１１ｂとを有し、第１スリット板１１ａのスリットを通過した光は平凸レンズからなる第１レンズ１２ａにて進路を略平行光に屈折させ、第２スリット板１１ｂのスリットを通過した光を凸レンズからなる第２レンズにて光検出部に向けて集光させ、この光検出部で光の強度を検出する。
ここで、光照射部から入射させる光の強度を入射光強度Ｉ０、光検出部にて検出された光の強度を検出光強度Ｉとすると、吸光度Ａは下記式（１）にて求められる。
被検体１に向けて光を照射する光照射部１３と、被検体から返ってくる光から散乱光によるノイズを除去するためのフィルター手段として第１スリット板１１ａと第２スリット板１１ｂとを有し、第１スリット板１１ａのスリットを通過した光は平凸レンズからなる第１レンズ１２ａにて進路を略平行光に屈折させ、第２スリット板１１ｂのスリットを通過した光を凸レンズからなる第２レンズにて光検出部に向けて集光させ、この光検出部で光の強度を検出する。
ここで、光照射部から入射させる光の強度を入射光強度Ｉ０、光検出部にて検出された光の強度を検出光強度Ｉとすると、吸光度Ａは下記式（１）にて求められる。
骨密度の計測原理を図２に基づいて説明する。
被検体に入射された光は、Ｓｋｉｎ（皮膚組織）を透過してＢｏｎｅ（骨）内まで通過し、返ってくる準直進光と皮膚組織層にて返ってくる散乱光を有する。
スリット板を設けることで図２（ａ）〜（ｃ）に示すようにスリット板のスリットから準直進光を通過させ、スリット板のスリット以外の部分にて皮膚組織層からの散乱光を除去することができる。
そこで、図２（ｄ）に示すようにスリット及びレンズの被検体からの距離Ｚを移動させると、光の入射位置から半径方向に減衰しながら返ってくる準直進光の強度分布を検出することができる。
ここで半径方向の減衰は骨密度が高い程、その傾きが大きいことから、その傾きの大きさを指標にして骨密度を計測することができる。
被検体に入射された光は、Ｓｋｉｎ（皮膚組織）を透過してＢｏｎｅ（骨）内まで通過し、返ってくる準直進光と皮膚組織層にて返ってくる散乱光を有する。
スリット板を設けることで図２（ａ）〜（ｃ）に示すようにスリット板のスリットから準直進光を通過させ、スリット板のスリット以外の部分にて皮膚組織層からの散乱光を除去することができる。
そこで、図２（ｄ）に示すようにスリット及びレンズの被検体からの距離Ｚを移動させると、光の入射位置から半径方向に減衰しながら返ってくる準直進光の強度分布を検出することができる。
ここで半径方向の減衰は骨密度が高い程、その傾きが大きいことから、その傾きの大きさを指標にして骨密度を計測することができる。
具体的に骨密度の計測に用いた骨密度計測装置２０の構成例を図３に示す。
ガンタイプとしたハンディ型の筐体に光学ユニット１０を組み込んだ計測装置２０となっている。
光学ユニット１０は光照射部１３にレーザーダイオードを配置し、その先端側に第１レンズ１２ａ及び第１スリット板１１ａを設け、このレーザーダイオードの後方に第２スリット板１１ｂと第２レンズ１２ｂを設け、集光された光がフォトダイオードからなる光検出部１４にて強度計測される。
光学ユニット１０は、アクチュエーター２１に取り付けられ、ステップモーター２２にて前進及び後退制御されている。
光検出部にて検出された信号の検出信号増幅器２３とステップモーターの制御部２４及びマイクロコンピューター２５等が組み込まれている。
ハンディタイプのグリップ部下側には、バッテリー２６が内蔵されている。
光学ユニットの移動距離データ、検出光強度データ等がパソコン２７にシリアル通信により取り込まれる。
ガンタイプとしたハンディ型の筐体に光学ユニット１０を組み込んだ計測装置２０となっている。
光学ユニット１０は光照射部１３にレーザーダイオードを配置し、その先端側に第１レンズ１２ａ及び第１スリット板１１ａを設け、このレーザーダイオードの後方に第２スリット板１１ｂと第２レンズ１２ｂを設け、集光された光がフォトダイオードからなる光検出部１４にて強度計測される。
光学ユニット１０は、アクチュエーター２１に取り付けられ、ステップモーター２２にて前進及び後退制御されている。
光検出部にて検出された信号の検出信号増幅器２３とステップモーターの制御部２４及びマイクロコンピューター２５等が組み込まれている。
ハンディタイプのグリップ部下側には、バッテリー２６が内蔵されている。
光学ユニットの移動距離データ、検出光強度データ等がパソコン２７にシリアル通信により取り込まれる。
灰化ウシ海綿骨をファントムとして用いた検証実験について以下説明する。
灰化ウシ海綿骨はウシ大腿骨から海綿骨を１辺、約４ｃｍのブロック状に切り出し、骨髄を煮沸除去し、６００℃で２４時間灰化処理したものを用いた。
試験体としては、灰化ウシ海綿骨を塩酸で溶解し、密度２５３ｍｇ／ｃｍ３、１７８ｍｇ／ｃｍ３、１５１ｍｇ／ｃｍ３のものを用いた。
光照射部の光源としてEdmond Optics社製のレーザーダイオード（＃５７−１０１）を用いた。
これは出力波長６５５ｎｍ、最大出力３ｍＷである。
吸光度Ａの測定として、レンズの集点距離となる被検体からの位置をＺ＝０として被検体の深部方向を正方向として計測した。
図４（ａ）は密度２５８ｍｇ／ｃｍ３ファントムを用いて異なる位置９点の移動量Ｚと吸光度Ａの計測グラフを示す。
各位置の吸光度分布は同様の変化を示した。
図４（ｂ）は密度２５８ｍｇ／ｃｍ３、１７８ｍｇ／ｃｍ３、１５１ｍｇ／ｃｍ３それぞれの９点測定平均値の推移グラフを示す。
このうちＺ＝１０ｍｍ以上の傾きに着目し、移動量Ｚ＝１１〜１８ｍｎの吸光度の減衰直線を図５（ａ）に示す。
ファントムの密度と図５（ａ）に示した減衰直線の傾きの関係を図５（ｂ）に示す。
これにより、吸光度の減衰の傾きを指標にして骨密度の計測が可能であることが分かる。
灰化ウシ海綿骨はウシ大腿骨から海綿骨を１辺、約４ｃｍのブロック状に切り出し、骨髄を煮沸除去し、６００℃で２４時間灰化処理したものを用いた。
試験体としては、灰化ウシ海綿骨を塩酸で溶解し、密度２５３ｍｇ／ｃｍ３、１７８ｍｇ／ｃｍ３、１５１ｍｇ／ｃｍ３のものを用いた。
光照射部の光源としてEdmond Optics社製のレーザーダイオード（＃５７−１０１）を用いた。
これは出力波長６５５ｎｍ、最大出力３ｍＷである。
吸光度Ａの測定として、レンズの集点距離となる被検体からの位置をＺ＝０として被検体の深部方向を正方向として計測した。
図４（ａ）は密度２５８ｍｇ／ｃｍ３ファントムを用いて異なる位置９点の移動量Ｚと吸光度Ａの計測グラフを示す。
各位置の吸光度分布は同様の変化を示した。
図４（ｂ）は密度２５８ｍｇ／ｃｍ３、１７８ｍｇ／ｃｍ３、１５１ｍｇ／ｃｍ３それぞれの９点測定平均値の推移グラフを示す。
このうちＺ＝１０ｍｍ以上の傾きに着目し、移動量Ｚ＝１１〜１８ｍｎの吸光度の減衰直線を図５（ａ）に示す。
ファントムの密度と図５（ａ）に示した減衰直線の傾きの関係を図５（ｂ）に示す。
これにより、吸光度の減衰の傾きを指標にして骨密度の計測が可能であることが分かる。
次に光照射部の光源にImatronic社製のレーザーダイオード（ＬＤＭ１１５Ｇ／８５０／１）を用いた。
これは波長８５０ｎｍ、最大出力１ｍＷである。
実施例１と同様に計測した移動量Ｚと吸光度の計測グラフを図６に示す。
実施例２では皮膚層の厚みの影響を調査すべく、その厚み０，０．２，０．５，１．０，２．０ｍｍの被検体を用いた。
皮膚層としては、２％イントラリピッド液を２枚のカバーガラスの間に封入し、その厚みを調整した模擬皮膚を用いた。
なお、２％イントラリピッドは、ヒトの皮膚と光学特性が等しいと報告されている（Troy,S.,et al.,Journal of biomedical optics,vol.6,pp.167-176,2001）。
Ｚ＝１８〜２５ｍｍの間の傾きを直線近似した傾きと骨密度の関係を図７に示す。
このグラフから皮膚組織が存在する場合に、その影響があることが推測された。
これは波長８５０ｎｍ、最大出力１ｍＷである。
実施例１と同様に計測した移動量Ｚと吸光度の計測グラフを図６に示す。
実施例２では皮膚層の厚みの影響を調査すべく、その厚み０，０．２，０．５，１．０，２．０ｍｍの被検体を用いた。
皮膚層としては、２％イントラリピッド液を２枚のカバーガラスの間に封入し、その厚みを調整した模擬皮膚を用いた。
なお、２％イントラリピッドは、ヒトの皮膚と光学特性が等しいと報告されている（Troy,S.,et al.,Journal of biomedical optics,vol.6,pp.167-176,2001）。
Ｚ＝１８〜２５ｍｍの間の傾きを直線近似した傾きと骨密度の関係を図７に示す。
このグラフから皮膚組織が存在する場合に、その影響があることが推測された。
次に図８に構造を模式的に示すように、骨密度計測用の照射源（光照射部）１３の他に皮膚の厚み情報を得るために、波長の異なる２つの照射源１３ａ，１３ｂを取り付け評価した。
具体的には、二つの円形スリット（スリット径：１０ｍｍ，スリット幅：２ｍｍ，厚さ５ｍｍ），二つの片凸レンズ（Edmund Optics, 48795, 及び48797），三つのレーザーダイオード［密度計測用（中央）：Egismos, H838501D（850nm，1mW）， 皮膚厚計測用：S638501D（850nm，1mW），H635151R（515nm，1mW）］、およびＰＤ（Hamamatsu, C12703-01）を１つの光学ユニットとしてまとめ、それをステッピングモーター（Orientalmoter, PKP225）により駆動されるアクチュエーター(Misumi, LX2005P-MX-B1-T2028-150)でＺ方向に移動できるようにした。
ＰＤより得られた検出信号は２次のアンチエイリアジング・フィルター回路を介してＡＤＣ素子（Microchip，MCP3208，12bit）によりＡＤＣ変換した。
モータやレーザー等の制御は、シングルボードコンピューター（Raspberry Pi，RS Components，Raspberry Pi 2 Model B）で行い、検出結果はＣＳＶテキストファイルで保存され、このデータに基づいて移動距離と吸光度の関係のグラフ、すなわち光強度分布がディスプレイ上に表示される。
また、実験に灰化ウシ海綿骨および灰化ウシ緻密骨、シリコーン（セメダイン株，シリコーンシーラント）を模擬皮膚として用いて実験を行った。
灰化ウシ海綿骨は、ウシ大腿骨から海綿骨を一辺４ｃｍほどのブロック状に切り出し、それを煮沸することで骨髄を除去、その後、電気炉内において６００℃で２４時間灰化処理することで９つのサンプルを作製した。
灰化ウシ海綿骨は場所により生体特有の密度のばらつきが存在する。
そのため、μＣＴ（Shimadzu，inspeXio SMX-90CT Plus）でスキャンし、各サンプルに対して異なる３つの位置で一辺１ｃｍの立方体範囲を選択し、骨形態計測用ソフト（ラトックシステムエンジニアリング株，TRI/3D-BON-FCS64）を用いて同範囲の骨密度を算出した。
緻密骨部は、ウシ大腿骨の骨梁部分より軸方向と平行に厚さ０．３ｍｍでスライスしたものを海綿骨と同様の方法で灰化処理を行ったものを用いた。
緻密骨厚さは、アメリカ国立医学図書館のデータベースであるVisible Human FTP Resourceの切片画像から橈骨遠位端の皮質骨厚さを読み取り、この値を参考に決定した。
また、皮膚の光学特性を模擬するために２％イントラリピッドの液体試料がよく用いられるが、本実験ではより扱い易い固体試料をシリコーンにより作製した。
なお、ホワイトとクリアを１：６で混ぜ合わせたシリコーンは２％イントラリピッドと同様の吸光度特性となることが確認されたことから、これを用い厚さ１〜２ｍｍの異なる厚さを持つ模擬皮膚を作製した。
これらの模擬皮膚，緻密骨，海綿骨を順に並べ、実験試料とした。
８５０ｎｍと５１５ｎｍの二波長を図８のように配置し、前述の模擬試料を用いて実験を行った。
図１２（ａ）は、８５０ｎｍならびに５１５ｎｍレーザーから得られる光強度分布のピークの位置を、Ｚを用いて表している。
なお、各プロットの値は１５回の計測の平均を示している。
同図の値において、８５０ｎｍから５１５ｎｍの値を引いたもの、すなわちピーク間距離をδとし、それと皮膚厚との関係を示したものが図１２（ｂ）である。
δは皮膚厚に対し負の相関を示し、この関係は、各波長の皮膚層での到達深度と光散乱特性の違いから生じたものだと考えられる。
一般的に、８５０ｎｍの光は５１５ｎｍのものに比べ皮膚内の奥深くまで侵入するため、皮膚厚に対してピークのシフトが大きく、逆に５１５ｎｍの光では、皮膚厚に対してあまり影響を受けない。
波長による皮膚厚の影響を図１０，図１１に示した。
そのため、各波長によるピークの差であるδが皮膚厚に対して変化したと考えられる。
δ，ｓｌｏｐｅ，およびＢＭＤの関係を示す近似平面を求めると次式のようになる。
なお、近似平面は、Python(Python 3.5)のＰｙＳｔａｎモジュールを用いてモンテカルロ法により求めた．
この近似式を使い、計測されるＳｌｏｐｅとδからＢＭＤを予測する式を導くと次のようになる．
ここで、ｓｄｐＢＭＤ(slope and delta predicted bone mineral density)は、Ｓｌｏｐｅとδから予測されるＢＭＤである。
図１３は、式（３）より算出された予測骨密度とμＣＴによる計測骨密度の関係を示している。
Ｓｌｏｐｅとδより予測された骨密度は、μＣＴで得られる骨密度とほぼ同様の値をとり、良好な正の相関（ｒ２＝０．７２５８１）を示した。
このことから、異なる波長を用いてＺ方向に移動させたその強度差から皮膚層の厚み情報を得ることで、さらに精度の高い骨密度を計測することができる。
具体的には、二つの円形スリット（スリット径：１０ｍｍ，スリット幅：２ｍｍ，厚さ５ｍｍ），二つの片凸レンズ（Edmund Optics, 48795, 及び48797），三つのレーザーダイオード［密度計測用（中央）：Egismos, H838501D（850nm，1mW）， 皮膚厚計測用：S638501D（850nm，1mW），H635151R（515nm，1mW）］、およびＰＤ（Hamamatsu, C12703-01）を１つの光学ユニットとしてまとめ、それをステッピングモーター（Orientalmoter, PKP225）により駆動されるアクチュエーター(Misumi, LX2005P-MX-B1-T2028-150)でＺ方向に移動できるようにした。
ＰＤより得られた検出信号は２次のアンチエイリアジング・フィルター回路を介してＡＤＣ素子（Microchip，MCP3208，12bit）によりＡＤＣ変換した。
モータやレーザー等の制御は、シングルボードコンピューター（Raspberry Pi，RS Components，Raspberry Pi 2 Model B）で行い、検出結果はＣＳＶテキストファイルで保存され、このデータに基づいて移動距離と吸光度の関係のグラフ、すなわち光強度分布がディスプレイ上に表示される。
また、実験に灰化ウシ海綿骨および灰化ウシ緻密骨、シリコーン（セメダイン株，シリコーンシーラント）を模擬皮膚として用いて実験を行った。
灰化ウシ海綿骨は、ウシ大腿骨から海綿骨を一辺４ｃｍほどのブロック状に切り出し、それを煮沸することで骨髄を除去、その後、電気炉内において６００℃で２４時間灰化処理することで９つのサンプルを作製した。
灰化ウシ海綿骨は場所により生体特有の密度のばらつきが存在する。
そのため、μＣＴ（Shimadzu，inspeXio SMX-90CT Plus）でスキャンし、各サンプルに対して異なる３つの位置で一辺１ｃｍの立方体範囲を選択し、骨形態計測用ソフト（ラトックシステムエンジニアリング株，TRI/3D-BON-FCS64）を用いて同範囲の骨密度を算出した。
緻密骨部は、ウシ大腿骨の骨梁部分より軸方向と平行に厚さ０．３ｍｍでスライスしたものを海綿骨と同様の方法で灰化処理を行ったものを用いた。
緻密骨厚さは、アメリカ国立医学図書館のデータベースであるVisible Human FTP Resourceの切片画像から橈骨遠位端の皮質骨厚さを読み取り、この値を参考に決定した。
また、皮膚の光学特性を模擬するために２％イントラリピッドの液体試料がよく用いられるが、本実験ではより扱い易い固体試料をシリコーンにより作製した。
なお、ホワイトとクリアを１：６で混ぜ合わせたシリコーンは２％イントラリピッドと同様の吸光度特性となることが確認されたことから、これを用い厚さ１〜２ｍｍの異なる厚さを持つ模擬皮膚を作製した。
これらの模擬皮膚，緻密骨，海綿骨を順に並べ、実験試料とした。
８５０ｎｍと５１５ｎｍの二波長を図８のように配置し、前述の模擬試料を用いて実験を行った。
図１２（ａ）は、８５０ｎｍならびに５１５ｎｍレーザーから得られる光強度分布のピークの位置を、Ｚを用いて表している。
なお、各プロットの値は１５回の計測の平均を示している。
同図の値において、８５０ｎｍから５１５ｎｍの値を引いたもの、すなわちピーク間距離をδとし、それと皮膚厚との関係を示したものが図１２（ｂ）である。
δは皮膚厚に対し負の相関を示し、この関係は、各波長の皮膚層での到達深度と光散乱特性の違いから生じたものだと考えられる。
一般的に、８５０ｎｍの光は５１５ｎｍのものに比べ皮膚内の奥深くまで侵入するため、皮膚厚に対してピークのシフトが大きく、逆に５１５ｎｍの光では、皮膚厚に対してあまり影響を受けない。
波長による皮膚厚の影響を図１０，図１１に示した。
そのため、各波長によるピークの差であるδが皮膚厚に対して変化したと考えられる。
δ，ｓｌｏｐｅ，およびＢＭＤの関係を示す近似平面を求めると次式のようになる。
なお、近似平面は、Python(Python 3.5)のＰｙＳｔａｎモジュールを用いてモンテカルロ法により求めた．
図１３は、式（３）より算出された予測骨密度とμＣＴによる計測骨密度の関係を示している。
Ｓｌｏｐｅとδより予測された骨密度は、μＣＴで得られる骨密度とほぼ同様の値をとり、良好な正の相関（ｒ２＝０．７２５８１）を示した。
このことから、異なる波長を用いてＺ方向に移動させたその強度差から皮膚層の厚み情報を得ることで、さらに精度の高い骨密度を計測することができる。
１ 被検体
１０ 光学ユニット
１１ａ 第１スリット板
１１ｂ 第２スリット板
１２ａ 第１レンズ
１２ｂ 第２レンズ
１３ 光照射部
１４ 光検出部
２０ 装置
２１ アクチュエーター
１０ 光学ユニット
１１ａ 第１スリット板
１１ｂ 第２スリット板
１２ａ 第１レンズ
１２ｂ 第２レンズ
１３ 光照射部
１４ 光検出部
２０ 装置
２１ アクチュエーター
Claims (6)
- 皮膚の下側に存在する骨密度の計測装置であって、
被検体に所定の強度の光を入射する光照射部と、
前記被検体の皮膚から反射された散乱光を除去するフィルター手段と、前記散乱光が除去された反射光の集光手段と、当該集光手段にて集光された光の強度を計測する光検出部とを有することを特徴とする骨密度計測装置。 - 前記光照射部、光検出部、フィルター手段及び集光手段を組み込んだ光学ユニットを有し、前記被検体から前記光学ユニットの距離を移動制御出来る移動制御手段と、前記光照射部と光検出部とから吸光度を演算する吸光度演算手段とを有することを特徴とする請求項１記載の骨密度計測装置。
- 前記集光手段は一対のレンズであり、前記フィルター手段はスリット板であることを特徴とする請求項１又は２記載の骨密度計測装置。
- 前記光照射部は波長の異なる複数の照射源を有し、
前記光検出部は前記複数の照射源間での前記光の強度差の検出手段を有していることを特徴とする請求項１〜３のいずれかに記載の骨密度計測装置。 - 請求項１〜３のいずれかに記載の骨密度計測装置を用いた骨密度計測方法であって、
前記光照射部と被検体からの距離を変化させ、前記光照射部から照射された入射光の入射強度に対する前記光検出部にて検出された検出強度から求められた光強度分布に基づいて骨密度を計測することを特徴とする骨密度計測方法。 - 請求項４に記載の骨密度計測装置を用いた骨密度計測方法であって、
前記光照射部と被検体からの距離を変化させ、前記光照射部から照射された入射光の入射強度に対する前記光検出部にて検出された検出強度から求められた光強度分布と前記光の強度差の検出手段にて求められた皮膚厚情報とに基づいて骨密度を計測することを特徴とする骨密度計測方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016039691 | 2016-03-02 | ||
| JP2016039691 | 2016-03-02 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2017159032A true JP2017159032A (ja) | 2017-09-14 |
| JP6951730B2 JP6951730B2 (ja) | 2021-10-20 |
Family
ID=59854498
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017026832A Active JP6951730B2 (ja) | 2016-03-02 | 2017-02-16 | 骨密度計測装置及び方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6951730B2 (ja) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1043163A (ja) * | 1996-03-26 | 1998-02-17 | Hologic Inc | 骨特性評価システム |
| US20060192965A1 (en) * | 2005-01-21 | 2006-08-31 | Tromberg Bruce J | Method for assessing the condition of bone in-vivo |
| WO2006115207A1 (ja) * | 2005-04-22 | 2006-11-02 | Kanazawa University | 骨密度計測装置 |
| JP2007007267A (ja) * | 2005-07-01 | 2007-01-18 | Kanazawa Univ | 骨密度計測装置 |
| JP2008155011A (ja) * | 2006-11-27 | 2008-07-10 | Kanazawa Univ | 密度計測装置およびその方法 |
| JP2010099414A (ja) * | 2008-10-27 | 2010-05-06 | Panasonic Electric Works Co Ltd | 骨密度計測装置 |
| JP2010175448A (ja) * | 2009-01-30 | 2010-08-12 | Kowa Co | 光学撮像装置 |
| JP2013527469A (ja) * | 2010-06-03 | 2013-06-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | ブルースター角を使用してビリルビンのような組織検体を測定する装置及び方法 |
| JP2015188590A (ja) * | 2014-03-28 | 2015-11-02 | カシオ計算機株式会社 | ダーモスコープ及びダーモスコープ用レンズユニット |
-
2017
- 2017-02-16 JP JP2017026832A patent/JP6951730B2/ja active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1043163A (ja) * | 1996-03-26 | 1998-02-17 | Hologic Inc | 骨特性評価システム |
| US20060192965A1 (en) * | 2005-01-21 | 2006-08-31 | Tromberg Bruce J | Method for assessing the condition of bone in-vivo |
| WO2006115207A1 (ja) * | 2005-04-22 | 2006-11-02 | Kanazawa University | 骨密度計測装置 |
| JP2007007267A (ja) * | 2005-07-01 | 2007-01-18 | Kanazawa Univ | 骨密度計測装置 |
| JP2008155011A (ja) * | 2006-11-27 | 2008-07-10 | Kanazawa Univ | 密度計測装置およびその方法 |
| JP2010099414A (ja) * | 2008-10-27 | 2010-05-06 | Panasonic Electric Works Co Ltd | 骨密度計測装置 |
| JP2010175448A (ja) * | 2009-01-30 | 2010-08-12 | Kowa Co | 光学撮像装置 |
| JP2013527469A (ja) * | 2010-06-03 | 2013-06-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | ブルースター角を使用してビリルビンのような組織検体を測定する装置及び方法 |
| JP2015188590A (ja) * | 2014-03-28 | 2015-11-02 | カシオ計算機株式会社 | ダーモスコープ及びダーモスコープ用レンズユニット |
Non-Patent Citations (1)
| Title |
|---|
| 三浦要ほか: "光を用いた骨密度計測装置の開発とモンテカルロ法による検証", 日本機械学会 北陸信越支部 第５３期総会・講演会 講演論文集, vol. 904, JPN6020041403, 1 March 2016 (2016-03-01), ISSN: 0004504814 * |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6951730B2 (ja) | 2021-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102531819B1 (ko) | 혈관내 광 간섭 단층촬영법을 사용한 자동 관상동맥 플라크 특징화 및 위험 평가를 위한 시스템들 및 방법들 | |
| Popescu et al. | Assessment of early demineralization in teeth using the signal attenuation in optical coherence tomography images | |
| JP5463545B2 (ja) | 濃度定量装置、濃度定量方法及びプログラム | |
| JP5201920B2 (ja) | 測定装置 | |
| EP2077753B1 (en) | Apparatus, probe and method for providing depth assessment in an anatomical structure | |
| JP5684443B2 (ja) | 生体成分測定装置 | |
| JP2009538156A (ja) | 組織のラマン分析 | |
| JP2011092780A (ja) | 生体情報イメージング装置、生体情報の解析方法、及び生体情報のイメージング方法 | |
| JPWO2007088709A1 (ja) | 光音響断層撮影装置及び光音響断層撮影方法 | |
| JP2010088873A (ja) | 生体情報イメージング装置 | |
| CA2895519C (en) | Apparatus for in-vitro imaging and analysis of dental samples | |
| JP2009018153A (ja) | 生体情報イメージング装置 | |
| KR20170041138A (ko) | 피검체 정보 취득장치 및 그 제어 방법 | |
| KR20080108918A (ko) | 초음파와 광의 상호 작용을 이용하여 피검체의 내부 정보를얻는 의료 기구 | |
| JP6461288B2 (ja) | 生体情報イメージング装置、生体情報の解析方法、及び生体情報のイメージング方法 | |
| JPWO2010131697A1 (ja) | 血管内壁分析装置及び血管内壁分析方法 | |
| JP6300615B2 (ja) | 光音響画像撮影装置及び光音響画像撮影方法 | |
| EP3797687A1 (en) | Blood vessel detection device and method therefor | |
| JP6486085B2 (ja) | 光音響波測定装置 | |
| JP6951730B2 (ja) | 骨密度計測装置及び方法 | |
| JP6431928B2 (ja) | 光音響顕微鏡及び光音響信号検出方法 | |
| JP2008155011A (ja) | 密度計測装置およびその方法 | |
| JP2010223770A (ja) | 散乱体内部観察装置および散乱体内部観察方法 | |
| JP2017006541A (ja) | 被検体情報取得装置および被検体情報取得方法 | |
| JP2008039635A (ja) | 弾性粘性測定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20191220 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20201023 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20201109 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20210105 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210302 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210519 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210708 |
|
| 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: 20210913 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210916 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6951730 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |