[go: up one dir, main page]

WO2018151150A1 - 散乱体濃度計測装置及びその方法 - Google Patents

散乱体濃度計測装置及びその方法 Download PDF

Info

Publication number
WO2018151150A1
WO2018151150A1 PCT/JP2018/005072 JP2018005072W WO2018151150A1 WO 2018151150 A1 WO2018151150 A1 WO 2018151150A1 JP 2018005072 W JP2018005072 W JP 2018005072W WO 2018151150 A1 WO2018151150 A1 WO 2018151150A1
Authority
WO
WIPO (PCT)
Prior art keywords
light
living body
irradiation
light intensity
wavelength
Prior art date
Application number
PCT/JP2018/005072
Other languages
English (en)
French (fr)
Inventor
一也 飯永
Original Assignee
メディカルフォトニクス株式会社
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 メディカルフォトニクス株式会社 filed Critical メディカルフォトニクス株式会社
Priority to JP2018568564A priority Critical patent/JP6894088B2/ja
Publication of WO2018151150A1 publication Critical patent/WO2018151150A1/ja

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration or pH-value ; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid or cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration or pH-value ; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid or cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated

Definitions

  • a step of detecting the received light intensity of 2 a step of calculating an absorption coefficient from the second received light intensity, a step of calculating a light scattering coefficient in the living body based on the first received light intensity and the absorption coefficient, and And calculating a scatterer concentration in the living body based on the light scattering coefficient.
  • the scatterer concentration measurement apparatus 1 includes a first irradiation unit 21 and a second irradiation unit that irradiate light from a predetermined irradiation position A of the living body E toward the inside of the living body E. 22, a first light intensity detection unit 31 and a second light intensity detection unit that detect the light intensity at a predetermined detection position B of the living body E, and a third light that detects the light intensity at the predetermined detection position C.
  • the light intensity detected by the intensity detector 33 and the first light intensity detector 31 is corrected using the light intensity detected by the second light intensity detector 32 and the third light intensity detector 33.
  • a control unit 4 for calculating the corrected to calculate the scattering coefficient mu s of the light in the light intensity based vivo, scatterers in vivo based on the scattering coefficient of the calculated optical mu s (lipid) concentration.
  • the second light intensity detection unit 32, and the third light intensity detection unit 33 are formed by increasing the distance between the first irradiation unit 21 and the second irradiation unit 22 and the first light intensity detection unit 31, the second light intensity detection unit 32, and the third light intensity detection unit 33. Since the optical path length becomes long, the number of collisions with lipids increases, and the detected light is greatly affected by scattering, so that it is easy to capture the influence of scattering that has been weak and difficult to detect.
  • the control unit 4 has an absorption coefficient calculated from the received light intensity detected by the first light intensity detector 31 and the received light intensity detected by the second light intensity detector 32 and the third light intensity detector 33. Based on the above, the light scattering coefficient ⁇ s in the living body is calculated. Note that the scattering coefficient ⁇ s in the embodiment is not limited to a numerical value obtained by quantifying the efficiency of a general scattering process, and is a value obtained by quantifying the influence of scattering under a certain condition in consideration of a scattering phenomenon. Is also included.
  • D is represents the diffusion coefficient
  • mu a is the absorption coefficient
  • the above equation (1) is expected to be applied to biological scatter measurement as a calculation formula that linearly approximates the measured light attenuation and derives the absorption coefficient from the intercept and the scattering coefficient from the slope.
  • Figure 6 shows the results of correction of absorption coefficient in phantom measurement. As a result of multiplying the correction coefficient due to water absorption at 950 nm, it was confirmed that the absorption coefficient could be calculated accurately.
  • the controller 4 causes the controller 4 to detect the light intensity of the first irradiation light detected by the first light intensity detector 31, the second light intensity detector 32, and Based on the absorption coefficient calculated from the received light intensity detected by the third light intensity detector 33, the light scattering coefficient ⁇ s in the living body is calculated.
  • the method for calculating the scattering coefficient ⁇ s has been described above.
  • the scatterer concentration can be calculated by acquiring the light intensity distribution emitted from the living body.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Optics & Photonics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

【課題】 吸収散乱複合体において、散乱係数を精度よく計測する。 【解決手段】 生体外から生体内に向けて、第1の照射強度で第1の波長の光を照射する第1の照射部と、生体外から生体内に向けて第2の照射強度で第2の波長の光を照射する第2の照射部と、第1の照射部による光の照射位置から第1の距離をあけた第1の検出位置における、生体から放出された第1の波長の光の第1の受光強度を検出する第1の光強度検出部と、第2の照射部による光の照射位置から第2の距離をあけた第2の検出位置における、生体から放出された第2の波長の光の第2の受光強度を検出する第2の光強度検出部と、第2の受光強度に基づき吸収係数を算出し、第1の受光強度と吸収係数に基づき生体内における光の散乱係数を算出し、算出された光の散乱係数に基づき生体内における散乱体濃度を算出する制御部とを有する。

Description

散乱体濃度計測装置及びその方法
  本発明は、散乱体濃度計測器及びその方法に関するものである。
 食後高脂血症の診断は、食後6~8時間の血中の脂質濃度変化を観測する必要がある。つまり、食後の高脂血状態を計測するためには、被験者を6~8時間拘束し、複数回の採血が必要である。そのため、食後高脂血症の診断は臨床研究の域を出ず、食後高脂血症の診断を臨床現場で実施することは、現実的ではなかった。
 このような課題を解決する手法が、特許文献1に開示されている。特許文献1の手法によれば、非侵襲脂質計測により、採血を無くすことができる。これにより医療機関のみならず家庭でも血中脂質を計測できるようになる。即時的なデータ取得を可能とすることで、時間的に連続した血中脂質を計測することが可能となる。
国際公開第2014/087825号公報
  特許文献1に記載された発明では、吸収係数および散乱係数の2変数が未知であっても、それぞれの値を正確に計測する。しかしながら、特許文献1の記載の計算式は、吸収係数と散乱係数が等しいか、散乱係数が吸収係数より大きいときに最大の効果を得ることができる。そのため、吸収散乱複合体のうち吸収が大きくなるような状態や箇所では、散乱係数を計測することが難しい場合もある。
  本発明は、このような問題点を解決するためになされたものであって、吸収散乱複合体において、散乱係数を精度よく計測することを目的とする。
 本発明の散乱体濃度計測装置は、生体外から生体内に向けて、第1の照射強度で第1の波長の光を照射する第1の照射部と、生体外から生体内に向けて第2の照射強度で第2の波長の光を照射する第2の照射部と、第1の照射部による光の照射位置から第1の距離をあけた第1の検出位置における、生体から放出された第1の波長の光の第1の受光強度を検出する第1の光強度検出部と、第2の照射部による光の照射位置から第2の距離をあけた第2の検出位置における、生体から放出された第2の波長の光の第2の受光強度を検出する第2の光強度検出部と、第2の受光強度に基づき吸収係数を算出し、第1の受光強度と吸収係数に基づき生体内における光の散乱係数を算出し、当該算出された光の散乱係数に基づき生体内における散乱体濃度を算出する制御部と、を有する。
 また、本発明の散乱体濃度計測方法は、生体外から生体内に向けて第1の照射強度で、第1の波長の光を照射する工程と、生体外から生体内に向けて第2の照射強度で、第2の波長の光を照射する工程と、第1の波長の光の照射位置から第1の距離をあけた第1の検出位置における、生体から放出された第1の波長の光の第1の受光強度を検出する工程と、第2の波長の光の照射位置から第2の距離をあけた第2の検出位置における、生体から放出された第2の波長の光の第2の受光強度を検出する工程と、第2の受光強度から吸収係数を算出する工程と、第1の受光強度と吸収係数に基づき生体内における光の散乱係数を算出する工程と、算出された光の散乱係数に基づき生体内における散乱体濃度を算出する工程とを有する。
 また、本発明の散乱体濃度計測装置は、生体外から生体内に向けて第1の照射強度で、第1の波長の光を照射する第1の照射部と、生体外から生体内に向けて第2の照射強度で、第2の波長の光を照射する第2の照射部と、第1の照射部による光の照射位置から第1の距離をあけた第1の検出位置における、生体から放出された第1の波長の光の第1の受光強度を検出する第1の光強度検出部と、第2の照射部による光の照射位置から第2の距離をあけた第2の検出位置における、生体から放出された第2の波長の光の第2の受光強度を検出する第2の光強度検出部と、を有するユーザー装置に、通信可能に接続される散乱体濃度計測装置であって、第2の受光強度から吸収係数を算出し、第1の受光強度と吸収係数に基づき生体内における光の散乱係数を算出し、算出された光の散乱係数に基づき生体内における散乱体濃度を算出する制御部を有する。
  本発明によれば、吸収散乱複合体において、散乱係数を精度よく計測することができる。
実施形態の散乱体濃度計測装置を示す図である。 血中脂質による光の散乱を示す図である。 実施形態の散乱体濃度計測装置のブロック図である。 モンテカルロシミュレーションの結果を示す図である。 式(1)で求めた吸収係数と理論上の吸収係数の相関を示す図である。 式(1)で求めた吸収係数と理論上の吸収係数の相関を示す図である。 実施形態の散乱体濃度計測のフローチャートである。 実施形態の散乱体濃度計測システムを示す図である。 実施形態の散乱体濃度計測装置のブロック図である。 式(1)をプロットした図である。
  以下、実施形態の散乱体濃度計測装置及びその方法について、図を参照して詳細に説明をする。また、本実施形態では、散乱体の例として、血中脂質を検出する場合について主に説明するが、これに限られず、散乱体一般についても適用可能である。
 図1は、実施形態の散乱体濃度計測装置の構成を示す図である。
  図1に示すように、実施形態の散乱体濃度計測装置1は、生体Eの所定の照射位置Aから生体Eの内部に向けて光を照射する第1の照射部21及び第2の照射部22と、生体Eの所定の検出位置Bにおける光強度を検出する第1の光強度検出部31及び第2の光強度検出部と、所定の検出位置Cにおける光強度を検出する第3の光強度検出部33と、第1の光強度検出部31により検出された光強度を、第2の光強度検出部32及び第3の光強度検出部33により検出された光強度を用いて補正し、補正された光強度基づき生体内における光の散乱係数μを算出し、算出された光の散乱係数μに基づき生体内における散乱体(脂質)濃度を算出する制御部4とを有する。
  第1の照射部21及び第2の照射部22は、生体外から生体に向けて、所定の照射位置Aに光を照射する。第1の照射部21及び第2の照射部22は、相異なる波長の光を照射する。第1の照射部21及び第2の照射部22は、光を照射するための光源を有し、照射する光の波長を自在に調整することができる。実施形態では、第1の照射部21は散乱係数の計測を目的とし、第2の照射部22は吸収係数の計測を目的とする。
 実施形態では、第1の照射部21及び第2の照射部22と、2つの照射部としているが、例えば、照射部に白色光を用い、光強度検出部側で可視光フィルターなどにより分光をすれば照射部を1つとするができる。また、照射部が白色LEDのような複数波長の発光体の複合体の場合には、その中で複数波長を交互点灯させるなどにより、照射部を1つとすることができる。
  第1の照射部21及び第2の照射部22は、後述する制御部4による散乱係数μの算出方法に応じて、光の連続的な照射や光のパルス状の照射等の光を照射する時間長さを任意に調整することができ、かつ照射する光の強度または光の位相を任意に変調することができる。
  第1の照射部21及び第2の照射部22は、波長が固定された光源を用いてもよく、複数の波長の光を混合したものであってもよい。第1の照射部21及び第2の照射部22は白色光などの広範囲のスぺクトルを有する光源でもよく、その場合、受光部で特定波長を検出するようにしても良い。
  第1の光強度検出部31、第2の光強度検出部32、及び第3の光強度検出部33は、光を受光してその受光強度を検出する。第1の光強度検出部31は、第1の照射部21が照射した光が血中脂質により散乱し、生体から生体外に放出される光を、検出位置Bで受光し、その受光強度を検出する。第2の光強度検出部32は、第2の照射部22が照射した光が血中脂質により散乱し、生体から生体外に放出される光を、検出位置Bで受光し、その受光強度を検出する。第3の光強度検出部33は、第2の照射部22が照射した光が血中脂質により散乱し、生体から生体外に放出される光を、検出位置Cで受光し、その受光強度を検出する。
 実施形態では、図1に示すように、第1の照射部21から所定の間隔ρで同一面上でかつ直線状に第1の光強度検出部31が並べられている。また、第2の照射部22から所定の間隔ρ、ρで同一面上でかつ直線状に第2の光強度検出部32及び第3の光強度検出部33が順に並べられている。
  図1に示すように、第1の照射部21の照射位置Aから第1の光強度検出部311の検出位置Bまでの距離を第1の照射検出間距離ρとし、第2の照射部22の照射位置Aから第2の光強度検出部31までの距離を第2の照射検出間距離ρとし、第2の照射部22の照射位置Aから第3の光強度検出部33までの距離を第3の照射検出間距離ρとする。
  このように、光を生体に照射する第1の照射部21及び第2の照射部22と、生体から放出される受光強度を検出する第1の光強度検出部31、第2の光強度検出部32及び第3の光強度検出部33との間に所定の距離を設けることにより、図2に示すように、照射した光が生体表面および表面近傍の散乱体Dにより反射して直接的に生体Eから放出される光の影響を抑制し、血液や脂質が存在する深さに達したのち、血中脂質によって光が反射することによる散乱を経て生体から放出される後方散乱光による受光強度を検出する。また、第1の照射部21及び第2の照射部22と第1の光強度検出部31、第2の光強度検出部32及び第3の光強度検出部33との距離を長くすることで、光路長は長くなるため、脂質との衝突回数が増え、検出される光は散乱の影響を多く受けることにより、これまでは弱く、検出しにくかった散乱の影響を捉えやすくしている。
  複数の光強度検出部を設ける場合の配列は、第1の照射部21及び第2の照射部22を中心として各々異なる距離に配置されるのであれば直線状に限定されるものではなく、円状、波状、ジグザグ状などでもよい。また、第1及び第2の照射検出間距離ρや第3の照射検出間距離ρは、一定の間隔でなくともよい。また、実施形態では、また、第1の照射検出間距離と第2の照射検出間距離をρとしたが、異なっていてもよい。第1の光強度検出部31、第2の光強度検出部32及び第3の光強度検出部33は、フォトダイオードやCCDやCMOS等の受光素子であってもよい。
 次に、散乱体濃度計測装置1の制御系の構成について説明する。図3は実施形態の散乱体濃度計測装置1のブロック図である。システムバス109を介して、CPU(Central Processing Unit)104、ROM(Read Only Memory)105、RAM(Random Access Memory)106、記憶部107、外部I/F(Interface)108、第1の照射部21、第2の照射部22、第1の光強度検出部31、第2の光強度検出部32、及び、第3の光強度検出部33が接続される。CPU104とROM105とRAM106とで制御部(コントローラー)4を構成する。
 ROM105は、CPU104により実行されるプログラムや閾値を予め記憶する。
 RAM106は、CPU104が実行するプログラムを展開するエリアと、プログラムによるデータ処理の作業領域となるワークエリアなどの様々なメモリエリア等を有する。
 記憶部107は、予め用意された散乱係数μsと散乱体血中脂質濃度との統計データを記憶する。記憶部107は、HDD(Hard Disk Drive)や、フラッシュメモリや、SSD(Solid State Drive)等の、不揮発性に記憶する内部メモリーでよい。
 外部I/F108は、例えばクライアント端末(PC)などの外部装置と通信するためのインターフェースである。外部I/F108は、外部装置とデータ通信を行うインターフェースであれば良く、たとえば、外部装置にローカルに接続する機器(USBメモリ等)であっても良いし、ネットワークを介して通信するためのネットワークインターフェイスであっても良い。
 制御部4は、第2の光強度検出部32及び第3の光強度検出部33により検出された受光強度から、傾きを求め吸収係数を算出する。吸収係数の算出法は、以下の通りである。ここでは、例えば、下記式1を用いる場合について示す。
 下記式1は、y = -ax + bと同義である。そこで、yはR(ρ)、ρ、及び、Sが分かれば算出でき、得られたyをρに対しプロットすれば、図10に示すプロットとなる。なお、R(ρ)は受光強度(光強度)であり、ρは入射?受光部間距離であり、Sは、入射光強度である。プロットの傾きが有効減衰係数μeffなので、有効減衰係数μeffを右の項に代入すると吸収係数μaが求まる。実施形態では、第2の光強度検出部32及び第3の光強度検出部33で得られた光強度を図10に示すようにプロットして一次関数に近似させ、当該一次関数の傾きからμeffを求め、これを式(1)に代入して吸収係数μaを求める。
 なお、実施形態では、第2の光強度検出部32及び第3の光強度検出部33により検出された受光強度から、吸収係数を算出しているが、例えば、吸収係数を事前に測定することにより図10に示す検量線等の吸収係数の統計データを取得し、当該吸収係数の統計データをROM105や記憶部107に等に格納しておき、当該検量線データと第2の光強度検出部32の受光強度から、吸収係数を算出することもできる。この場合には、光強度検出部を第2の光強度検出部32の1つとすることもできる。
 制御部4は、第1の光強度検出部31により検出された受光強度と、第2の光強度検出部32及び第3の光強度検出部33により検出された受光強度から算出された吸収係数とに基づき生体内における光の散乱係数μを算出する。なお、実施形態における散乱係数μは、一般的な散乱過程の効率を数値化したものに限定されるものではなく、散乱現象を考慮して散乱の影響を一定の条件下で数値化したものも含むものである。
  実施形態の制御部4は、第1の照射部21により連続光を照射するとともに、第1の光強度検出手段31により検出された受光強度R(ρ)と照射検出間距離ρとを、下記式(1)「数1」および式(2)「数2」に代入することで散乱係数μを算出する。
Figure JPOXMLDOC01-appb-M000003
Figure JPOXMLDOC01-appb-M000004
  ここで、μは吸収係数、μeffは有効減衰係数(Effective Attenuation Coefficient)、Sは照射部により照射された光の光強度である。
  なお、上記式(1)および式(2)は以下のように導き出される。
  まず、図1に示すように、生体外から生体内に向けて所定の光強度Sを有する光を連続光として照射するとともに、第1の照射部21の照射位置Aから第1の光強度検出部31の検出位置Bまでの距離を照射検出間距離ρとすると、その後方散乱光によって生体外に放射される光の分布は、以下の式(3)で表される。
Figure JPOXMLDOC01-appb-M000005
  ここで、zは光源の深さ、つまり散乱が開始する深さであって、下記式(4)で表される。
Figure JPOXMLDOC01-appb-M000006
  ここで、μは散乱係数を表している。
  また、μeffは有効減衰係数であり、下記式(5)で表される。
Figure JPOXMLDOC01-appb-M000007
  ここで、Dは拡散係数、μは吸収係数をそれぞれ表している。
  また、皮膚表面や表面近傍の血管での散乱を想定すると、照射検出間距離ρと光源の深さzとの関係は下記式(6)のように近似する事ができる。
Figure JPOXMLDOC01-appb-M000008
  さらに、本実施形態における計測対象は、上述のとおり血中の脂質であり、血中脂質による散乱は吸収よりも大きいと考えられる。そのため、有効減衰係数μeffは下記式(7)のように近似する事ができる。
Figure JPOXMLDOC01-appb-M000009
  以上の式(6)および式(7)を式(3)に代入すると、下記式(8)の近似式となる。
Figure JPOXMLDOC01-appb-M000010
  ここで、照射検出間距離ρと有効減衰係数μeffとに関して、下記式(9)のような関係を有する場合、式(8)は下記式(10)のように表される。ここでμeff=5.77mm(μ=1/mm、μ=0.01/mm)とする。
Figure JPOXMLDOC01-appb-M000011
Figure JPOXMLDOC01-appb-M000012
  そして、上記式(10)を対数表示させると、上記式(1)が導き出される。
  また、照射検出間距離ρと有効減衰係数μeffとに関して、下記式(11)のような関係を有する場合、式(8)は下記式(12)のように表される。
Figure JPOXMLDOC01-appb-M000013
Figure JPOXMLDOC01-appb-M000014
  そして、上記式(12)を対数表示させると、上記式(2)が導き出される。
  なお、制御部4における散乱係数算出は、本実施形態のように上記式(1)および式(2)によるものに限定されるものではなく、適宜選択されるものであり、例えば、検出された受光強度R(ρ)と散乱係数μとが単純に比例しているものとしてもよい。
  また、制御部4における散乱係数算出は、第1の光強度検出部31の検出位置Aが一点のものに限定されるものではない。実際の計測に置いては、多くの計測ノイズが発生することが想定される。そのような場合は、検出位置を多数設置し、照射検出間距離ρに応じた連続的な受光強度から散乱係数を導くこともできる。つまり、制御部4における散乱係数算出において、計測点が少数である各計測データのノイズが相対的に大きくなる場合、検出位置を増やすことで、実測で想定されるノイズの影響を軽減させることが可能である。
 上記式(1)においては、計測した光の減衰度を線形近似し、切片から吸収係数、傾きから散乱係数を導き出す計算式として、生体の散乱計測などに応用が期待されている。
 しかしながら、上記式(1)の導出においては、拡散方程式の展開の過程で、散乱係数μs’>>吸収係数μaという仮定を設定しているが、具体的な計算条件とはいいがたい。(なお、「>>」は、非常に大であることを示す記号である。)
 例えば、モンテカルロシミュレーションの結果では、図4に示した通り、吸収係数の増加に伴い式(1)で導かれた散乱係数は、増加するという傾向が得られた。また、式(1)で散乱係数を算出する工程で、理論値の吸収係数を代入した場合は、目標値である散乱係数1.5を導き出すことが分かった。
 すなわち、上記式(1)においては、切片から吸収係数を算出する工程に改良の余地があると考えられる。
 図5では、上記式(1)で求めた吸収係数と理論上の吸収係数の相関であるが、相関係数はr2=0.99以上であり、申し分のない相関であるが、傾きを持っていることがわかる。
 一見すると、この傾きを代入することで補正が可能に見えるが、異なる散乱係数で同様の分析を行うと、傾きに違いがみられる。
 すなわち、散乱係数が既知でなければ、一定水準以上の吸収を持つ吸収散乱複合体の散乱係数を導出することは困難である。
 これらの課題は、上記式(1)の導出における、散乱係数μs’>>吸収係数μaという条件において、この関係の具体性が不明瞭であるため生じたものである。
 ここで吸収係数μa>>散乱係数μs’という仮定を置くと、式(1)で計測される傾きの有効減衰係数μeffは、吸収係数μaに相当することとなる。すなわち以下の式である。
Figure JPOXMLDOC01-appb-M000015
 つまり未知物質であっても、波長により散乱係数μs’>>吸収係数μaと、吸収係数μa>>散乱係数μs’が混在する吸収散乱複合体は、それぞれに特徴的な波長域を用いることができれば、1回の計測で、計算により2つの未知数すなわち、散乱係数と吸収係数の絶対値を精度良く導き出すことができる。
 なお、生体における吸収散乱同時分析を行うためには、散乱係数μs’>>吸収係数μaの領域の照射光の波長は750nm以上1600nm以下が望ましく、吸収係数μa>>散乱係数μs’の領域の照射光の波長は、350nm以上750nm以下, 又は、1600nm以上2000nm以下が望ましい。
 波長350nm以上750nm以下は、ヘモグロビンやタンパクなどの成分依存が多く、波長1600nm以上2000nm以下は水の吸収を対象としている。水の吸収を活用する場合は、900nm以上の波長でもよい。これは、水の吸収が広範囲に分布しているため、測定対象により計測精度向上のための散乱係数と吸収係数の強度バランスが異なるためである。
 精度向上のために、吸収係数は短波長領域と長波長領域の2種を用いてもよい。
 方法としては、まず吸収係数μa>>散乱係数μs’の領域の照射光の波長範囲にて吸収係数μaを求め、散乱係数μs’>>吸収係数μaの領域の照射光の波長範囲で、式(1)に実測の吸収係数μaを代入する。
 しかしながら、吸収係数μa>>散乱係数μs’の領域と、散乱係数μs’>>吸収係数μaの領域では波長が異なるため、実測の吸収係数μaをそのまま用いることはできない。
 そこで、例えば、吸収係数μa800と吸収係数μa650の換算係数を用いることで、吸収係数μa>>散乱係数μs’の領域の吸収係数μaから、散乱係数μs’>>吸収係数μaの領域における吸収係数μaを求め、式(1)から散乱係数μs’を求めることが可能となる。また、有効減衰係数μeff650=2μa650となる。有効減衰係数μeff650の値から、散乱係数を計測するための吸収係数μaを以下のように計算式を選択することで、散乱係数の絶対値を算出することが可能となる。
 μa650>0.5のとき (計算式a) μa800 = A・μa650
 0.5≧μa650≧0.1のとき (計算式b) μa800 =B・μa650
 μa650<0.1のとき (計算式c) μa800 =C・μa650
 A,B,Cは、それぞれの補正係数である。また、それぞれの係数は、計算式であってもよく、対数やべき乗を含んでもよく、三角関数やπ、粒子の体積補正などを含んでもよい。換算係数は、ROM105や記憶部107に格納されればよい。
 図6に、ファントム計測における吸収係数の補正の結果を示す。950nmにおける水の吸収による補正係数を乗じた結果、吸収係数を正確に算出できることが確認できた。
  制御部4は、散乱係数μに基づいて血中脂質の濃度を算出する。散乱係数μと脂質濃度とは相関があり、散乱係数μの値に基づいて脂質濃度を算出することができる。実施形態では、散乱係数μと散乱体濃度との関係についてあらかじめ統計データを取り、ROM105や記憶部107に等に格納して、算出された散乱係数μと、格納された統計データとを比較することにより、実際の散乱体濃度を算出する。
  例えば、特定の生体A氏の血中脂質濃度を計測対象とする場合は、A氏の血中脂質濃度を採血などの他の血中脂質濃度計測方法等により事前に計測した計測結果をROM105や記憶部107に格納して、この計測結果と算出された散乱係数μとを比較して、A氏個人の統計データを作成して、濃度を算出できるようにすることができる。
  若しくは、A氏の血中脂質の濃度を他の血中脂質の濃度の測定方法等により事前に測定した測定結果をROM105や記憶部107等に格納して、この測定結果と検出された光強度より得られた濃度の測定結果とを比較して、その比較により得られた濃度と、一般的な生体の場合の上記統計データにおける濃度との誤差を算出し、その誤差を修正するキャリブレーションすることで、A氏個人の統計データを作成してもよい。
  なお、統計データの形式は特に限定されるものではなく、例えば、性別、身長、体重、BMI等で分類されていてもよく、表やグラフ、関数式等を用いて算出できるようにしてもよい。
  また、臨床現場において、濃度と濁度とは同義で使われることがあり、本発明における濃度には濁度の概念も含まれる。よって、脂質濃度算出手段は、その算出結果として、濃度のみならず単位量当たりの粒子数やホルマジン濁度とすることができる。
 次に、本実施形態の散乱体濃度計測方法について説明する。図7は、実施形態の散乱体濃度計測のフローチャートである。
 照射工程(S101)では、第1の照射部21及び第2の照射部22が、生体の照射位置に対して第1の波長の第1の照射光及び第2の波長の第2の照射光を照射する。
 第1の波長は、750nm以上1600nm以下であってよく、第2の波長は、350nm以上750nm以下, 又は、1600nm以上2000nm以下であってよい。
 光強度検出工程(S102)では、第1の光強度検出部31、第2の光強度検出部32及び第3の光強度検出部33が、生体から放出される相異なる波長の第1の照射光及び第2の照射光を各々受光して、各々光強度を検出する。光強度検出工程で検出された光強度は、吸収係数算出工程へと送られる。
 吸収係数算出工程(S103)では、制御部4は、第2の光強度検出部32及び第3の光強度検出部33により検出された受光強度から、傾きを求め、吸収係数を算出する。なお、光強度検出部を第2の光強度検出部32及び第3の光強度検出部33の2つとした場合の吸収係数の算出法、及び、光強度検出部を第2の光強度検出部32の1つとした場合の吸収係数の算出法については上述した。
 散乱係数算出工程(S104)では、制御部4は、制御部4は、第1の光強度検出部31により検出された第1の照射光の光強度と、第2の光強度検出部32及び第3の光強度検出部33により検出された受光強度から算出された吸収係数と、に基づき生体内における光の散乱係数μsを算出する。散乱係数μの算出法については上述した。
 散乱体濃度算出工程(S105)では、制御部4は、散乱係数μに基づいて、血中の散乱体濃度を算出する。なお、血中の散乱体濃度の算出法については、上述した。散乱体濃度算出工程では、散乱係数を算出した後、脂質濃度(散乱体濃度)を算出してもよい。
 以上説明したように、本実施形態の散乱体濃度計測装置及び方法によれば、生体から放出される光強度分布を取得することで、散乱体濃度の算出が可能となる。
 次に、本発明の他の実施形態の散乱体濃度計測装置について説明をする。なお、本発明の他の実施形態の散乱体濃度計測装置の構成は、上記実施形態の散乱体濃度計測装置の構成と共通する部分もあるため、相違する部分を主に説明する。
 上記実施形態では、第1の照射部21及び第2の照射部22と、第1の光強度検出部31、第2の光強度検出部32及び第3の光強度検出部33と、制御部4とを一体として構成した例を示したが、これに限られず、第1の照射部21及び第2の照射部22と、第1の光強度検出部31、第2の光強度検出部32及び第3の光強度検出部33をユーザー装置として構成し、制御部4を、ユーザー装置に接続したサーバー装置に設けたシステムとしてもよい。
 図8は、実施形態の散乱体濃度計測システムの構成を示す図である。システムは、散乱体濃度計測装置200と、アクセスポイント300と、ユーザー装置400とを有する。
 散乱体濃度計測装置200は、ユーザー装置400から送信された光強度に基づいて所定の処理を行い、散乱体濃度を算出するための装置であり、具体的には、パーソナルコンピュータや、装置の台数や送受信するデータ量によってはサーバー装置が適宜用いられる。
 ユーザー装置400は、ユーザーが所持する装置であり、単独の装置である場合もあり、スマートフォン、携帯電話、腕時計等に搭載される場合もある。また、第1の照射部421と、第2の照射部422と、第1の光強度検出部431と、第2の光強度検出部432と、第3の光強度検出部433と、通信部404として、スマートフォンや携帯電話に備わるカメラや照明、通信機能等を使用してもよい。
 ユーザー装置400は、光を照射する第1の照射部421と、第2の照射部422と、第1の光強度検出部431と、第2の光強度検出部432と、第3の光強度検出部433と、通信部404とを有する。通信部404は、第1の光強度検出部431、第2の光強度検出部432、及び、第3の光検出部433とで各々検出された光強度を送信する。第1の照射部421と、第2の照射部422と、第1の光強度検出部431と、第2の光強度検出部432と、第3の光強度検出部433の構成や機能・動作は、上述した第1の照射部21と、第2の照射部22と、第1の光強度検出部31と、第2の光強度検出部32と、第3の光強度検出部33と各々同様である。
 散乱体濃度計測装置200は、通信部204と制御部203とを有する。通信部204は、通信部404から送信された光強度をアクセスポイント300を介して受信し、制御部203へ送信する。
 散乱体濃度計測装置200は、通信部204と制御部203とを有する。通信部204は、通信部404から送信された光強度をアクセスポイント300を介して受信し、制御部203へ送信する。制御部203の機能・動作は、上述した制御部103と同様である。
 次に、散乱体濃度計測装置200の制御系の構成について説明する。図9は、実施形態の散乱体濃度計測装置200のブロック図である。システムバス209を介して、CPU(Central Processing Unit)204、ROM(Read Only Memory)205、RAM(Random Access Memory)206、記憶部207、及び、通信部(外部I/F(Interface))208が接続される。CPU204とROM205とRAM206とで制御部(コントローラー)203を構成する。
 ROM205は、CPU204により実行されるプログラムや閾値を予め記憶する。
 RAM206は、CPU204が実行するプログラムを展開するエリアと、プログラムによるデータ処理の作業領域となるワークエリアなどの様々なメモリエリア等を有する。
 記憶部207は、予め用意された、静的パラメータ及び動的パラメータの適切な数値範囲のデータを記憶する。記憶部207は、HDD(Hard Disk Drive)や、フラッシュメモリや、SSD(Solid State Drive)等の、不揮発性に記憶する内部メモリーでよい。
 通信部(外部I/F)208は、例えばクライアント端末(PC)などの外部装置と通信するためのインターフェースである。外部I/F208は、外部装置とデータ通信を行うインターフェースであれば良く、たとえば、外部装置にローカルに接続する機器(USBメモリ等)であっても良いし、ネットワークを介して通信するためのネットワークインターフェイスであっても良い。
 なお、本実施形態では、ユーザー装置400から散乱体計測装置200へ、アクセスポイント300を介して光強度を送信したが、これに限られず、ユーザー装置400と散乱体計測装置200とが、ネットワークNを介さずに直接接続し、有線通信や無線通信等の手段により光強度を送信してもよい。

Claims (8)

  1.  生体外から生体内に向けて、第1の照射強度で第1の波長の光を照射する第1の照射部と、
     生体外から生体内に向けて、第2の照射強度で第2の波長の光を照射する第2の照射部と、
     前記第1の照射部による光の照射位置から第1の距離をあけた第1の検出位置における、前記生体から放出された前記第1の波長の光の第1の受光強度を検出する第1の光強度検出部と、
     前記第2の照射部による光の照射位置から第2の距離をあけた第2の検出位置における、前記生体から放出された前記第2の波長の光の第2の受光強度を検出する第2の光強度検出部と、
     前記第2の受光強度に基づき吸収係数を算出し、
     前記第1の受光強度と前記吸収係数に基づき生体内における光の散乱係数を算出し、
     当該算出された光の散乱係数に基づき生体内における散乱体濃度を算出する、制御部と、
      を有する散乱体濃度計測装置。
  2.  前記第2の照射部による光の照射位置から第3の距離をあけた第3の検出位置における、前記生体から放出された前記第2の波長の光の第3の受光強度を検出する第3の光強度検出部をさらに備え、
     前記制御部は、前記第2の受光強度及び第3の受光強度に基づき前記吸収係数を算出する、ことを特徴とする請求項1に記載の散乱体濃度計測装置。
  3.  前記第1の波長は、750nm以上1600nm以下であることを特徴とする請求項1または2に記載の散乱体濃度計測計測器。
  4.  前記第2の波長は、350nm以上750nm以下, 又は、1600nm以上2000nm以下であることを特徴とする請求項1から3のいずれかに記載の散乱体濃度計測計測器。
  5.   前記制御部は、前記第1の受光強度R(ρ)と前記第1の距離ρと前記吸収係数μaを、下記式(1)及び式(2)に代入して散乱係数μを算出する請求項1から4のいずれかに記載の散乱体濃度計測装置。
    Figure JPOXMLDOC01-appb-M000001
    Figure JPOXMLDOC01-appb-M000002
  6.  前記制御部は、前記第2の波長における吸収係数と前記第1の波長における吸収係数の換算係数により、前記第2の波長における吸収係数を、前記第1の波長における吸収係数に換算する、ことを特徴とする請求項1から5のいずれかに記載の散乱体濃度計測装置。
  7.   生体外から生体内に向けて第1の照射強度で、第1の波長の光を照射する工程と、
     生体外から生体内に向けて第2の照射強度で、第2の波長の光を照射する工程と、
      前記第1の波長の光の照射位置から第1の距離をあけた第1の検出位置における、前記生体から放出された前記第1の波長の光の第1の受光強度を検出する工程と、
     前記第2の波長の光の照射位置から第2の距離をあけた第2の検出位置における、前記生体から放出された前記第2の波長の光の第2の受光強度を検出する工程と、
     前記第2の受光強度に基づき吸収係数を算出する工程と、
     前記第1の受光強度と前記吸収係数に基づき生体内における光の散乱係数を算出する工程と、
     当該算出された光の散乱係数に基づき生体内における散乱体濃度を算出する工程と、
      を有する散乱体濃度計測方法。
  8.   生体外から生体内に向けて第1の照射強度で、第1の波長の光を照射する第1の照射部と、生体外から生体内に向けて第2の照射強度で、第2の波長の光を照射する第2の照射部と、前記第1の照射部による光の照射位置から第1の距離をあけた第1の検出位置における、前記生体から放出された前記第1の波長の光の第1の受光強度を検出する第1の光強度検出部と、前記第2の照射部による光の照射位置から第2の距離をあけた第2の検出位置における、前記生体から放出された前記第2の波長の光の第2の受光強度を検出する第2の光強度検出部と、を有するユーザー装置に、通信可能に接続される散乱体濃度計測装置であって、
     前記第2の受光強度に基づき吸収係数を算出し、前記第1の受光強度と前記吸収係数に基づき生体内における光の散乱係数を算出し、当該算出された光の散乱係数に基づき生体内における散乱体濃度を算出する制御部を有する散乱体濃度計測装置。
     
PCT/JP2018/005072 2017-02-14 2018-02-14 散乱体濃度計測装置及びその方法 WO2018151150A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018568564A JP6894088B2 (ja) 2017-02-14 2018-02-14 散乱体濃度計測装置及びその方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017025311 2017-02-14
JP2017-025311 2017-02-14

Publications (1)

Publication Number Publication Date
WO2018151150A1 true WO2018151150A1 (ja) 2018-08-23

Family

ID=63169349

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/005072 WO2018151150A1 (ja) 2017-02-14 2018-02-14 散乱体濃度計測装置及びその方法

Country Status (2)

Country Link
JP (1) JP6894088B2 (ja)
WO (1) WO2018151150A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020138225A1 (ja) * 2018-12-26 2020-07-02 メディカルフォトニクス株式会社 吸収カロリー計測装置、吸収カロリー計測方法、及び、吸収カロリー計測プログラム
WO2020246455A1 (ja) * 2019-06-04 2020-12-10 メディカルフォトニクス株式会社 無侵襲計測装置、方法およびプログラム
WO2022262221A1 (zh) * 2021-06-18 2022-12-22 中国科学院深圳先进技术研究院 一种红细胞比容的检测方法、装置、终端和可读存储介质
JP2024508717A (ja) * 2021-02-11 2024-02-28 先陽科技有限公司 生体組織成分の検出方法、装置及びウェアラブル機器

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004517300A (ja) * 2000-07-28 2004-06-10 アボット・ラボラトリーズ 組織の水分補給効果の補償による分析物濃度の非観血的決定方法
JP2013070822A (ja) * 2011-09-27 2013-04-22 Seiko Epson Corp 濃度定量装置、光吸収係数算出方法、等価散乱係数算出方法、濃度定量方法、光吸収係数の算出を行うプログラム及び濃度の算出を行うプログラム
WO2014087825A1 (ja) * 2012-12-06 2014-06-12 国立大学法人北海道大学 非侵襲型生体脂質濃度計測器、非侵襲型生体脂質代謝機能計測器、非侵襲による生体脂質濃度計測方法および非侵襲による生体脂質代謝機能検査方法
WO2016041073A1 (en) * 2014-09-17 2016-03-24 2352409 Ontario Inc. Device and method for monitoring fat balance

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5595534B2 (ja) * 2007-05-15 2014-09-24 キヤノン株式会社 生体情報イメージング装置、生体情報の解析方法、及び生体情報のイメージング方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004517300A (ja) * 2000-07-28 2004-06-10 アボット・ラボラトリーズ 組織の水分補給効果の補償による分析物濃度の非観血的決定方法
JP2013070822A (ja) * 2011-09-27 2013-04-22 Seiko Epson Corp 濃度定量装置、光吸収係数算出方法、等価散乱係数算出方法、濃度定量方法、光吸収係数の算出を行うプログラム及び濃度の算出を行うプログラム
WO2014087825A1 (ja) * 2012-12-06 2014-06-12 国立大学法人北海道大学 非侵襲型生体脂質濃度計測器、非侵襲型生体脂質代謝機能計測器、非侵襲による生体脂質濃度計測方法および非侵襲による生体脂質代謝機能検査方法
WO2016041073A1 (en) * 2014-09-17 2016-03-24 2352409 Ontario Inc. Device and method for monitoring fat balance

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020138225A1 (ja) * 2018-12-26 2020-07-02 メディカルフォトニクス株式会社 吸収カロリー計測装置、吸収カロリー計測方法、及び、吸収カロリー計測プログラム
JP2020106318A (ja) * 2018-12-26 2020-07-09 メディカルフォトニクス株式会社 吸収カロリー計測装置、吸収カロリー計測方法、及び、吸収カロリー計測プログラム
WO2020246455A1 (ja) * 2019-06-04 2020-12-10 メディカルフォトニクス株式会社 無侵襲計測装置、方法およびプログラム
JP2024508717A (ja) * 2021-02-11 2024-02-28 先陽科技有限公司 生体組織成分の検出方法、装置及びウェアラブル機器
WO2022262221A1 (zh) * 2021-06-18 2022-12-22 中国科学院深圳先进技术研究院 一种红细胞比容的检测方法、装置、终端和可读存储介质

Also Published As

Publication number Publication date
JPWO2018151150A1 (ja) 2020-01-09
JP6894088B2 (ja) 2021-06-30

Similar Documents

Publication Publication Date Title
JP5463545B2 (ja) 濃度定量装置、濃度定量方法及びプログラム
WO2018151150A1 (ja) 散乱体濃度計測装置及びその方法
JP6108705B2 (ja) 被検体情報取得装置、および、被検体情報取得方法
WO2012153769A1 (ja) 光トモグラフィ装置
US9259486B2 (en) Method and system for calculating a quantification indicator for quantifying a dermal reaction on the skin of a living being
CN103222880B (zh) 被检物信息获取装置
WO2019225612A1 (ja) 血管検知装置及びその方法
US20210401334A1 (en) Device for measuring blood lipid concentration and method therefor
US9860455B2 (en) Object information acquiring apparatus and signal processing method
EP2732756B1 (en) Object information acquisition apparatus
JP5924658B2 (ja) 濃度定量装置、光吸収係数算出方法、等価散乱係数算出方法、濃度定量方法、光吸収係数の算出を行うプログラム及び濃度の算出を行うプログラム
WO2015174273A1 (ja) 乳房計測方法及び計測装置
JP2008155011A (ja) 密度計測装置およびその方法
JP6894087B2 (ja) 脂質計測装置及びその方法
JP5881058B2 (ja) 小動物の体脂肪測定方法
US20210059575A1 (en) Lipid concentration measurement device and method therefor
RU2413930C1 (ru) Способ определения оптических характеристик однородного рассеивающего вещества
WO2020246455A1 (ja) 無侵襲計測装置、方法およびプログラム
JP2022149431A (ja) 総蛋白濃度計測装置及びその方法
JP2022100907A (ja) 脂質濃度計測装置、プログラム、及び、方法
JP2014006063A (ja) 光吸収係数算出装置、濃度定量装置、光吸収係数算出方法、濃度定量方法、光吸収係数の算出を行うプログラム及び濃度の算出を行うプログラム
TW200914814A (en) Method and apparatus for estimating the content of an analyte in a multi-layer medium

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18754929

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2018568564

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18754929

Country of ref document: EP

Kind code of ref document: A1