JP2022061939A - 分析物質の濃度推定装置及び方法、ならびに信号測定装置 - Google Patents

分析物質の濃度推定装置及び方法、ならびに信号測定装置 Download PDF

Info

Publication number
JP2022061939A
JP2022061939A JP2021084569A JP2021084569A JP2022061939A JP 2022061939 A JP2022061939 A JP 2022061939A JP 2021084569 A JP2021084569 A JP 2021084569A JP 2021084569 A JP2021084569 A JP 2021084569A JP 2022061939 A JP2022061939 A JP 2022061939A
Authority
JP
Japan
Prior art keywords
concentration
wavelength
light
path length
spectrum
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.)
Pending
Application number
JP2021084569A
Other languages
English (en)
Inventor
ホ リ,ジュン
Jun Ho Lee
モ アン,スン
Sung Mo Ahn
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.)
Samsung Electronics Co Ltd
Original Assignee
Samsung Electronics 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 Samsung Electronics Co Ltd filed Critical Samsung Electronics Co Ltd
Publication of JP2022061939A publication Critical patent/JP2022061939A/ja
Pending legal-status Critical Current

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, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14546Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring analytes not otherwise provided for, e.g. ions, cytochromes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0062Arrangements for scanning
    • A61B5/0066Optical coherence imaging
    • 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
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0075Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
    • 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, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14532Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
    • 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, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/0209Low-coherence interferometers
    • G01B9/02091Tomographic interferometers, e.g. based on optical coherence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/0205Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/45Interferometric spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06Investigating concentration of particle suspensions
    • G01N15/075Investigating concentration of particle suspensions by optical means
    • 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
    • G01N21/59Transmissivity
    • G01N21/5907Densitometers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/10Arrangements of light sources specially adapted for spectrometry or colorimetry
    • G01J2003/102Plural sources
    • 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
    • G01N2021/178Methods for obtaining spatial resolution of the property being measured
    • G01N2021/1785Three dimensional
    • G01N2021/1787Tomographic, i.e. computerised reconstruction from projective measurements
    • 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
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N2021/3125Measuring the absorption by excited molecules
    • 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
    • G01N21/59Transmissivity
    • G01N21/5907Densitometers
    • G01N2021/5915Processing scan data in densitometry
    • G01N2021/5953Processing scan data in densitometry for detecting a spatial spectrum
    • 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/75Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
    • G01N21/77Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator
    • G01N2021/7769Measurement method of reaction-produced change in sensor
    • G01N2021/7779Measurement method of reaction-produced change in sensor interferometric

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Animal Behavior & Ethology (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Radiology & Medical Imaging (AREA)
  • Emergency Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Dispersion Chemistry (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

【課題】非侵襲的に分析物質の濃度を推定する分析物質の濃度推定装置を提供する。【解決手段】濃度推定装置は、被検体にOCT信号を照射し、被検体から反射または散乱されたOCT信号を受信するOCT装置と、被検体に分光器信号を照射し、被検体から反射または散乱された分光器信号に基づいてスペクトルを獲得する分光器と、を含む信号測定部と、OCT信号に基づいて所定の波長範囲で各波長ごとの経路長の分布を予測し、予測された経路長の分布及び獲得されたスペクトルに基づいて分析物質の濃度を推定するプロセッサとを含む。

Description

本発明は、非侵襲的に分析物質の濃度を推定する装置及び方法、ならびに信号測定装置に関する。
糖尿病は、多様な合併症を起こす難治性慢性疾患であって、規則的に血糖をチェックして合併症の発生を予防しなければならない。また、インスリンを投与する場合には、低血糖を回避して、インスリン投与量を調節するために血糖をチェックしなければならない。一般的に、血糖を測定するために侵襲的な方式が用いられている。侵襲的に血糖を測定する方式は、測定の信頼性が高いとはいえ、注射を用いるために、血液採取の苦痛や煩わしさがあり、かつ疾病感染の危険が存在する。最近、血液を直接採取せず、散乱分光学(diffuse spectroscopy)に基づいて非侵襲的に分析物の濃度を推定する方法が研究されている。
本発明は、光の散乱による光路長の分布を考慮して非侵襲的に分析物質の濃度を推定する装置及び方法、ならびに信号測定装置を提供することを目的とする。
一態様によれば、分析物質の濃度推定装置は、被検体にOCT(optical coherence tomography)信号を照射し、被検体から反射または散乱されたOCT信号を受信するOCT装置と、被検体に分光器(spectrometer)信号を照射し、被検体から反射または散乱された分光器信号に基づいてスペクトルを獲得する分光器と、を含む信号測定部と、OCT信号に基づいて所定の波長範囲で各波長ごとの経路長の分布を予測し、予測された経路長の分布及び獲得されたスペクトルに基づいて分析物質の濃度を推定するプロセッサと、を含みうる。
OCT装置は、第1光を照射する第1光源と、第1光を第1分割光と第2分割光とに分割するビームスプリッタと、第1分割光を反射させる基準ミラーと、第1光源によって照射されて基準ミラーで反射した第1分割光と被検体によって散乱された第2分割光とによって生成された干渉光を検出してOCT信号に転換する第1光感知部とを含みうる。
分光器は、第2光を照射する第2光源と、第2光源によって照射されて被検体から散乱された第2光を検出してスペクトルを獲得する第2光感知部と、を含みうる。
プロセッサは、波長ごとの散乱係数及びg因子(g-factor)に基づいて、前記波長ごとの経路長の分布を予測しうる。
プロセッサは、OCT信号に基づいて、被検体に対する分光器信号の浸透深さごとの散乱係数及びg因子を獲得しうる。
プロセッサは、散乱係数関数式に基づいて、浸透深さごとの散乱係数から波長ごとの散乱係数を推定しうる。
プロセッサは、モンテカルロシミュレーション(Monte Carlo simulation)を用いて経路長の分布DBを構築しうる。
プロセッサは、推定された波長ごとの散乱係数及びg因子に基づいて、構築された経路長の分布DBから波長ごとの経路長の分布を獲得しうる。
プロセッサは、波長ごとの経路長の分布及びスペクトルに基づいて吸収係数スペクトルを獲得し、獲得された吸収係数スペクトルに基づいて分析物質の濃度を推定しうる。
プロセッサは、波長ごとの経路長の分布と波長ごとの吸収係数とに基づいて獲得されるスペクトルが、被検体から獲得されたスペクトルに収斂するように波長ごとの吸収係数を決定しうる。
プロセッサは、獲得された吸収係数スペクトルから温度を含むノイズ成分を除去しうる。
プロセッサは、主成分分析(Principal Component Analysis)及び特異値分解(Singular Value Decomposition)のうちの少なくとも1つに基づいて、獲得された吸収係数スペクトルからノイズ成分ベクトルを抽出しうる。
プロセッサは、最小二乗法(Least square method)を含んだノイズ除去技法を用いて吸収係数スペクトルからノイズを除去しうる。
分析物質は、グルコース(glucose)、尿素(urea)、乳酸(lactate)、中性脂肪(triglyceride)、総タンパク質(total protein)、コレステロール(cholesterol)及びエタノール(ethanol)のうちの少なくとも1つを含みうる。
一態様によれば、分析物質の濃度推定方法は、被検体からOCT信号を獲得する段階と、被検体からスペクトルを獲得する段階と、OCT信号に基づいて、所定の波長範囲で各波長ごとの経路長の分布を予測する段階と、予測された経路長の分布及びスペクトルに基づいて分析物質の濃度を推定する段階とを含みうる。
経路長の分布を予測する段階は、波長ごとの散乱係数及びg因子に基づいて、波長ごとの経路長の分布を予測しうる。
経路長の分布を予測する段階は、OCT信号に基づいて浸透深さごとの散乱係数及びg因子を獲得する段階を含みうる。
経路長の分布を予測する段階は、散乱係数関数式に基づいて、浸透深さごとの散乱係数から波長ごとの散乱係数を推定する段階をさらに含みうる。
経路長の分布を予測する段階は、モンテカルロシミュレーションを用いて経路長の分布DBを構築する段階をさらに含みうる。
経路長の分布を予測する段階は、波長ごとの散乱係数及びg因子に基づいて、経路長の分布DBから波長ごとの経路長の分布を獲得する段階をさらに含みうる。
分析物質の濃度を推定する段階は、波長ごとの経路長の分布及びスペクトルに基づいて吸収係数スペクトルを獲得し、獲得された吸収係数スペクトルに基づいて分析物質の濃度を推定しうる。
分析物質の濃度を推定する段階は、波長ごとの経路長の分布と波長ごとの吸収係数とに基づいて獲得されるスペクトルを、被検体から獲得されたスペクトルに収斂させる波長ごとの吸収係数を決定することができる。
分析物質の濃度を推定する段階は、獲得された吸収係数スペクトルから温度を含むノイズ成分を除去する段階をさらに含みうる。
ノイズ成分を除去する段階は、主成分分析及び特異値分解のうちの少なくとも1つに基づいてノイズ成分ベクトルを抽出しうる。
ノイズ成分を除去する段階は、最小二乗法を含んだノイズ除去技法を用いて、吸収係数スペクトルからノイズを除去しうる。
一態様によれば、信号測定装置は、第1光を照射する第1光源、第1光を第1分割光と第2分割光とに分割するビームスプリッタ、第1分割光を反射させる基準ミラー、及び第1光源によって照射されて基準ミラーに反射した第1分割光と被検体によって散乱された第2分割光とによって生成された干渉光を検出してOCT信号に転換する第1光感知部を含むOCT装置と、第2光を照射する第2光源と、第2光源によって照射されて被検体から散乱された第2光を検出してスペクトルを獲得する第2光感知部を含む分光器とを含みうる。
OCT装置は、第1光感知部から出力されたOCT(光干渉断層映像)信号を収集するDAQ装置と、DAQ装置に生成されたアナログ信号をデジタル信号に変換するA/D変換部とをさらに含みうる。
本発明によれば、光の散乱による光路長の分布を考慮して非侵襲的に分析物質の濃度を推定することができる。
一実施形態に係る濃度推定装置のブロック図である。 一実施形態に係る信号測定部のブロック図である。 深さごとの散乱係数を抽出することを説明する図面である。 深さ別の散乱係数を抽出することを説明する図面である。 深さごとの経路長の分布を獲得することを説明する図面である。 深さごとの経路長の分布を獲得することを説明する図面である。 深さごとの経路長の分布を獲得することを説明する図面である。 波長ごとの吸収係数を抽出することを説明する図面である。 波長ごとの吸収係数を抽出することを説明する図面である。 抽出された吸収係数スペクトルの例である。 他の実施形態に係る濃度推定装置のブロック図である。 一実施形態に係る濃度推定方法のフローチャートである。 波長ごとの経路長の分布の予測段階(830)の一実施形態である。 濃度推定装置を含む電子装置のブロック図である。 濃度推定装置を含む電子装置の時計型の実施形態を図示する図面である。 濃度推定装置を含む電子装置のモバイル型の実施形態を図示する図面である。 濃度推定装置を含む電子装置のイヤ(Ear)ウェアラブル型の実施形態を図示する図面である。
分析物質の濃度を推定するための実施形態について、以下図面を参照して詳細に説明する。実施形態の優位点及び特徴、ならびに、それらを実現するための手段は、図面とともに実施形態を参照することで、明確になるであろう。明細書の全般にわたって同じ参照符号は、同じ構成要素を意味する。
第1、第2などの用語は、多様な構成要素の説明に使われるが、構成要素は、用語によって限定されるものではない。このような用語は、1つの構成要素を他の構成要素から区別する目的でのみ使われる。単数の表現は、文脈上、特に明示しない限り、複数をも含む。また、ある部分がある構成要素を「含む」と記載されている場合、特記しない限り、他の構成要素を除外するものではなく、他の構成要素をさらに含みうることを意味する。また、明細書に記載の「...部」、「モジュール」などの用語は、少なくとも1つの機能や動作を処理する単位を意味し、これは、ハードウェアまたはソフトウェアとして具現されるか、ハードウェアとソフトウェアとの組み合わせによっても構成可能である。
図1は、一実施形態に係る濃度推定装置のブロック図である。
図1を参照すれば、濃度推定装置100は、信号測定部110及びプロセッサ120を含む。
信号測定部110は、散乱媒質である被検体から深さごと散乱係数とg因子とを獲得するためのOCT信号を測定する第1測定部111と、被検体から散乱または反射したスペクトルを測定する第2測定部112と、を含みうる。第1測定部111と第2測定部112は、互いに結合され、1つのレンズを共有して被検体の同一部位を測定するように構成される。各測定部111、112により測定された信号は、互いに同期されて保存部に保存されるか、プロセッサ120へ出力される。信号測定部110は、プロセッサ120から物理的に分離された独立した別個装置として構成され、分析物質の濃度推定アルゴリズムを含んだ外部機器に有線または無線通信を通じて測定された信号を伝送しうる。第1測定部111は、例えばOCTスキャナーまたはOCT装置であり、第2測定部112は、例えば分光器である。
図2は、一実施形態に係る信号測定部のブロック図である。
図2を参照すれば、第1測定部111は、光を照射する第1光源211、被検体からOCT信号測定のための干渉光を生成する干渉計、干渉光を検出してOCT信号に転換する第1光感知部215、DAQ(data acquisition)装置216、及びA/D変換部217を含みうる。干渉計は、ビームスプリッタ212、レンズ213a、213b、及び基準ミラー214を含みうる。
第1光源211は、OCT信号測定に使われる波長光を照射しうる。OCT信号測定に使われる波長は、散乱係数が吸収係数よりも十分大きい条件を満足する環境でOCT信号が測定されるように、水の吸収が相対的に少ない850nm以下または1300nm以下の領域である。OCT信号の測定に使われる光の波長においては、光の散乱係数は、光の吸収係数よりも現在値だけ大きい。ただし、前述した波長領域に限定されるものではなく、2つ以上の波長領域の光を照射するように構成されうる。
ビームスプリッタ212は、第1光源211によって照射された光(以下、「第1光」とも称する)を分割し、分割された一部の光(以下、「第1分割光」とも称する)を基準ミラー214に入射させ、分割された残りの光(以下、「第2分割光」とも称する)を被検体(SMP)の測定部位に入射させる。ビームスプリッタ212は、第1光源211からの光を損失なしに一定の比率、例えば、50:50に分割することができ、光カプラーを含みうる。ビームスプリッタ212によって分割された第1分割光は、第1レンズ213aを通過して基準ミラー214に入射され、第2分割光は、第2レンズ213bを通過して被検体SMPに入射される。
ビームスプリッタ212によって分割された第1分割光は、基準ミラー214で反射し、第2分割光は、被検体の表面や被検体の内部組織によって透過、散乱または反射というように多様に反応する。基準ミラー214で反射した第1分割光と被検体から反射、透過及び散乱など多様に反応した(以下、「散乱」と通称する)第2分割光は、再びビームスプリッタ212に入射され、ビームスプリッタ212は、入射された第1分割光と第2分割光とから干渉光を生成して第1光感知部215に伝達しうる。
第1光感知部215は、第1分割光と第2分割光とから生成された干渉光を検出し、電気的な信号に変換しうる。第1光感知部215は、1つ以上のピクセルからなり、各ピクセルは、フォトダイオード(Photo Diode)、フォトトランジスタ(Photo Transistor)、フォトゲート(Photogate)、またはピンドフォトダイオード(Pinned Photo Diode)などで構成可能である。OCT信号は、光の浸透深さことの強度プロファイルを含みうる。強度プロファイルは、ピクセルに受光された光の強度(intensity)を示しうる。
DAQ装置216は、第1光感知部215からのOCT信号を収集し、OCT信号を増幅しうる。A/D変換部217は、アナログのOCT信号をデジタル信号に変換して出力しうる。DAQ216は、OCT信号を保存するためのメモリ及び/またはOCT信号を増幅するための増幅器を含みうる。
図2を参照すれば、第2測定部112は、被検体SMPに光(以下、「第2光」とも称する)を照射する第2光源221と、被検体SMPから散乱された第2光を検出し、検出された第2光に基づいてスペクトルを獲得する第2光感知部222と、を含みうる。
第2光源221は、被検体SMPに所定の波長光を照射することができる。第2光源221によって照射された第2光は、第1測定部111の第2レンズ213bを共用して第1測定部111で測定する被検体部位と同じ部位に入射される。第2光源221は、例えば、近赤外線(Near Infrared Ray、NIR)、中赤外線(Mid Infrared Ray、MIR)、レーザ光などを照射することができるが、これに制限されず、光の波長は、測定目的や分析物質の種類などによって変わりうる。第2光源221は、発光ダイオード(light emitting diode、LED)、レーザダイオード(laser diode)、または蛍光体などの発光体で形成されうるが、これらに限定されない。
第2光源221は、必ずしも単一の発光体で構成される必要はなく、複数の発光体の集合で構成してもよい。第2光源221の各発光体は、時分割で順次駆動されるか、同時駆動される。この際、複数の発光体は、それぞれ互いに異なる浸透深さに対するスペクトルを測定するために、互いに異なる波長光を照射するように形成されうる。例えば、短波長である緑色波長の発光体を用いて浸透深さが浅い部位でスペクトルを獲得することができ、相対的に長波長である赤外波長の発光体を用いて浸透深さが深い部位でスペクトルを獲得することができる。ただし、例示された緑色波長、赤外波長に限定されず、青色波長、赤色波長など分析物質、測定部位の深さなどによって多様な波長の発光体で構成することができる。
第2光感知部222は、被検体SMPから散乱されて第2レンズ213bを通じて帰って来る光を検出してスペクトルを獲しうる。第2光感知部222は、フォトダイオード、フォトトランジスタ、CMOS(Complementary Metal-Oxide Semiconductor)イメージセンサー、CCD(Charge-Coupled Device)イメージセンサーなどで形成されうる。第2光感知部222は、必ずしも1つの素子で構成される必要はなく、複数の素子が集まってアレイ形態で構成することもできる。第2光感知部222は、被検体SMPから散乱された光を分光するためにプリズムや回折格子などを含みうる。
第2光源221と第2光感知部222は、一定の距離を有するように所定の間隔で離隔して配置される。第2光源221と第2光感知部222は、互いに異なる2つ以上の光路上のスペクトルを測定できるように互いに異なる2つ以上の隔離距離を有するように配置される。例えば、1つのLEDとそのLEDから互いに異なる距離上に2つ以上のPD(photodiode)とが配置される。または、複数のLEDと各LEDから互いに異なる距離上に1つ以上のPDとが配置される。互いに短い距離(例えば、所定の第1距離よりもさらに短い距離)上に離隔配置されたLEDとPDとの組み合わせを用いて、浸透深さが浅い部位(例えば、毛細血管)からスペクトルを獲得することができ、互いに長い距離(例えば、所定の第2距離よりもさらに長い距離)上に離隔配置されたLEDとPDとの組み合わせを用いて、浸透深さが深い部位(例えば、細動脈)からスペクトルを測定することができる。
再び図1を参照すれば、プロセッサ120は、信号測定部110を制御しうる。プロセッサ120は、信号測定部110からOCT信号とスペクトルとを受信し、受信されたOCT信号とスペクトルとを用いて分析物質の濃度を推定することができる。例えば、プロセッサ120は、OCT信号に基づいて波長ごとの経路長の分布を予測することができる。また、このように予測された波長ごとの経路長の分布を考慮してスペクトルから分析物質の濃度を推定することができる。
以下、図3Aないし図6を参照して、波長ごと経路長の分布を獲得し、分析物質の濃度を推定する実施形態について説明する。
図3A及び図3Bは、深さごとの散乱係数の抽出を説明する図面である。
図3A及び図3Bを参照して、深さごとの散乱係数及びg因子を抽出する実施形態について説明する。
図3Aは、信号測定部110が測定したOCT信号の浸透深さごとの強度プロファイル(intensity profile)、すなわち、A-スキャンデータを例示したものである。ここで、区間Aは、皮膚の表皮(epidermis)区間、区間Bは、表皮と真皮(dremis)との間の中間区間、区間Cは、皮膚の真皮区間を示す。図3Bは、信号測定部110によって測定されたOCT信号のB-スキャンデータを例示したものである。
プロセッサ120は、このように信号測定部110によって測定されたOCT信号の浸透深さごとの強度プロファイルに基づいて、各浸透深さごとの散乱係数及び/またはg因子を抽出しうる。例えば、プロセッサ120は、ヘテロダイン効率方程式(Heterodyne efficiecy equation)を用いて浸透深さごとに散乱係数及びg因子を抽出することができる。プロセッサ120は、特定の深さに対してヘテロダイン効率方程式に特定の範囲内の任意の散乱係数とg因子とを入力し、ヘテロダイン効率方程式の出力値と信号測定部110によって測定された強度プロファイルとの間の誤差を最小にする散乱係数とg因子とを抽出することができる。ここで、g因子は、異方性係数(anisotropic factor)を示し、散乱角θのコサインの平均であって、単一散乱方向を媒介変数として使用する確率密度を有する値を示す。プロセッサ120は、ヘテロダイン効率方程式(下記数1及び数5)を信号測定部110によって検出された光に適用することができる。
Figure 2022061939000001
ここで、I(z)は、信号測定部110によって測定された浸透深さごとの光の強度プロファイルで各浸透深さ(z)に対する強度を意味する。Iは、初期入射光の強度を示す。
Figure 2022061939000002
は、特定の範囲に属す散乱係数
Figure 2022061939000003
とg因子(g)とを用いて浸透深さ(z)に対してモデリングしたプロファイルを示す。ここで、
Figure 2022061939000004
は、ヘテロダイン効率因子(Heterodyne efficiency factor)であって、下記の数式を通じて求めうる。
Figure 2022061939000005
ここで、A、B及びfは、レンズの特性であって、Aは、1に設定しうる。Bは、浸透深さ(z)をサンプルの屈折率(n)で割った値をレンズとサンプルとの距離に加えて算出される。fは、レンズの焦点距離を示す。zは、深さ、nは、サンプルの屈折率を示し、w0は、レファレンスの1/e強度半径(intensity radius)を示す。θrmsは、g=cosθrmsの関係式を通じて求めうる。ここで、散乱係数
Figure 2022061939000006
とg因子(g)は、特定の範囲内の散乱係数とg因子とのうちから入力された散乱係数とg因子とを示す。
プロセッサ120は、このように特定の範囲内で散乱係数とg因子とを変更しながら、初期入射光(I0)を用いてプロファイルをモデリングし、モデリングされたプロファイル
Figure 2022061939000007
を測定されたプロファイル(I(z))に収斂させる散乱係数とg因子とを浸透深さ(z)に対する散乱係数とg因子として抽出することができる。このような方式であらゆる深さに対して散乱係数とg因子とを求めうる。
プロセッサ120は、このように浸透深さごとの散乱係数とg因子とが抽出されれば、抽出された浸透深さごとの散乱係数とg因子とを用いて波長ごとの経路長の分布を獲得することができる。例えば、深さごとの散乱係数とg因子とをミー散乱方程式(Mie scattering equation)を用いて波長ごとの散乱係数を獲得することができる。
図4Aないし図4Cは、深さごとの経路長の分布の獲得を説明する図面である。
プロセッサ120は、波長(λ1)から獲得されたOCT信号を用いて浸透深さ(z)に対する散乱係数(μ,1)が抽出されれば、波長(λ1)と散乱係数(μ,1)とを用いてミー散乱方程式を通じて、下記の数式のような関係式を導出することができ、この関係式は、図4Aのグラフのように表現される。プロセッサ120は、このように導出された関係式を用いて求めようとする波長(λ2)に対する散乱係数(μ,2)を求めうる。
Figure 2022061939000008
ここで、
Figure 2022061939000009
は、求めようとする波長(λ)に対する散乱係数を示す。Bは、被検体によってあらかじめ設定される値であって、例えば、-1.5に設定しうる。Aは、OCT信号の波長(λ1)と散乱係数(μ,1)とを通じて求められる値である。この際、OCT信号を2つ以上の波長で測定する場合、より正確な関係式を導出することができる。
また、プロセッサ120は、特定の範囲内の散乱係数とg因子とを用いてモンテカルロシミュレーションを通じて波長ごとの経路長の分布DBを構築することができる。図4Bは、光源と光感知部との第1距離に対してモンテカルロシミュレーションを通じて獲得した経路長の分布(1)と、第1距離と異なる光源と光感知部との第2距離に対してモンテカルロシミュレーションを通じて獲得した経路長の分布(2)と、を例示したものである。プロセッサ120は、モンテカルロシミュレーションを通じて経路長の分布DBが構築されれば、波長ごとの散乱係数及びg因子を用いて経路長の分布DBから波長ごとの経路長の分布を抽出することができる。図4Cは、このような方式で1550nm、1600nm、1650nmに対して抽出された経路長の分布を例示したものである。
図5Aないし図5Bは、波長ごとの吸収係数の抽出を説明する図面である。図6は、抽出された吸収係数スペクトルの例である。
図5Aは、一般的に光源によって照射された光が被検体SMPの組織特性によって多様な浸透深さを通過する複数の経路長の分布を有することを示す。図5Bは、光源と光感知部との距離による散乱スペクトルを例示したものである。示したように、光源と光感知部との距離(r)、浸透深さ(z)、入射光の強度などによって生体組織から吸収される吸光度(A)が異なり、これは、多様な光路長によって散乱スペクトル(I(r))が変わりうるということが分かる。
プロセッサ120は、測定されたスペクトルを用いて分析物質の濃度推定のための吸収係数スペクトルを獲得することができ、この際、図5A及び図5Bを通じて分かるように、多様な光路長の分布によってスペクトルの強度が変わりうるので、波長別の経路長の分布を考慮して波長ごとの吸収係数を獲得しうる。
例えば、下記数式は、波長ごとの吸収係数を算出する式の一例であって、プロセッサ120は、下記数式を満足させる吸収係数を各波長ごとに獲得することができる。すなわち、プロセッサ120は、各波長ごとに特定の範囲内の吸収係数を用いて獲得したスペクトルを被検体SMPから測定されたスペクトルに収斂させる波長ごとの吸収係数を求めうる。図6は、各波長ごとに特定の範囲(LB、UB)内の吸収係数を用いて獲得した吸収係数スペクトル(AS)を例示している。
Figure 2022061939000010
ここで、I(λ)は、測定されたスペクトルの波長(λ)での強度を示す。Iは、入射光の強度を示す。
Figure 2022061939000011
は、波長(λ)に対する経路長の分布を示す。
Figure 2022061939000012
は、経路長を示す。
Figure 2022061939000013
は、求めようとする波長(λ)に対する吸収係数を示す。
プロセッサ120は、吸収係数スペクトルから温度を含むノイズ成分を除去しうる。例えば、プロセッサ120は、主成分分析及び/または特異値分解技法を用いてノイズ成分ベクトルを抽出し、最小二乗法を含むノイズ除去方法を用いてノイズ成分ベクトルに基づいて吸収係数スペクトルからノイズを除去しうる。
プロセッサ120は、このように吸収係数スペクトルが獲得されれば、吸収係数スペクトルを用いて分析物質の濃度を推定することができる。この際、分析物質は、グルコース、尿素、尿酸(Uric acid)、乳酸、中性脂肪、タンパク質(protein)、コレステロール、及びエタノールなどを含みうるが、これらに限定されるものではない。プロセッサ120は、吸収係数スペクトルを用いて線形回帰分析、PLS(Partial Least Square)、CLS(Classical Least Square)などを用いて分析物質の濃度を推定することができる。この際、波長ごとの吸収係数と分析物質の濃度との相関関係を定義した濃度推定モデルが線形回帰分析などを用いてあらかじめ定義される。
上述の実施形態によれば、散乱媒質で経路長の分布を考慮して波長ごとの吸収係数を獲得して分析物質の濃度を推定することにより、あらゆる波長の経路長の分布が同一であると仮定して分析する一般的なランベルト・ベールの法則(Beer-Lambert Law)を用いた解析に基づく分析物質の濃度推定に比べて正確性を向上させうる。
図7は、他の実施形態に係る濃度推定装置のブロック図である。
図7を参照すれば、濃度推定装置700は、信号測定部710、プロセッサ720、出力部730、保存部740、及び通信部750を含みうる。信号測定部710及びプロセッサ720は、図1の信号測定部110及びプロセッサ120と同じ構成なので、詳しい説明は省略する。
出力部730は、プロセッサ720の処理結果を使用者に供給することができる。一例として、出力部730は、ディスプレイを含み、ディスプレイを通じて分析物質、例えば、血糖推定値を表示することができる。この際、血糖推定値が正常範囲を外れれば、使用者が容易に認識できるように色や線の太さなどを調節するか、正常範囲を併せて表示することにより、使用者に警告情報を提供しうる。また、出力部730は、スピーカー、ハプティックモジュール(haptic module)などの非視覚的出力モジュールを含み、視覚的表示と共に、または単独で音声出力モジュール、ハプティックモジュールを用いて音声、振動、触感などの非視覚的な方式で使用者に推定結果を提供しうる。
保存部740は、分析物質の濃度の推定に必要な基準情報、信号測定部710及び/またはプロセッサ720の処理結果を保存しうる。この際、基準情報は、使用者の年齢、性別、健康状態のような使用者特性情報を含みうる。また、基準情報は、濃度推定モデルを含みうる。また、基準情報は、散乱係数、g因子獲得、モンテカルロシミュレーション、波長ごとの経路長の分布を獲得するためのアルゴリズムなどの情報を含みうる。
保存部740は、フラッシュメモリタイプ(flash memory type)、ハードディスクタイプ(hard disk type)、マルチメディアカードマイクロタイプ(multimedia card micro type)、カードタイプのメモリ(例えば、SDまたはXDメモリなど)、RAM(Random Access Memory)、SRAM(Static Random Access Memory)、ROM(Read-Only Memory)、EEPROM(Electrically Erasable Programmable Read-Only Memory)、PROM(Programmable Read-Only Memory)、磁気メモリ、磁気ディスク、光ディスクなどの記録媒体を含みうるが、これらに限定されるものではない。
通信部750は、外部機器と通信して分析物質の濃度推定と関連した各種のデータを送受信することができる。外部機器は、スマートフォン、タブレットPC、デスクトップPC、ノート型パソコンなどの情報処理装置を含みうる。例えば、血糖推定の結果を使用者のスマートフォンなどに送信して使用者が相対的に性能に優れた機器を通じて使用者の血糖を管理及びモニタリングさせうる。ただし、これらに限定されるものではない。通信部750は、ブルートゥース(Bluetooth)通信、BLE(Bluetooth Low Energy)通信、近距離無線通信(Near Field Communication、NFC)、WLAN通信、ジグビー(Zigbee)通信、赤外線(Infrared Data Association、IrDA)通信、WFD(Wi-Fi Direct)通信、UWB(ultra-wideband)通信、Ant+通信、Wi-Fi通信、RFID(Radio Frequency Identification)通信、3G通信、4G通信及び5G通信などを含む多様な有線または無線通信技術を用いて外部機器と通信することができる。ただし、これらに限定されるものではない。
図8は、一実施形態に係る濃度推定方法のフローチャートである。
図8は、図1及び図7の濃度推定装置100、700によって実行される濃度推定方法の一実施形態である。重複説明を避けるために、既述した内容については、簡略化して以下説明する。
まず、濃度推定装置は、第1測定部を通じて被検体SMPからOCT信号を獲得し(810)、第2測定部を通じてOCT信号が獲得された被検体SMPの同じ部位に光を照射し、被検体からの散乱光を検出してスペクトルを獲得しうる(820)。
次いで、濃度推定装置は、OCT信号を用いて波長ごとの経路長の分布を予測しうる(830)。
図9は、波長ごとの経路長の分布の予測段階(830)の一実施形態である。
図9を参照して、波長ごとの経路長の分布を予測する段階(830)の一実施形態について説明する。
濃度推定装置は、OCT信号を用いて浸透深さごとの散乱係数及びg因子を抽出しうる(910)。浸透深さごとの散乱係数及びg因子は、数式上記数1及び数5に例示されたように、ヘテロダイン効率方程式を用いて獲得することができる。例えば、特定の範囲内の散乱係数とg因子とを変更しながら、プロファイルをモデリングし、モデリングされたプロファイルを測定された強度プロファイルに収斂させる散乱係数とg因子とを抽出することができる。
次いで、段階(910)において獲得された浸透深さごとの散乱係数に基づいて波長ごとの散乱係数を獲得しうる(920)。例えば、ミー散乱方程式に基づいて導出された数8を通じて波長ごとの散乱係数を獲得することができる。
次いで、特定の範囲の散乱係数及びg因子を用いてモンテカルロシミュレーションを用いて経路長の分布DBを構築しうる(930)。波長ごとの散乱係数を獲得する段階(920)及び経路長の分布DBを構築する段階(930)は、先後関係が存在するものではなく、何れか1つが先行されるか、同時に行われうる。
次いで、波長ごとの散乱係数及びg因子を用いて経路長の分布DBから各波長ごとの経路長の分布を抽出しうる(940)。
再び図8を参照すれば、段階(830)で波長ごとの経路長の分布が予測されれば、波長ごとの経路長の分布及び被検体SMPから獲得されたスペクトルに基づいて分析物質の濃度を推定しうる(840)。濃度推定装置は、波長ごとの経路長の分布を用いてスペクトルから波長ごとの吸収係数を獲得することができる。この際、上記の数10のように、波長ごとの経路長の分布と求めようとする吸収係数を用いて獲得したスペクトルと段階(820)から測定されたスペクトルとを最小にする吸収係数を波長ごとに求めうる。このように吸収係数スペクトルが獲得されれば、温度のようなノイズベクトルを抽出して、吸収係数スペクトルからノイズを除去することができ、ノイズが除去された吸収係数スペクトルを用いて線形回帰分析などを通じて分析物質の濃度を推定することができる。
図10は、濃度推定装置を含む電子装置のブロック図である。
図10を参照すれば、電子装置1001は、プロセッサ1020、メモリ1030、入力装置1050、音響出力装置1055、表示装置1060、オーディオモジュール1070、センサーモジュール1076、インターフェース1077、ハプティックモジュール1079、カメラモジュール1080、電力管理モジュール1088、バッテリ1089、通信モジュール1090、加入者識別モジュール1096、及び/またはアンテナモジュール1097を含みうる。この構成要素の一部は、省略され、他の構成要素が追加されうる。
濃度推定装置100、700は、統合された回路として構成されて、電子装置1001のセンサーモジュール1076に搭載されるか、互いに異なる構成要素に分散される。例えば、濃度推定装置100、700の信号測定部110、710は、センサーモジュール1076に含まれ、プロセッサ120、720は、プロセッサ1020に含まれうる。また、出力部730は、音響出力装置1055、表示装置1060、オーディオモジュール1070などに分散され、保存部740は、メモリ1030として構成され、通信部750は、通信モジュール1090に含まれうる。
プロセッサ1020は、プログラム1040などを実行してプロセッサ1020に連結された電子装置1001の構成要素を制御し、多様なデータ処理または演算を行いうる。例えば、プロセッサ1020は、データ処理または演算の一例として他の構成要素、例えば、センサーモジュール1076や通信モジュール1090などから受信された命令及び/またはデータを揮発性メモリ1032にロードし、揮発性メモリ1032に保存された命令及び/またはデータを処理して、結果データを不揮発性メモリ1034に保存することができる。
プロセッサ1020は、メインプロセッサ1021(例:中央処理装置、アプリケーションプロセッサなど)及びこれと独立して、または共に運用可能な補助プロセッサ1023(例:グラフィック処理装置、イメージシグナルプロセッサ、センサーハブプロセッサ、コミュニケーションプロセッサなど)を含みうる。補助プロセッサ1023は、メインプロセッサ1021よりも電力を少なく使用し、特化された機能を実行することができる。補助プロセッサ1023は、メインプロセッサ1021がインアクティブ状態(スリープ状態)にある間にメインプロセッサ1021の代わりに、またはメインプロセッサ1021がアクティブ状態(アプリケーション実行状態)にある間にメインプロセッサ1021と共に電子装置1001の構成要素のうち、一部構成要素、例えば、表示装置1060、センサーモジュール1076、通信モジュール1090などと関連した機能及び/または状態を制御しうる。補助プロセッサ1023、例えば、イメージシグナルプロセッサやコミュニケーションプロセッサなどは、機能的に関連がある他の構成要素、例えば、カメラモジュール1080や通信モジュール1090などの一部として構成されることもある。
プロセッサ1020は、分析物質の濃度推定に対する使用者の要請に応じて濃度推定装置100、700に制御信号を伝達しうる。
メモリ1030は、電子装置1001の構成要素が必要とする多様なデータ、例えば、ソフトウェア、及びこれと関連した命令に対する入力データ及び/または出力データなどを保存しうる。メモリ1030は、揮発性メモリ1032及び/または不揮発性メモリ1034を含みうる。
プログラム1040は、メモリ1030にソフトウェアとして保存され、オペレーションシステム(OS)1042、ミドルウェア1044及び/またはアプリケーション1046を含みうる。
入力装置1050は、電子装置1001の構成要素に使われる命令及び/またはデータを使用者などから受信しうる。入力装置1050は、マイク、マウス、キーボード、及び/またはデジタルペン(スタイラスペンなど)を含みうる。
音響出力装置1055は、音響信号を電子装置1001の外部に出力しうる。音響出力装置1055は、スピーカー及び/またはレシーバーを含みうる。スピーカーは、マルチメディア再生または録音再生のように一般的な用途として使われ、レシーバーは、着信電話を受信するために使われる。レシーバーは、スピーカーの一部として結合されているか、または独立した別途の装置として構成可能である。
表示装置1060は、電子装置1001の外部に情報を視覚的に提供しうる。表示装置1060は、ディスプレイ、ホログラム装置、またはプロジェクター及びこれらの装置を制御するための制御回路を含みうる。表示装置1060は、タッチを感知するように設定されたタッチ回路(Touch Circuitry)及び/またはタッチによって発生する力の強度を測定するように設定されたセンサー回路(圧力センサーなど)を含みうる。
オーディオモジュール1070は、音を電気信号に変換させるか、逆に電気信号を音に変換させうる。オーディオモジュール1070は、入力装置1050を通じて音を獲得するか、音響出力装置1055及び/または電子装置1001と直接または無線で連結された他の電子装置1002、1004のスピーカー及び/またはヘッドホンを通じて音を出力することができる。
センサーモジュール1076は、電子装置1001の作動状態(電力、温度など)、または外部の環境状態(使用者状態など)を感知し、感知された状態に対応する電気信号及び/またはデータ値を生成しうる。センサーモジュール1076は、ジェスチャセンサー、ジャイロセンサー、気圧センサー、マグネチックセンサー、加速度センサー、グリップセンサー、近接センサー、カラーセンサー、IR(Infrared)センサー、生体センサー、温度センサー、湿度センサー、及び/または照度センサーを含みうる。濃度推定装置100、700は、センサーモジュール1076に含まれた生体センサーの1つである。
インターフェース1077は、電子装置1001が他の電子装置1002、1004と直接または無線で連結されるために使われる1つ以上の指定されたプロトコルを支援しうる。インターフェース1077は、HDMI(High Definition Multimedia Interface)、USB(Universal Serial Bus)インターフェース、SDカードインターフェース、及び/またはオーディオインターフェースを含みうる。
連結端子1078は、電子装置1001が他の電子装置1002、1004と物理的に連結されるコネクタを含みうる。連結端子1078は、HDMIコネクタ、USBコネクタ、SDカードコネクタ、及び/またはオーディオコネクタ(ヘッドホンコネクタなど)を含みうる。
ハプティックモジュール1079は、電気的信号を使用者が触覚または運動感覚を通じて認知することができる機械的な刺激(振動、動きなど)または電気的な刺激に変換しうる。ハプティックモジュール1079は、モータ、圧電素子、及び/または電気刺激装置を含みうる。
カメラモジュール1080は、静止画像及び動画を撮影しうる。カメラモジュール1080は、1つ以上のレンズを含むレンズアセンブリー、イメージセンサー、イメージシグナルプロセッサ、及び/またはフラッシュを含みうる。カメラモジュール1080に含まれたレンズアセンブリーは、イメージ撮影の対象である被写体からの放出光を収集することができる。
電力管理モジュール1088は、電子装置1001に供給される電力を管理しうる。電力管理モジュール1088は、PMIC(Power Management Integrated Circuit)の一部として構成可能である。
バッテリ1089は、電子装置1001の構成要素に電力を供給しうる。バッテリ1089は、再充電不可能な1次電池、再充電可能な2次電池及び/または燃料電池を含みうる。
通信モジュール1090は、電子装置1001とネットワーク環境1000内にある他の電子装置1002、1004やサーバ1008との間の直接(有線)通信チャンネル及び/または無線通信チャンネルの確立、及び確立された通信チャンネルを通じた通信遂行を支援しうる。通信モジュール1090は、プロセッサ1020から独立して運用され、直接通信及び/または無線通信を支援する1つ以上のコミュニケーションプロセッサを含みうる。通信モジュール1090は、無線通信モジュール1092(例:セルラー通信モジュール、近距離無線通信モジュール、GNSS(Global Navigation Satellite Systemなど)通信モジュール)及び/または有線通信モジュール1094(例:LAN(Local Area Network)通信モジュール、電力線通信モジュールなど)を含みうる。通信モジュール1090は、第1ネットワーク1098(例:ブルートゥース、Wi-Fi DirectまたはIrDA(Infrared Data Association)のような近距離通信ネットワーク)または第2ネットワーク1099(例:セルラーネットワーク、インターネット、またはコンピュータネットワーク(LAN、WANなど)のような遠距離通信ネットワーク)を通じて他の電子装置と通信することができる。このように多種の通信モジュール1090は、単一チップなどで統合されるか、互いに個別の複数チップとして構成可能である。無線通信モジュール1092は、加入者識別モジュール1096に保存された加入者情報(例:国際モバイル加入者識別子(IMSI)など)を用いて第1ネットワーク1098及び/または第2ネットワーク1099のような通信ネットワーク内で電子装置1001を確認及び認証することができる。
アンテナモジュール1097は、信号及び/または電力を外部に送信するか、外部から受信しうる。アンテナは、基板(例:PCBなど)上に形成された導電性パターンからなる放射体を含みうる。アンテナモジュール1097は、1つまたは複数のアンテナを含みうる。複数のアンテナが含まれた場合、通信モジュール1090によって複数のアンテナのうちから第1ネットワーク1098及び/または第2ネットワーク1099のような通信ネットワークで使われる通信方式に適したアンテナが選択されうる。選択されたアンテナを通じて通信モジュール1090と他の電子装置との間に信号及び/または電力が送信または受信される。アンテナの以外に他の部品(例:RFICなど)がアンテナモジュール1097の一部として含まれうる。
構成要素の一部は、周辺機器間の通信方式(バス、GPIO(General Purpose Input and Output)、SPI(Serial Peripheral Interface)、MIPI(Mobile Industry Processor Interface)など)を通じて互いに連結され、命令やデータなどを相互に交換しうる。命令またはデータは、第2ネットワーク1099に連結されたサーバ1008を通じて電子装置1001と外部の電子装置1004との間で送信または受信される。他の電子装置1002、1004は、電子装置1001と同一または他種の装置である。電子装置1001で実行される動作の全部または一部は、他の装置1002、1004、1008のうち、1つ以上の装置で実行されうる。例えば、電子装置1001が所定の機能やサービスの実行が要求されている時、当該機能またはサービスを実行する代わりに、1つ以上の他の電子装置に、その機能またはそのサービスの一部または全体の実行を要請することができる。要請を受信した1つ以上の他の電子装置は、要請と関連した追加機能またはサービスを実行し、その実行の結果を電子装置1001に伝達することができる。そのために、クラウドコンピュータ、分散コンピューティング、及び/またはクライアント-サーバコンピューティング技術が用いられる。
図11ないし図13は、濃度推定装置が搭載された電子装置を例示的に示す図面である。
図11を参照すれば、電子装置1001は、時計タイプのウェアラブル装置1001aとして構成され、本体と手首ストラップとを含みうる。本体の前面には、ディスプレイが設けられて、時間情報、受信メッセージ情報などを含む多様なアプリケーション画面が表示される。本体の裏面には、濃度推定装置100、700が配されて、時計タイプのウェアラブル装置1001aが使用者の手首に着用される時、使用者の手首に光を照射し、手首から散乱または反射した光に基づいて血糖などの分析物質に対する濃度を推定しうる。
図12を参照すれば、電子装置1001は、スマートフォン(Smart Phone)のようなモバイル装置1001bとして実現可能である。
モバイル装置1001bは、ハウジング及びディスプレイパネルを含みうる。ハウジングは、電子装置1001bの外觀を形成しうる。ハウジングの第1面には、ディスプレイパネル及びカバーグラス(Cover Glass)が順に配置され、ディスプレイパネルは、カバーグラスを通じて外部に露出される。ハウジングの第2面には、濃度推定装置100、700、カメラモジュール、及び/または赤外線センサーなどが配置される。使用者が電子装置1001bに搭載されたアプリケーションなどを実行して分析物質の濃度情報を要請する場合、濃度推定装置100を用いて分析物質の濃度を推定し、使用者に濃度情報を映像及び/または音響として提供しうる。
図13を参照すれば、電子装置1001は、イヤウェアラブル装置1001cとしても構成可能である。
イヤウェアラブル装置1001cは、本体とイヤストラップ(Ear Strap)とを含みうる。使用者は、イヤストラップを耳殻に引っ掛けて電子装置1001cを着用することができる。イヤストラップは、イヤウェアラブル装置1001cの形態によって省略可能である。本体は、使用者の外耳道(External Auditory Meatus)に挿入される。本体には、濃度推定装置100、700が搭載される。電子装置1001cは、濃度推定の結果を使用者に音響として提供するか、本体内部に設けられた通信モジュールを通じて外部機器、例えば、モバイル、タブレット、PCなどに伝送しうる。
一方、本実施形態は、コンピュータで読み取り可能な記録媒体にコンピュータで読み取り可能なコードとして構成することが可能である。コンピュータで読み取り可能な記録媒体は、コンピュータシステムによって読み取れるデータが保存されるあらゆる種類の記録装置を含む。
コンピュータで読み取り可能な記録媒体の例としては、ROM、RAM、CD-ROM、磁気テープ、フロッピーディスク、光データ保存装置などがあり、また、キャリアウェーブ(例えば、インターネットを介した伝送)の形態で具現するものを含む。また、コンピュータで読み取り可能な記録媒体は、ネットワークで連結されたコンピュータシステムに分散されて、分散方式でコンピュータで読み取り可能なコードとして保存されて実行可能である。そして、本実施形態を具現するための機能的な(functional)プログラム、コード及びコードセグメントは、当該技術分野のプログラマーによって容易に推論される。
当業者ならば、実施形態の技術的思想や必須な特徴を変更せずとも、他の具体的な形態で実施可能であることを理解できるであろう。したがって、前述した実施形態は、あらゆる面で例示的なものであり、限定的ではないということを理解せねばならない。
100 濃度推定装置、
110 信号測定部、
111 第1測定部、
112 第2測定部、
120 プロセッサ、
211 第1光源、
212 ビームスプリッタ、
213a、213b レンズ、
214 基準ミラー、
215 光感知部、
221 第2光源、
222 分光部、
700 濃度推定装置、
710 信号測定部、
720 プロセッサ、
730 出力部、
740 保存部、
750 通信部。

Claims (27)

  1. 被検体にOCT信号を照射し、被検体から反射または散乱されたOCT信号を受信するOCT装置と、被検体に分光器信号を照射し、前記被検体から反射または散乱された分光器信号に基づいてスペクトルを獲得する分光器と、を含む信号測定部と、
    前記OCT信号に基づいて所定の波長範囲で各波長ごとの経路長の分布を予測し、予測された経路長の分布及び前記獲得されたスペクトルに基づいて分析物質の濃度を推定するプロセッサと、
    を含む、分析物質の濃度推定装置。
  2. 前記OCT装置は、
    第1光を照射する第1光源と、
    前記第1光を第1分割光と第2分割光とに分割するビームスプリッタと、
    前記第1分割光を反射させる基準ミラーと、
    基準ミラーで反射した第1分割光と被検体によって散乱された第2分割光とによって生成された干渉光を検出してOCT信号に転換する第1光感知部と、
    を含む、請求項1に記載の分析物質の濃度推定装置。
  3. 前記分光器は、
    第2光を照射する第2光源と、
    第2光源によって照射されて被検体から散乱された第2光を検出してスペクトルを獲得する第2光感知部と、
    を含む、請求項1に記載の分析物質の濃度推定装置。
  4. 前記プロセッサは、
    波長ごとの散乱係数及びg因子に基づいて、前記波長ごとの経路長の分布を予測する、請求項1に記載の分析物質の濃度推定装置。
  5. 前記プロセッサは、
    前記OCT信号に基づいて、前記被検体に対する前記分光器信号の浸透深さごとの散乱係数及び前記g因子を獲得する、請求項4に記載の分析物質の濃度推定装置。
  6. 前記プロセッサは、
    散乱係数関数式に基づいて、前記浸透深さごとの散乱係数を推定する、請求項5に記載の分析物質の濃度推定装置。
  7. 前記プロセッサは、
    モンテカルロシミュレーションを用いて経路長の分布DBを構築する、請求項5に記載の分析物質の濃度推定装置。
  8. 前記プロセッサは、
    前記推定された波長ごとの散乱係数及びg因子に基づいて、前記構築された経路長の分布DBから前記波長ごとの経路長の分布を獲得する、請求項7に記載の分析物質の濃度推定装置。
  9. 前記プロセッサは、
    前記波長ごとの経路長の分布及び前記スペクトルに基づいて吸収係数スペクトルを獲得し、該獲得された吸収係数スペクトルに基づいて分析物質の濃度を推定する、請求項1に記載の分析物質の濃度推定装置。
  10. 前記プロセッサは、
    前記波長ごとの経路長の分布と波長ごとの吸収係数とに基づいて獲得されるスペクトルが、前記被検体から獲得されたスペクトルに収斂するように、前記波長ごとの吸収係数を決定する、請求項9に記載の分析物質の濃度推定装置。
  11. 前記プロセッサは、
    前記獲得された吸収係数スペクトルから温度を含むノイズ成分を除去する、請求項9に記載の分析物質の濃度推定装置。
  12. 前記プロセッサは、
    主成分分析及び特異値分解のうちの少なくとも1つに基づいて、前記獲得された吸収係数スペクトルから前記ノイズ成分ベクトルを抽出する、請求項11に記載の分析物質の濃度推定装置。
  13. 前記プロセッサは、
    最小二乗法を含んだノイズ除去技法を用いて、前記吸収係数スペクトルからノイズを除去する、請求項11に記載の分析物質の濃度推定装置。
  14. 前記分析物質は、
    グルコース、尿素、乳酸、中性脂肪、総タンパク質、コレステロール及びエタノールのうちの少なくとも1つを含む、請求項1に記載の分析物質の濃度推定装置。
  15. 被検体からOCT信号を獲得する段階と、
    被検体からスペクトルを獲得する段階と、
    前記OCT信号に基づいて、所定の波長範囲で各波長ごとの経路長の分布を予測する段階と、
    前記予測された経路長の分布及び前記スペクトルに基づいて分析物質の濃度を推定する段階と、
    を含む、分析物質の濃度推定方法。
  16. 前記経路長の分布を予測する段階は、
    波長ごとの散乱係数及びg因子に基づいて、前記波長ごとの経路長の分布を予測する、請求項15に記載の分析物質の濃度推定方法。
  17. 前記経路長の分布を予測する段階は、
    前記OCT信号に基づいて浸透深さごとの散乱係数及び前記g因子を獲得する段階を含む、請求項16に記載の分析物質の濃度推定方法。
  18. 前記経路長の分布を予測する段階は、
    散乱係数関数式に基づいて、前記浸透深さごとの散乱係数から前記波長ごとの散乱係数を推定する段階をさらに含む、請求項17に記載の分析物質の濃度推定方法。
  19. 前記経路長の分布を予測する段階は、
    モンテカルロシミュレーションを用いて経路長の分布DBを構築する段階をさらに含む、請求項17に記載の分析物質の濃度推定方法。
  20. 前記経路長の分布を予測する段階は、
    前記波長ごとの散乱係数及びg因子に基づいて、前記経路長の分布DBから前記波長ごとの経路長の分布を獲得する段階をさらに含む、請求項19に記載の分析物質の濃度推定方法。
  21. 前記分析物質の濃度を推定する段階は、
    前記波長ごとの経路長の分布及び前記スペクトルに基づいて吸収係数スペクトルを獲得し、該獲得された吸収係数スペクトルに基づいて分析物質の濃度を推定する、請求項15に記載の分析物質の濃度推定方法。
  22. 前記分析物質の濃度を推定する段階は、
    前記波長ごとの経路長の分布と波長ごとの吸収係数とに基づいて獲得されたスペクトルを、前記被検体から獲得されたスペクトルに収斂させる前記波長ごとの吸収係数を決定する、請求項21に記載の分析物質の濃度推定方法。
  23. 前記分析物質の濃度を推定する段階は、
    前記獲得された吸収係数スペクトルから温度を含むノイズ成分を除去する段階をさらに含む、請求項21に記載の分析物質の濃度推定方法。
  24. 前記ノイズ成分を除去する段階は、
    主成分分析及び特異値分解のうちの少なくとも1つに基づいて、前記ノイズ成分ベクトルを抽出する、請求項23に記載の分析物質の濃度推定方法。
  25. 前記ノイズ成分を除去する段階は、
    最小二乗法を含んだノイズ除去技法を用いて、前記吸収係数スペクトルからノイズを除去する、請求項23に記載の分析物質の濃度推定方法。
  26. 第1光を照射する第1光源、前記第1光を第1分割光と第2分割光とに分割するビームスプリッタ、前記第1分割光を反射する基準ミラー、及び第1光源によって照射されて基準ミラーで反射した第1分割光と被検体によって散乱された第2分割光とによって生成された干渉光を検出してOCT信号に転換する第1光感知部を含むOCT装置と、
    第2光を照射する第2光源と、第2光源によって照射されて被検体から散乱された第2光を検出してスペクトルを獲得する第2光感知部を含む分光器と、
    を含む、信号測定装置。
  27. 前記OCT装置は、
    前記第1光感知部から出力されたOCT信号を収集するDAQ装置と、
    前記DAQ装置に生成されたアナログ信号をデジタル信号に変換するA/D変換部と、
    をさらに含む、請求項26に記載の信号測定装置。
JP2021084569A 2020-10-07 2021-05-19 分析物質の濃度推定装置及び方法、ならびに信号測定装置 Pending JP2022061939A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR10-2020-0129269 2020-10-07
KR1020200129269A KR20220046168A (ko) 2020-10-07 2020-10-07 분석 물질의 농도 추정 장치 및 방법과, 신호 측정 장치

Publications (1)

Publication Number Publication Date
JP2022061939A true JP2022061939A (ja) 2022-04-19

Family

ID=77207105

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021084569A Pending JP2022061939A (ja) 2020-10-07 2021-05-19 分析物質の濃度推定装置及び方法、ならびに信号測定装置

Country Status (5)

Country Link
US (1) US20220104735A1 (ja)
EP (1) EP3981329A1 (ja)
JP (1) JP2022061939A (ja)
KR (1) KR20220046168A (ja)
CN (1) CN114295562A (ja)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3588061B1 (en) * 2014-12-23 2023-04-19 Apple Inc. Optical inspection system and method including accounting for variations of optical path length within a sample
KR102615932B1 (ko) 2015-09-01 2023-12-21 애플 인크. 물질의 비접촉 감지를 위한 레퍼런스 스위치 아키텍처
CN114719977A (zh) 2016-04-21 2022-07-08 苹果公司 用于参考切换的光学系统
WO2019160949A1 (en) 2018-02-13 2019-08-22 Masseta Technologies Llc Integrated photonics device having integrated edge outcouplers
US12059239B2 (en) 2018-05-08 2024-08-13 Know Labs, Inc. Electromagnetic shielding in non-invasive analyte sensors
US11903689B2 (en) 2019-12-20 2024-02-20 Know Labs, Inc. Non-invasive analyte sensor device
US20210174960A1 (en) * 2019-12-09 2021-06-10 Medtronic Minimed, Inc. Generative modeling methods and systems for simulating sensor measurements
US12023151B2 (en) 2020-02-20 2024-07-02 Know Labs, Inc. Non-invasive analyte sensing and notification system with decoupled transmit and receive antennas
US11832926B2 (en) 2020-02-20 2023-12-05 Know Labs, Inc. Non-invasive detection of an analyte and notification of results
US12089927B2 (en) 2020-02-20 2024-09-17 Know Labs, Inc. Non-invasive analyte sensing and notification system with decoupled and inefficient transmit and receive antennas
US12092589B1 (en) 2020-04-03 2024-09-17 Know Labs, Inc. In vitro analyte sensor using one or more detector arrays that operate in radio/microwave frequency bands
US12019034B2 (en) 2020-09-09 2024-06-25 Know Labs, Inc. In vitro sensing methods for analyzing in vitro flowing fluids
US11764488B2 (en) 2020-09-09 2023-09-19 Know Labs, Inc. Methods for determining variability of a state of a medium
WO2022056142A1 (en) 2020-09-09 2022-03-17 Apple Inc. Optical system for noise mitigation
US11510597B2 (en) 2020-09-09 2022-11-29 Know Labs, Inc. Non-invasive analyte sensor and automated response system
US12007338B2 (en) 2020-09-09 2024-06-11 Know Labs Inc. In vitro sensor for analyzing in vitro flowing fluids
US11689274B2 (en) 2020-09-09 2023-06-27 Know Labs, Inc. Systems for determining variability in a state of a medium
US20230355140A1 (en) 2022-05-05 2023-11-09 Know Labs, Inc. High performance glucose sensor
US11802843B1 (en) 2022-07-15 2023-10-31 Know Labs, Inc. Systems and methods for analyte sensing with reduced signal inaccuracy
US12033451B2 (en) 2022-08-15 2024-07-09 Know Labs, Inc. Systems and methods for analyte-based access controls
KR20240039440A (ko) * 2022-09-19 2024-03-26 부산대학교 산학협력단 근적외선 광원을 이용한 동맥경화반 성분분석이 가능한 광단층 내시경 시스템 및 이의 제어 방법
US11696698B1 (en) 2022-10-03 2023-07-11 Know Labs, Inc. Analyte sensors with position adjustable transmit and/or receive components
US11903701B1 (en) 2023-03-22 2024-02-20 Know Labs, Inc. Enhanced SPO2 measuring device

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005058154A1 (en) * 2003-12-16 2005-06-30 Medeikon Corporation Method for monitoring of analytes in biological samples using low coherence interferometry
US9464983B2 (en) * 2010-07-12 2016-10-11 Seiko Epson Corporation Concentration determination apparatus, probe, concentration determination method, and program
US11096608B2 (en) * 2013-12-03 2021-08-24 Imec Vzw Device and method for non-invasive measuring of analytes
US11026608B2 (en) * 2017-01-09 2021-06-08 Vox Biomedical Llc Cerebral oximetry using time-gated direct sequence spread spectrum
US12089928B2 (en) * 2017-05-18 2024-09-17 Wear2B Ltd. Device, system and method for non-invasive monitoring of physiological measurements
US11026604B2 (en) * 2017-07-13 2021-06-08 Cercacor Laboratories, Inc. Medical monitoring device for harmonizing physiological measurements
US10952616B2 (en) * 2018-03-30 2021-03-23 Canon U.S.A., Inc. Fluorescence imaging apparatus
KR20200047981A (ko) * 2018-10-29 2020-05-08 삼성전자주식회사 생체성분 추정 장치 및 방법
KR102143484B1 (ko) * 2018-11-12 2020-08-12 한국기초과학지원연구원 복합 현미경 시스템

Also Published As

Publication number Publication date
CN114295562A (zh) 2022-04-08
KR20220046168A (ko) 2022-04-14
EP3981329A1 (en) 2022-04-13
US20220104735A1 (en) 2022-04-07

Similar Documents

Publication Publication Date Title
JP2022061939A (ja) 分析物質の濃度推定装置及び方法、ならびに信号測定装置
KR102574088B1 (ko) 분석 물질의 농도 추정 장치 및 방법과, 농도 추정 모델 생성 장치 및 방법
EP3272280B1 (en) Apparatus and method for estimating blood sugar based on heterogeneous spectrums
CN110811635B (zh) 估计浓度的装置和方法以及生成浓度估计模型的装置
KR102693439B1 (ko) 스펙트럼 측정 장치, 스펙트럼의 광원 온도 변화 보정 방법, 농도 추정 장치 및 방법
KR102418120B1 (ko) 복수의 발광부 및 복수의 수광부를 포함하는 전자 장치
US11596331B2 (en) Apparatus and method for estimating analyte concentration
US10663405B2 (en) Raman probe and bio-component analyzing apparatus using the same
US11439327B2 (en) Apparatus and method for measuring Raman spectrum
US20210137445A1 (en) Apparatus and method for estimating skin barrier function
EP4184147B1 (en) Apparatus and method for estimating bio-information
KR20210022319A (ko) 생체정보 추정 장치 및 방법
US11883166B2 (en) Apparatus and method for estimating component of analyte
KR20210050967A (ko) 분석 물질의 농도 추정 장치 및 방법과, 캘리브레이션 방법
US20230119697A1 (en) Method of evaluating quality of bio-signal and apparatus for estimating bio-information
US20240225493A1 (en) Apparatus and method for estimating bio-information
US20220233149A1 (en) Apparatus and method for estimating body component
US20230309844A1 (en) Apparatus and method for estimating blood pressure
US11872032B2 (en) Apparatus and method for estimating analyte concentration
US20240172970A1 (en) Apparatus and method for estimating antioxidant component
KR20220025378A (ko) 대상 신호 스펙트럼 획득 장치 및 방법
KR20220029061A (ko) 대상 성분 추정 장치 및 방법
KR20200116751A (ko) 분석 물질의 농도 추정 장치 및 방법
KR20200133525A (ko) 생체 정보 추정 모델의 유효성 판단 장치 및 방법

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20231219