JP3931638B2 - 生体成分の定量装置 - Google Patents
生体成分の定量装置 Download PDFInfo
- Publication number
- JP3931638B2 JP3931638B2 JP2001350746A JP2001350746A JP3931638B2 JP 3931638 B2 JP3931638 B2 JP 3931638B2 JP 2001350746 A JP2001350746 A JP 2001350746A JP 2001350746 A JP2001350746 A JP 2001350746A JP 3931638 B2 JP3931638 B2 JP 3931638B2
- Authority
- JP
- Japan
- Prior art keywords
- skin
- light
- infrared spectrum
- biological component
- emitting
- 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.)
- Expired - Lifetime
Links
- 238000002329 infrared spectrum Methods 0.000 claims description 85
- 238000005259 measurement Methods 0.000 claims description 57
- 239000000523 sample Substances 0.000 claims description 47
- 238000011002 quantification Methods 0.000 claims description 26
- 238000000034 method Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000010521 absorption reaction Methods 0.000 claims description 7
- 238000002835 absorbance Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 5
- 238000004611 spectroscopical analysis Methods 0.000 claims description 5
- 238000012014 optical coherence tomography Methods 0.000 claims description 4
- 210000004003 subcutaneous fat Anatomy 0.000 claims description 2
- 239000000470 constituent Substances 0.000 claims 3
- 230000036548 skin texture Effects 0.000 claims 1
- 238000006467 substitution reaction Methods 0.000 claims 1
- 210000003491 skin Anatomy 0.000 description 132
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 89
- 239000008103 glucose Substances 0.000 description 89
- 239000008280 blood Substances 0.000 description 73
- 210000004369 blood Anatomy 0.000 description 73
- 210000001519 tissue Anatomy 0.000 description 28
- 210000004207 dermis Anatomy 0.000 description 24
- 239000013307 optical fiber Substances 0.000 description 13
- 239000000835 fiber Substances 0.000 description 12
- 210000002615 epidermis Anatomy 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 11
- 238000000513 principal component analysis Methods 0.000 description 11
- 206010033675 panniculitis Diseases 0.000 description 10
- 210000004304 subcutaneous tissue Anatomy 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 8
- 210000000245 forearm Anatomy 0.000 description 7
- 238000000491 multivariate analysis Methods 0.000 description 7
- 238000004445 quantitative analysis Methods 0.000 description 7
- 238000003325 tomography Methods 0.000 description 5
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 4
- 241000282414 Homo sapiens Species 0.000 description 4
- 210000000577 adipose tissue Anatomy 0.000 description 4
- 210000001736 capillary Anatomy 0.000 description 3
- 241000282412 Homo Species 0.000 description 2
- 238000000862 absorption spectrum Methods 0.000 description 2
- 238000011088 calibration curve Methods 0.000 description 2
- HVYWMOMLDIMFJA-DPAQBDIFSA-N cholesterol Chemical compound C1C=C2C[C@@H](O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2 HVYWMOMLDIMFJA-DPAQBDIFSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229910052736 halogen Inorganic materials 0.000 description 2
- 150000002367 halogens Chemical class 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000010238 partial least squares regression Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 210000000434 stratum corneum Anatomy 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 102000009027 Albumins Human genes 0.000 description 1
- 108010088751 Albumins Proteins 0.000 description 1
- 102000006395 Globulins Human genes 0.000 description 1
- 108010044091 Globulins Proteins 0.000 description 1
- 102000017011 Glycated Hemoglobin A Human genes 0.000 description 1
- 102000011782 Keratins Human genes 0.000 description 1
- 108010076876 Keratins Proteins 0.000 description 1
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- LEHOTFFKMJEONL-UHFFFAOYSA-N Uric Acid Chemical compound N1C(=O)NC(=O)C2=C1NC(=O)N2 LEHOTFFKMJEONL-UHFFFAOYSA-N 0.000 description 1
- TVWHNULVHGKJHS-UHFFFAOYSA-N Uric acid Natural products N1C(=O)NC(=O)C2NC(=O)NC21 TVWHNULVHGKJHS-UHFFFAOYSA-N 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 210000001124 body fluid Anatomy 0.000 description 1
- 239000010839 body fluid Substances 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 235000012000 cholesterol Nutrition 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- IXZISFNWUWKBOM-ARQDHWQXSA-N fructosamine Chemical compound NC[C@@]1(O)OC[C@@H](O)[C@@H](O)[C@@H]1O IXZISFNWUWKBOM-ARQDHWQXSA-N 0.000 description 1
- 150000002303 glucose derivatives Chemical class 0.000 description 1
- 108091005995 glycated hemoglobin Proteins 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 210000002751 lymph Anatomy 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 210000000944 nerve tissue Anatomy 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000036559 skin health Effects 0.000 description 1
- 230000036558 skin tension Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 229940116269 uric acid Drugs 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/1495—Calibrating or testing of in-vivo probes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/14532—Measuring 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 for measuring glucose, e.g. by tissue impedance measurement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/1455—Measuring 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0233—Special features of optical sensors or probes classified in A61B5/00
- A61B2562/0242—Special features of optical sensors or probes classified in A61B5/00 for varying or adjusting the optical path length in the tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1075—Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Optics & Photonics (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Emergency Medicine (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は生体成分を定量する装置に関するものであり、特に、近赤外領域における光の吸収を利用した分光分析装置に関するものである。さらには、生体組織中あるいは体液中の化学成分を分光分析する分光分析装置に関するものである。具体的な応用例としては、皮膚組織中のグルコース濃度の定量分析を行うことにより血糖値測定を行うための装置に応用されるものである。
【０００２】
【従来の技術】
分光分析装置の従来例として回折格子を用いた分光分析装置を示す。本従来例は生体組織中のグルコース濃度の定量分析を行うもので、本出願人が特願平１０−３０８８１７号「グルコース濃度の定量方法及びその装置」として出願したものであり、その概略図を図１に示す。
【０００３】
このグルコース濃度の定量装置は、ハロゲンランプからなる光源１と、光源１からの光を集光する集光レンズ２と、集光レンズ２を通過した光を被測定物１０に照射すると共に被測定物１０を透過あるいは拡散反射した光を受光するための受発光プローブ３と、受光後の光を分光する回折格子を収めた回折格子ユニット４、前記回折格子ユニット４で分光された光を検出するためのＩｎＧａＡｓアレイ型受光素子ユニット５、前記ＩｎＧａＡｓアレイ型受光素子ユニット５で得られた信号をもとにグルコース濃度を演算する演算ユニット６とから構成されている。
【０００４】
実際の動作を図１に基づいて説明すると、スペクトル測定に際して前記光源１からの光が前記受発光プローブ３を経由してセラミック板などの標準板を照射すると共に標準板を反射した光（参照信号）が受光される。続いて、接触位置を一定化するための位置決め治具を用い、前記受発光プローブ３の測定端部と皮膚表面の接触圧力を９．８〜４９ｋＰａ（１００〜５００ｇｆ／ｃｍ2）に設定して当接させ、同様にして、皮膚組織内を透過あるいは拡散反射した光（生体信号）が受光される。そして、得られた参照信号と生体信号をもとに前記演算ユニット６よりグルコース濃度が演算される。
【０００５】
グルコース濃度の定量には予め作成した検量式が用いられる。この検量式の作成には通常、多変量解析によるスペクトル解析手法が用いられ、標準手段で定量されたグルコース濃度を目的変量に、前記分光分析装置による測定で得た生体組織スペクトルを説明変数として多変量解析が行われる。多変量解析手段としては、重回帰分析、ＰＬＳ回帰分析、ニューラルネットワーク等の解析手段が用いられている。
【０００６】
皮膚に当接させ近赤外スペクトルを測定する受発光プローブ３の受発光間隔、すなわち、受発光プローブ３に設けた発光ファイバー７と受光ファイバー８の間隔Ｌは２ｍｍ以下に設定されている。人間を含む生物の皮膚組織は通常、図２に示すように角質層を含む表皮層１２、真皮層１３、皮下組織層１４の３層から形成されている。尚、表皮層１２の厚みは約０．２〜０．４ｍｍ、真皮層１３の厚みは約０．５〜２ｍｍ、皮下組織層１４の厚みは約１〜３ｍｍである。真皮層１３は毛細血管等が発達しており、血中グルコースに応じた物質移動が速やかに起こり、血中グルコース濃度（血糖値）に対して真皮層１３（真皮組織）中のグルコース濃度は追随して変化すると考えられている。皮下組織層１４は脂肪組織が中心であり、水溶性のグルコースは皮下組織層１４（皮下組織）中に均一に存在しにくい。
【０００７】
従って、血中グルコース濃度（血糖値）を精度よく測定するには、真皮層１３の近赤外スペクトルを選択的に測定する必要がある。我々は、真皮層１３の近赤外スペクトルを選択的に測定する手法として、２ｍｍ以下に受発光間隔を設定した受発光プローブ３を用いる方法を用いた。
【０００８】
この方法は図２に示すように、受発光プローブ３の発光ファイバー７から照射された近赤外光が皮膚組織中を透過、拡散し、後方散乱した光が受発光プローブ３の受光ファイバー８に到達する特性を利用したものである。受光ファイバー８に到達する近赤外光の経路Ａは、図２に示したＵ字型の経路Ａをとると考えられ、その経路の皮膚組織への到達深さは受発光間隔により変化する。従って、前記受発光プローブ３の受発光間隔を適切に設定することにより、相構造である皮膚組織の真皮層１３の近赤外スペクトルを選択的に測定することが可能となるのである。
【０００９】
【発明が解決しようとする課題】
従来の技術で示した真皮層１３（真皮組織）の近赤外スペクトルを選択的に測定する手法のみならず、公知になっている光学的に生体成分濃度を定量分析する手法は、個体差が大きい人間や生物の皮膚組織中の生体成分濃度を測定する場合に、通常、生体成分の定量を行う個体毎、あるいは測定部分毎に検量式を作成して定量分析が行われている。
【００１０】
しかしながら、個体毎、あるいは測定部分毎に検量式を作成する場合、測定対象である人や生物を長時間拘束し近赤外スペクトルおよび実測したデータを多変量解析することで検量式を作成する必要があり、測定対象の人や生物の負担が大きい。さらに、測定者にも検量のための式を作成する等の複雑な作業が必要となり、負担が大きい。
【００１１】
そこで、個体差が大きい人間や生物の皮膚組織における生体成分濃度の定量、あるいはその成分濃度を代用特性とした特性値の定量を、個体毎、あるいは測定部分毎に検量式を作成することなく定量する定量装置を提供することが本発明の課題である。
【００１２】
また、個体差を考慮する際に、その個体の属性を判断し最適な測定条件や検量式を選択することで測定精度を向上させることも本発明の課題である。
【００１３】
【課題を解決するための手段】
本発明の請求項１に係る生体成分の定量装置は、皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニット６により、目的とする生体成分濃度を定量する定量装置において、前記演算ユニット６により、皮膚厚さを指標として皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有し、前記皮膚厚さを指標として皮膚性状の分類を行うにあたって、皮膚の近赤外スペクトルの形状を実際の皮膚厚さの代用特性とすることを特徴とするものである。
【００１９】
また、本発明の請求項２に係る生体成分の定量装置は、請求項１に加えて、前記皮膚の近赤外スペクトルの形状を実際の皮膚厚さの代用特性として行う手法として、皮下脂肪の特異吸収波長である１７００〜１８００ｎｍの波長範囲から少なくとも一つの波長を選択し、その吸光度から皮膚性状の分類を行うことを特徴とするものである。
【００２１】
また、本発明の請求項３に係る生体成分の定量装置は、皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニット６により、目的とする生体成分濃度を定量する定量装置において、前記演算ユニット６により、皮膚厚さを指標として皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有し、前記皮膚の皮膚厚さを指標として皮膚性状の分類を行うにあたって、光コヒーレンストモグラフィを利用することを特徴とするものである。
【００２２】
本発明の請求項４に係る生体成分の定量装置は、皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニット６により、目的とする生体成分濃度を定量する定量装置において、前記演算ユニット６により、皮膚の近赤外スペクトルで皮膚性状を分類する段階、分類した皮膚性状に対応する近赤外スペクトルの測定条件を設定する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有することを特徴とするものである。
【００２４】
本発明の請求項５に係る生体成分の定量装置は、受発光間隔を２ｍｍ以下に設定した受発光プローブ３と、前記受発光プローブ３に近赤外光を供給する光源１と、前記受発光プローブ３で受光された近赤外光を分光する分光手段と、前記分光手段で分光されたスペクトル光を受光する受光手段と、前記受光手段からの信号を用い生体成分の定量を行う演算手段、皮膚性状に対応する検量式を複数記憶させる記憶手段からなる生体成分の定量装置であり、生体成分を定量する際に、少なくとも、前記受発光プローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、前記受発光プローブを当接させた部分の皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階を経て目的とする生体成分濃度を定量することを特徴とするものである。
【００２５】
本発明の請求項６に係る生体成分の定量装置は、受発光間隔を２ｍｍ以下に設定した受発光プロ−ブ３と、前記受発光間隔を選択できる受発光間隔選択手段と、前記受発光プローブに近赤外光を供給する光源１と、前記受発光プローブ３で受光された近赤外光を分光する分光手段と、前記分光手段で分光されたスペクトル光を受光する受光手段と、前記受光手段からの信号を用い生体成分の定量を行う演算手段、皮膚性状に対応する検量式を複数記憶させる記憶手段からなる生体成分の定量装置であり、生体成分を定量する際に、少なくとも、前記受発光プローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、前記受発光プローブを当接させた部分の皮膚性状を分類する段階、分類した皮膚性状に対応する前記プローブの受発光間隔を選択する段階、分類した皮膚性状と受発光間隔に対応する検量式を選択する段階、設定した受発光間隔のプローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階を経て目的とする生体成分濃度を定量することを特徴とするものである。
【００２６】
【発明の実施の形態】
以下に本発明の具体的な実施の形態を示す。
【００２７】
本発明の生体成分の定量装置は、外観上、図１に示す従来のものと同様に形成されるものであって、ハロゲンランプ等からなる光源１と、光源１からの光を集光するための集光レンズ２と、集光レンズ２を通過した光を被測定物１０に照射すると共に被測定物１０を透過あるいは拡散反射した光を受光するための受発光プローブ３と、受光後の光を分光する回折格子を収めた回折格子ユニット４、前記回折格子ユニット４で分光された光（スペクトル光）を検出するためのＩｎＧａＡｓアレイ型受光素子ユニット５、前記ＩｎＧａＡｓアレイ型受光素子ユニット５で得られた信号をもとにグルコース濃度を演算するノート型パーソナルコンピュータ等の演算ユニット６とから構成されている。すなわち、上記の回折格子ユニット４が受発光プローブ３で受光された近赤外光を分光する分光手段として、上記のＩｎＧａＡｓアレイ型受光素子ユニット５が分光手段で分光されたスペクトル光を受光する受光手段として、上記の演算ユニット６が受光手段からの信号を用い生体成分の定量を行う演算手段及び皮膚性状に対応する検量式を複数記憶させる記憶手段として、それぞれ形成されている。そして、上記の演算ユニット６を使用することにより、皮膚性状を判断して分類する段階、分類した皮膚性状に対応する検量式を選択する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階がそれぞれ行なわれるものである。
【００２８】
また、皮膚に当接させ近赤外スペクトルを測定する受発光プローブ３は上記と同様に形成されるものであって、受発光プローブ３の受発光間隔Ｌ、すなわち、受発光プローブ３に設けた発光ファイバー７と受光ファイバー８の間隔Ｌは２ｍｍ以下に設定されている。
【００２９】
また、図１５に示す他の受発光プローブ３は、中心の受光ファイバー８の同心円状に受発光間距離（受発光間隔Ｌ１）３００μｍに複数本（６本）の発光ファイバー７ａと受発光間距離（受発光間隔Ｌ２）６５０μｍに複数本（１２本）の発光ファイバー７ｂを配置したものを用いた。受発光間距離３００μｍと６５０μｍからの発光は切り替え手段により切り替えられるようになっている。この切り替え手段が受発光間隔を選択できる受発光間隔選択手段として形成されるものであって、具体的には、光源１からの光をコリメートレンズにより平行光とした後、ハーフミラーにより３００μｍと６５０μｍの光ファイバー７ａ、７ｂへ１：２の割合で光量が分割されるように調整する。前記ハーフミラーで分割された光はそれぞれ集光ミラーや集光レンズで光ファイバー７ａ、７ｂの測定端と反対の端部に集光される。また、受発光間隔Ｌ１、Ｌ２の選択はこのハーフミラーと集光ミラーや集光レンズの間に配置した２つのメカニカルシャッターの開閉で行う。すなわち、３００μｍの受発光間隔Ｌ１を選択する時は３００μｍ側のメカニカルシャッターを開き、６５０μｍ側を閉じる。また、６５０μｍの受発光間隔Ｌ２を選択する時は６５０μｍ側のメカニカルシャッターを開き、３００μｍ側を閉じる。尚、この受発光プローブ３においても上記の受発光間距離（受発光間隔Ｌ１、Ｌ２）は３００μｍや６５０μｍに限定されるものではなく、上記と同様に２ｍｍ以下で皮膚性状（皮膚厚み）に応じて適宜に設定することができる。また、上記の受発光間隔Ｌは０．２ｍｍ以上に設定するのが好ましい。
【００３０】
一人の被験者（Ｐ１）から得られた皮膚（皮膚組織）の近赤外スペクトルを示すグラフを図３に、検量式の回帰係数を示すグラフを図４に示す。図３の近赤外スペクトルは図１の装置を用いて測定した。図４に示す検量式の回帰係数は経口でグルコース負荷を行い人為的に血糖値を変動させる実験を６回行い、実験で得られた皮膚の近赤外スペクトルと血液から測定した血糖値からなるデータを多変量解析することにより求めた。多変量解析はＰＬＳ回帰分析手法を用いている。また、解析に用いた近赤外スペクトルは微分処理あるいはＭＳＣ処理等の前処理をほどこさない生データを用いたが、これに限るものではない。被験者の性別は男、年齢は４１才、健常人である。皮膚の近赤外スペクトルは左手前腕内側部で行い、測定を行った部分の皮膚厚さは表皮層１２と真皮層１３の厚さの合計は約１．５ｍｍであった。皮膚厚さの測定はＣＯＲＴＥＸ社製（デンマーク）の超音波断層測定装置ＤｅｒｍａＳｃａｎ Ｃ Ｖｅｒ．３を用いて行った。
グルコース濃度、すなわち血糖値測定に用いられる検量式は
ａｎ：波長ｎにおける回帰係数
ｘｎ：波長ｎにおける吸光度
ｂ：定数
で表される。したがって、この検量式の回帰係数によりグルコース濃度の定量に重要な波長が推定できる。図４に検量線の回帰係数のグラフを示す。図４の回帰係数は１６００ｎｍ付近に特徴的な正のピークを有することがわかる。代表的な生体成分の吸光スペクトルを図５に示す。図５に示されるようにグルコース分子が１６００ｎｍ付近に他の生体成分に比較して特異な吸収を有することがわかる。このことから、上記実験で得られた検量式の回帰係数が生体皮膚組織中のグルコース分子の特異吸収を利用してグルコース濃度を定量していることがわかる。
【００３１】
上記の生体皮膚組織中グルコース濃度（血糖値）測定は上記従来例とほぼ同様であって、図６に示したような中心に配置した受光ファイバー８と、その受光ファイバー８を中心とする半径６５０μｍの円周上に等間隔で配置した１２本の発光ファイバー７からなる受発光プローブ３を皮膚面に９．８〜４９ｋＰａ（１００〜５００ｇｆ／ｃｍ2）、好ましくは２９．４ｋＰａ（３００ｇｆ／ｃｍ2）の一定加重で当接させ、近赤外スペクトルを測定することで行う。光源１からの光は円周上に配置された１２本の発光ファイバー７より照射され、皮膚組織を透過した光の一部が中心の受光ファイバー８で受光される。受光ファイバー８で受光される光の経路はＵ字型の経路Ａをとっていると推定でき、皮膚組織への到達深さは受光ファイバー８と発光ファイバー７間の距離で決まる。
【００３２】
人間の皮膚組織構造は、角質層を含む表皮層（表皮組織）１２、真皮層（真皮組織）１３、皮下組織層（皮下組織）１４の３層で構成される。真皮層１３には毛細血管、リンパ、神経組織が発達している。皮下組織層１４は主に脂肪組織で構成される。グルコースは水溶性のため脂肪組織中には存在しにくい。また、表皮層１２には毛細血管が発達していないため、血液で輸送されるグルコースが到達しにくい。そのため、血糖値の測定を皮膚の近赤外スペクトル測定を行うためには、真皮層１２の近赤外スペクトルを選択的に測定する必要がある。
【００３３】
実験に用いた受発光プローブ３は、受発光間隔Ｌを６５０μｍに設定したもので、表皮層１２と真皮層１３をあわせた皮膚厚さが１．５ｍｍである前記被験者から作成した検量式である。ある個体に対して作成された検量式は、その個体に適切であっても、皮膚厚さ等の皮膚性状に個体差があるため、そのまま他の個体に適用できるとは限らない。従来の他研究において個人毎の検量線が使われたのはそのためである。しかしながら、個体差があっても定量を目的とする生体成分（上記実験ではグルコース）は個体によらず共通であることから、その生体成分の変動（上記実験ではグルコースの変動）を捉えた検量式を用いれば、ある程度一般化された検量式を作成できることが期待できる。本発明に用いた検量式の一つは図４に示すような回帰係数形状を有する。この回帰係数を持つ検量式は、定量を目的とするグルコース分子の特異吸収波長である１６００ｎｍ付近の吸収を正のピークとして利用していることから、皮膚組織内のグルコース分子の変動を定量していることが期待できる。
【００３４】
次に、この検量式を用いて被験者（Ｐ１）の自身の血糖値変動を測定した結果を図７に示す。
【００３５】
同じ検量式を用いて他の４人の被験者（Ｐ２〜Ｐ５）の血糖値変動を測定する実験を合わせて実施した。被験者（Ｐ１）で作成した検量式を用いて、他の４人の被験者（Ｐ２〜Ｐ５）の血糖値変動を推定した場合の結果を図８に示す。また、この実験を行った被験者を含む被験者５名の属性を示す表を表１に示す。
【００３６】
【表１】
【００３７】
図７のグラフから、従来の技術と同様に、被験者（Ｐ１）から作成した検量式を用いて、同じ被験者（Ｐ１）の血糖値変動を推定する場合は相関係数０．９程度の良好な推定が可能であることがわかる。
【００３８】
次に、図８からわかるように、被験者（Ｐ１）で作成した検量式を用いて他の被験者の血糖値変動を推定した場合、被験者Ｐ２の血糖値変動については良好な測定が可能であるが、他の被験者Ｐ３、Ｐ４、Ｐ５についての血糖値測定は困難であることがわかった。図９に被験者（Ｐ１）で作成した検量式を用いて５名の被験者（Ｐ１〜Ｐ５）の血糖値変動を推定した結果のグラフを示す。この場合の相関係数は０．６９であった。
【００３９】
次に、被験者Ｐ２〜Ｐ５の被験者についても同様な６回のグルコース負荷実験を実施し、被験者毎に検量式を作成した。被験者（Ｐ１）で作成した検量式を含む６個の検量式を用いて、５名の被験者相互に血糖値変動の推定を行った。この結果を表２に示す。
【００４０】
【表２】
【００４１】
表２は作成した検量式に他の被験者毎のグルコース負荷実験で得られた近赤外スペクトルデータを代入し、血糖値変動を推定し、実測した血糖値と比較した結果をまとめたものである。この比較において相関係数（６回の当てはめの相関係数の平均値）が０．７以上の場合を○印、０．６以上０．７未満を△印、０．６未満を×印で示している。表２から分かるように、５名の被験者を用いたこの実験では、大きく２つの属性（Ａ群およびＢ群）に分かれ、その群内においてお互いに作成した検量式で血糖値の変動が測定できることが判明した。
【００４２】
表２と表１の結果を比較すると２群は皮膚厚さで層別できると考えられ、被験者の属性を分類する指標の一つとして皮膚厚さがあることが分かる。表１の皮膚厚さを求めた超音波断層写真を図１０に示す。皮膚厚さの値は、この断層写真より任意に選んだ３箇所の皮膚厚さ（表皮層と真皮層を合わせたもの）の平均値で算出した。
【００４３】
近赤外分光分析において精度よい検量式を作成するには、通常、多くのデータ蓄積をする必要がある。したがって、上記の血糖値測定においても検量式の作成には数日におよび被験者を拘束しデータ収集を行う必要があった。非侵襲的に、かつ、迅速に血糖値が測定できる近赤外線による定量手法の恩恵ははかり知れないが、そのための検量式作成に利用者の数日におよぶ拘束が必要とするのであれば、実用上の大きな問題と言わざるをえない。上記実験結果に示すように、我々は、利用者の個人毎で検量式を作成するのでなく、生体成分の測定対象である患者や被験者をある特性毎に分類し、分類した群に対応する検量式を前もって用意することで煩雑な検量式の作成プロセスを必要としない生体成分の定量方法および装置を提供できることがわかった。また、その分類は、皮膚厚さが大きな要因であることがわかった。また、予め複数の検量式を用意し、分類した群毎に適切な検量式を用いることで測定精度の向上がはかれることがわかった。
【００４４】
被験者の属性を分類する指標としては皮膚厚さ以外に、皮膚の角質水分量、皮膚組織の水分濃度、皮膚組織の密度、皮膚の色、皮膚の表面粗さ、性別、年齢、人種を使用することができる。
【００４５】
上記説明では定量分析を行う生体成分として主にグルコースあるいは血糖値を記述しているが、本発明はこれに限らず、組織水分量、中性脂肪、コレステロール、ＨｂＡ１ｃ（糖化ヘモグロビン）、フルクトサミン、アルブミン、グロブリン、尿酸等の生体成分の定量分析に利用できる。また、これらを代用特性とする皮膚の健康度や皮膚年齢、老化度、皮膚の張力（張り）等の測定に利用できる。
【００４６】
【実施例】
−参考例−
本参考例は超音波断層撮影装置を利用し、被験者の皮膚厚さを分類した例である。皮膚厚さの測定はＣＯＲＴＥＸ社製（デンマーク）の超音波断層測定装置ＤｅｒｍａＳｃａｎＣ Ｖｅｒ．３を用いて行った。皮膚厚さの測定結果を図１０に示す。
【００４７】
本参考例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。近赤外スペクトルから血糖値を推定するための演算手段では、表皮層と真皮層との合計の厚さが１．２ｍｍ以上の被験者の血糖値を推定する検量式と、１．２ｍｍ未満の被験者の血糖値を推定する検量式の２種類の検量式が準備されている。
【００４８】
表皮層と真皮層との合計の厚さが１．２ｍｍ以上の被験者の血糖値を推定する検量式は、
である。
【００４９】
また、表皮層と真皮層との合計の厚さが１．２ｍｍ未満の被験者の血糖値を推定する検量式は、
である。
【００５０】
尚、ａｎとｃｎは波長ｎにおける回帰係数、ｘｎは波長ｎにおける吸光度、ｂとｄは定数をそれぞれ示す。
【００５１】
皮膚の近赤外スペクトルによる血糖値の推定と、この測定値の測定精度を判断するための実測値の測定は以下の手順で行った。
（１）超音波断層撮影装置を利用し、被験者の皮膚厚さを測定する。
【００５２】
皮膚厚さの測定は、近赤外スペクトルを測定する同一部位である左手の前腕内側で行った。
（２）測定した皮膚厚さに基づいて、被験者の分類を行う。
（３）被験者の分類に応じた検量式を選択する。
（４）血糖値の実測を行う。
【００５３】
この段階は実際の測定では必要ないが、本参考例での測定精度を検証するために実施した。
（５）皮膚の近赤外スペクトルを測定する。
（６）測定した近赤外スペクトルを選択した検量式に代入し、血糖値を推定する。
【００５４】
上記手順は被験者が一人に固定される場合、２回目以降は皮膚厚さを測定する段階は省略し、分類された群の測定条件、検量式をそのまま用いてもよい。
【００５５】
本参考例での血糖定量の結果として、皮膚厚さ１．２ｍｍ以上と判断され、血糖値を推定された群の測定結果を図１１に、皮膚厚さ１．２ｍｍ未満と推定され、血糖値を推定された群の測定結果を図１２に示す。また、対照例として、同一データを被験者を分類することなく、一律に皮膚厚さ１．２ｍｍの血糖値推定に用いた検量式で血糖値推定を行った場合の結果を図９に示す。被験者を分類せず血糖値推定を行った場合の相関係数は０．６９で、被験者の分類を行った場合の相関係数は０．９３と０．８５であることから、本参考例による定量精度の向上が確認できた。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
【００５６】
本参考例では、皮膚厚さを実測する手段として超音波断層装置を用いたが、本発明はこれに限定するものではなく、断層像として計測しない超音波測定手法を用いても構わない。また、本発明は光コヒーレンストモグラフィを用いて皮膚厚さの実測を行うようにする。
【００５７】
超音波による皮膚厚さの測定は画像を用いて行う手法と、画像化せずに超音波の反射により距離を直接求める手法がある。従って、上記参考例のように断層画像を用いなくても超音波の反射信号を用いて皮膚厚さを求めることができる。
【００５８】
また、光コヒーレンストモグラフィーの代表的なものとして、低コヒーレンス干渉法と光周波数走査法がある。いずれの手法も数十μｍオーダーの分解能で生体の断層測定を可能としている。そして、近赤外レーザー光を前腕内側部に照射し、得られた皮膚断層画像から皮膚厚さを測定することができるものである。
−実施例１−
本実施例は測定した近赤外スペクトルより被験者の分類を行った例である。本実施例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。近赤外スペクトルの測定は、左手の前腕内側で行った。また、測定に用いた受発光プローブは従来の技術に示したものと同じ受発光ファイバー間が６５０μｍのものを用いた。
【００５９】
近赤外スペクトルを用いて行う被験者の分類は、主成分分析により行った。図１３に主成分分析により被験者の分類を行った結果を示す。図１３は、表皮層と真皮層との合計の厚さが１．２ｍｍ以上の被験者と、１．２ｍｍ未満の被験者が第２主成分（ＰＣ２）により明確に分類できることを示している（図１３の上部と下部）。
【００６０】
図１４に第２主成分のローディングプロットを示す。第２主成分のローディングプロットは波長１７２８ｎｍと１７６９ｎｍに特徴的なピークを有し、このピークは図５に示した中性脂肪のピークと一致する。したがって、第２主成分の寄与の大きさは皮下組織（脂肪組織）への光の到達度と解釈でき、皮膚厚さの代用特性として主成分分析の分類が利用できることを示している。すなわち、皮膚厚さの異なる個体差は皮下組織への光の到達量の違いとなり、最終的に近赤外スペクトルの形状差として表れる。この近赤外スペクトルの形状差を主成分分析等の多変量解析を用いて分類することで、被験者の分類が可能である。多変量解析による分類は主成分分析に限らず、クラスタ分析、判別分析等の定性分析手法を利用することができる。
【００６１】
血糖値を推定する検量式は、参考例と同様に２種類準備されている。皮膚の近赤外スペクトルによる血糖値の推定と、この測定値の測定精度を判断するための実測値の測定は以下の手順で行った。
（１）近赤外スペクトルの測定を行う。
（２）測定した近赤外スペクトルの主成分分析を行い、被験者の分類を行う。
（３）被験者の分類に応じた検量式を選択する。
（４）測定した近赤外スペクトルを選択した検量式に代入し、血糖値を推定する。
（５）血糖値の実測を行う。
【００６２】
この段階は実際の測定では必要ないが、本実施例での測定精度を検証するために実施した。
【００６３】
上記手順は被験者が一人に固定される場合、２回目以降は主成分分析する段階は省略し、分類された群の測定条件、検量式をそのまま用いてもよい。
【００６４】
本実施例での血糖定量の結果として、皮膚厚さ１．２ｍｍ以上と判断され、血糖値を推定された群の測定結果を図１１に、皮膚厚さ１．２ｍｍ未満と推定され、血糖値を推定された群の測定結果を図１２に示す。また、対照例として、同一データを被験者を分類することなく、一律に皮膚厚さ１．２ｍｍの血糖値推定に用いた検量式で血糖値推定を行った場合の結果を図９に示す。被験者を分離せず血糖値推定を行った場合の相関係数は０．６９で、被験者の分類を行った場合の相関係数は０．８５と０．９３であることから、本発明による定量精度の向上が確認できた。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
−実施例２−
本実施例は予備的に測定した近赤外スペクトルより被験者の分類を行った例である。
【００６５】
本実施例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。予備的な近赤外スペクトルの測定は、血糖値を推定するために近赤外スペクトルを測定する同一部位である左手の前腕内側で行った。また、測定に用いた受発光プローブは図１５に示すように、中心の受光ファイバー８の同心円状に受発光間距離３００μｍに６本の発光ファイバー７ａと受発光間距離６５０μｍに１２本の発光ファイバー７ｂを配置したものを用いた。受発光間距離３００μｍと６５０μｍからの発光は切り替え手段により切り替えられるようになっており、予備的な近赤外スペクトルの測定は受発光間距離６５０μｍを用いて行った。
【００６６】
予備的な近赤外スペクトルを用いて行う被験者の分類は、実施例１と同様に主成分分析により行った。本実施例では、主成分分析によって表皮層と真皮層との合計の厚さが１．２ｍｍ未満と判断した場合は、発光の切り替え手段によって受発光間隔３００μｍの近赤外スペクトルを測定する。これは、皮膚厚さが薄い群については真皮層を適切に透過する光路を設定した方がより血糖値の推定精度が向上するためで、受発光間距離を狭めることで皮膚の深さ方向への到達度を浅くする効果がある。
【００６７】
血糖値を推定する検量式は、他の実施例と同様に２種類準備されている。皮膚の近赤外スペクトルによる血糖値の推定と、この測定値の測定精度を判断するための実測値の測定は以下の手順で行った。
▲１▼ 予備的な皮膚の近赤外スペクトルの測定を受発光間隔６５０μｍで行う。
▲２▼ 測定した皮膚厚さに基づいて、被験者の分類を行う。
▲３▼ 被験者の分類に応じた検量式を選択する。
▲４▼ 血糖値の実測を行う。
【００６８】
この段階は実際の測定では必要ないが、本実施例での測定精度を検証するために実施した。
▲５▼ 皮膚厚さを１．２ｍｍ未満と判断した場合、受発光間距離３００μｍで皮膚の近赤外スペクトルを測定する。皮膚厚さを１．２ｍｍ以上と判断した場合、▲１▼で測定した皮膚の近赤外スペクトルをそのまま利用する。
▲６▼ 測定した近赤外スペクトルを選択した検量式に代入し、血糖値を推定する。
上記手順は被験者が一人に固定される場合、２回目以降は予備的な皮膚の近赤外スペクトルを測定する段階は省略し、分類された群の測定条件、検量式をそのまま用いてもよい。
【００６９】
本発明により従来の技術のような単一の測定条件でなく、個体差に応じて受発光ファイバー間隔を設定し、その分類に応じた検量式を準備することで定量精度の向上が可能である。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
−実施例３−
本実施例は測定した近赤外スペクトルより被験者の分類を行った例である。本実施例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。近赤外スペクトルの測定は、左手の前腕内側で行った。また、測定に用いた受発光プローブは従来の技術に示したものと同じ受発光ファイバー間が６５０μｍのものを用いた。
【００７０】
近赤外スペクトルを用いて行う被験者の分類は、近赤外スペクトルの１７２８ｎｍでの吸光度（Ａｂ１７２８）と１６７０ｎｍでの吸光度（Ａｂ１６７０）の比を用いて行った。すなわち、
Ａｂ１７２８／Ａｂ１６７０＞αの時は皮膚厚さ１．２ｍｍ未満
Ａｂ１７２８／Ａｂ１６７０＜α，Ａｂ１７２８／Ａｂ１６７０＝αの時は皮膚厚さ１．２ｍｍ以上
ここでαは定数
と判断し、皮膚性状の分類を行い、他の実施例と同様に、事前に準備した皮膚性状に対応する検量式に近赤外スペクトルを代入して血糖値の算出を行う。
【００７１】
本発明により、近赤外スペクトルによる皮膚性状の分類を脂肪成分由来の特異吸収波長を利用して行その分類に応じた検量式を準備することで定量精度の向上が可能である。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
【００７２】
【発明の効果】
従来の技術では、個体差が大きい人間や生物の皮膚組織中の生体成分濃度を測定する場合に、通常、生体成分の定量を行う個体毎、あるいは測定部分毎に検量式を作成し定量分析を行うことが行われるが、本発明のように、被験者を皮膚性状に分類し、その分類毎に検量式を準備することで、個体差が大きい人間や生物の皮膚組織における生体成分濃度の定量、あるいはその成分濃度を代用特性とした特性値の定量を、個体毎、あるいは測定部分毎に検量式を作成することなくできるようになった。それにより、個体毎あるいは測定部分毎に検量線を作成する場合での、測定対象である人や生物を長時間拘束する必要がなくなり、測定対象の人や生物の負担を軽減することができた。さらに、測定者にも検量のための式を作成する等の複雑な作業の必要がなくなり、負担を軽減することができた。
【００７３】
また、個体差を考慮する際に、その個体の属性を判断し最適な測定条件や検量式を選択することで測定精度を向上させることもができた。
【図面の簡単な説明】
【図１】 本発明及び従来例を示す概略図である。
【図２】 同上の皮膚組織の断面模式図及び真皮のスペクトル測定を示す概念図である。
【図３】 同上の前腕内側部で測定した皮膚の近赤外スペクトルを示すグラフである。
【図４】 同上の血糖値定量を行う検量式の回帰係数を示すグラフである。
【図５】 同上の生体成分の吸収スペクトルを示すグラフである。
【図６】 同上の受発光プローブを示す正面図である。
【図７】 従来技術による血糖値測定結果を示すグラフである。
【図８】 本発明の一つの検量式での複数の被験者の血糖値測定を行った結果を示すグラフである。
【図９】 同上の一つの検量式での複数の被験者の血糖値測定を行った結果を示すグラフである。
【図１０】 同上の被験者の皮膚厚さを求めた超音波断層写真である。
【図１１】 同上の皮膚厚さ１．２ｍｍ以上の群に対して血糖値変動を求めた結果を示すグラフである。
【図１２】 同上の皮膚厚さ１．２ｍｍ未満の群に対して血糖値変動を求めた結果を示すグラフである。
【図１３】 同上の皮膚の近赤外スペクトルを主成分分析した結果を示すグラフである。
【図１４】 同上の主成分分析の第２主成分のローディングプロットである。
【図１５】 同上の実施例２に用いた受発光プローブを示す正面図である。
【符号の説明】
１ 光源
３ 受発光プローブ
６ 演算ユニット
【発明の属する技術分野】
本発明は生体成分を定量する装置に関するものであり、特に、近赤外領域における光の吸収を利用した分光分析装置に関するものである。さらには、生体組織中あるいは体液中の化学成分を分光分析する分光分析装置に関するものである。具体的な応用例としては、皮膚組織中のグルコース濃度の定量分析を行うことにより血糖値測定を行うための装置に応用されるものである。
【０００２】
【従来の技術】
分光分析装置の従来例として回折格子を用いた分光分析装置を示す。本従来例は生体組織中のグルコース濃度の定量分析を行うもので、本出願人が特願平１０−３０８８１７号「グルコース濃度の定量方法及びその装置」として出願したものであり、その概略図を図１に示す。
【０００３】
このグルコース濃度の定量装置は、ハロゲンランプからなる光源１と、光源１からの光を集光する集光レンズ２と、集光レンズ２を通過した光を被測定物１０に照射すると共に被測定物１０を透過あるいは拡散反射した光を受光するための受発光プローブ３と、受光後の光を分光する回折格子を収めた回折格子ユニット４、前記回折格子ユニット４で分光された光を検出するためのＩｎＧａＡｓアレイ型受光素子ユニット５、前記ＩｎＧａＡｓアレイ型受光素子ユニット５で得られた信号をもとにグルコース濃度を演算する演算ユニット６とから構成されている。
【０００４】
実際の動作を図１に基づいて説明すると、スペクトル測定に際して前記光源１からの光が前記受発光プローブ３を経由してセラミック板などの標準板を照射すると共に標準板を反射した光（参照信号）が受光される。続いて、接触位置を一定化するための位置決め治具を用い、前記受発光プローブ３の測定端部と皮膚表面の接触圧力を９．８〜４９ｋＰａ（１００〜５００ｇｆ／ｃｍ2）に設定して当接させ、同様にして、皮膚組織内を透過あるいは拡散反射した光（生体信号）が受光される。そして、得られた参照信号と生体信号をもとに前記演算ユニット６よりグルコース濃度が演算される。
【０００５】
グルコース濃度の定量には予め作成した検量式が用いられる。この検量式の作成には通常、多変量解析によるスペクトル解析手法が用いられ、標準手段で定量されたグルコース濃度を目的変量に、前記分光分析装置による測定で得た生体組織スペクトルを説明変数として多変量解析が行われる。多変量解析手段としては、重回帰分析、ＰＬＳ回帰分析、ニューラルネットワーク等の解析手段が用いられている。
【０００６】
皮膚に当接させ近赤外スペクトルを測定する受発光プローブ３の受発光間隔、すなわち、受発光プローブ３に設けた発光ファイバー７と受光ファイバー８の間隔Ｌは２ｍｍ以下に設定されている。人間を含む生物の皮膚組織は通常、図２に示すように角質層を含む表皮層１２、真皮層１３、皮下組織層１４の３層から形成されている。尚、表皮層１２の厚みは約０．２〜０．４ｍｍ、真皮層１３の厚みは約０．５〜２ｍｍ、皮下組織層１４の厚みは約１〜３ｍｍである。真皮層１３は毛細血管等が発達しており、血中グルコースに応じた物質移動が速やかに起こり、血中グルコース濃度（血糖値）に対して真皮層１３（真皮組織）中のグルコース濃度は追随して変化すると考えられている。皮下組織層１４は脂肪組織が中心であり、水溶性のグルコースは皮下組織層１４（皮下組織）中に均一に存在しにくい。
【０００７】
従って、血中グルコース濃度（血糖値）を精度よく測定するには、真皮層１３の近赤外スペクトルを選択的に測定する必要がある。我々は、真皮層１３の近赤外スペクトルを選択的に測定する手法として、２ｍｍ以下に受発光間隔を設定した受発光プローブ３を用いる方法を用いた。
【０００８】
この方法は図２に示すように、受発光プローブ３の発光ファイバー７から照射された近赤外光が皮膚組織中を透過、拡散し、後方散乱した光が受発光プローブ３の受光ファイバー８に到達する特性を利用したものである。受光ファイバー８に到達する近赤外光の経路Ａは、図２に示したＵ字型の経路Ａをとると考えられ、その経路の皮膚組織への到達深さは受発光間隔により変化する。従って、前記受発光プローブ３の受発光間隔を適切に設定することにより、相構造である皮膚組織の真皮層１３の近赤外スペクトルを選択的に測定することが可能となるのである。
【０００９】
【発明が解決しようとする課題】
従来の技術で示した真皮層１３（真皮組織）の近赤外スペクトルを選択的に測定する手法のみならず、公知になっている光学的に生体成分濃度を定量分析する手法は、個体差が大きい人間や生物の皮膚組織中の生体成分濃度を測定する場合に、通常、生体成分の定量を行う個体毎、あるいは測定部分毎に検量式を作成して定量分析が行われている。
【００１０】
しかしながら、個体毎、あるいは測定部分毎に検量式を作成する場合、測定対象である人や生物を長時間拘束し近赤外スペクトルおよび実測したデータを多変量解析することで検量式を作成する必要があり、測定対象の人や生物の負担が大きい。さらに、測定者にも検量のための式を作成する等の複雑な作業が必要となり、負担が大きい。
【００１１】
そこで、個体差が大きい人間や生物の皮膚組織における生体成分濃度の定量、あるいはその成分濃度を代用特性とした特性値の定量を、個体毎、あるいは測定部分毎に検量式を作成することなく定量する定量装置を提供することが本発明の課題である。
【００１２】
また、個体差を考慮する際に、その個体の属性を判断し最適な測定条件や検量式を選択することで測定精度を向上させることも本発明の課題である。
【００１３】
【課題を解決するための手段】
本発明の請求項１に係る生体成分の定量装置は、皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニット６により、目的とする生体成分濃度を定量する定量装置において、前記演算ユニット６により、皮膚厚さを指標として皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有し、前記皮膚厚さを指標として皮膚性状の分類を行うにあたって、皮膚の近赤外スペクトルの形状を実際の皮膚厚さの代用特性とすることを特徴とするものである。
【００１９】
また、本発明の請求項２に係る生体成分の定量装置は、請求項１に加えて、前記皮膚の近赤外スペクトルの形状を実際の皮膚厚さの代用特性として行う手法として、皮下脂肪の特異吸収波長である１７００〜１８００ｎｍの波長範囲から少なくとも一つの波長を選択し、その吸光度から皮膚性状の分類を行うことを特徴とするものである。
【００２１】
また、本発明の請求項３に係る生体成分の定量装置は、皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニット６により、目的とする生体成分濃度を定量する定量装置において、前記演算ユニット６により、皮膚厚さを指標として皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有し、前記皮膚の皮膚厚さを指標として皮膚性状の分類を行うにあたって、光コヒーレンストモグラフィを利用することを特徴とするものである。
【００２２】
本発明の請求項４に係る生体成分の定量装置は、皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニット６により、目的とする生体成分濃度を定量する定量装置において、前記演算ユニット６により、皮膚の近赤外スペクトルで皮膚性状を分類する段階、分類した皮膚性状に対応する近赤外スペクトルの測定条件を設定する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有することを特徴とするものである。
【００２４】
本発明の請求項５に係る生体成分の定量装置は、受発光間隔を２ｍｍ以下に設定した受発光プローブ３と、前記受発光プローブ３に近赤外光を供給する光源１と、前記受発光プローブ３で受光された近赤外光を分光する分光手段と、前記分光手段で分光されたスペクトル光を受光する受光手段と、前記受光手段からの信号を用い生体成分の定量を行う演算手段、皮膚性状に対応する検量式を複数記憶させる記憶手段からなる生体成分の定量装置であり、生体成分を定量する際に、少なくとも、前記受発光プローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、前記受発光プローブを当接させた部分の皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階を経て目的とする生体成分濃度を定量することを特徴とするものである。
【００２５】
本発明の請求項６に係る生体成分の定量装置は、受発光間隔を２ｍｍ以下に設定した受発光プロ−ブ３と、前記受発光間隔を選択できる受発光間隔選択手段と、前記受発光プローブに近赤外光を供給する光源１と、前記受発光プローブ３で受光された近赤外光を分光する分光手段と、前記分光手段で分光されたスペクトル光を受光する受光手段と、前記受光手段からの信号を用い生体成分の定量を行う演算手段、皮膚性状に対応する検量式を複数記憶させる記憶手段からなる生体成分の定量装置であり、生体成分を定量する際に、少なくとも、前記受発光プローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、前記受発光プローブを当接させた部分の皮膚性状を分類する段階、分類した皮膚性状に対応する前記プローブの受発光間隔を選択する段階、分類した皮膚性状と受発光間隔に対応する検量式を選択する段階、設定した受発光間隔のプローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階を経て目的とする生体成分濃度を定量することを特徴とするものである。
【００２６】
【発明の実施の形態】
以下に本発明の具体的な実施の形態を示す。
【００２７】
本発明の生体成分の定量装置は、外観上、図１に示す従来のものと同様に形成されるものであって、ハロゲンランプ等からなる光源１と、光源１からの光を集光するための集光レンズ２と、集光レンズ２を通過した光を被測定物１０に照射すると共に被測定物１０を透過あるいは拡散反射した光を受光するための受発光プローブ３と、受光後の光を分光する回折格子を収めた回折格子ユニット４、前記回折格子ユニット４で分光された光（スペクトル光）を検出するためのＩｎＧａＡｓアレイ型受光素子ユニット５、前記ＩｎＧａＡｓアレイ型受光素子ユニット５で得られた信号をもとにグルコース濃度を演算するノート型パーソナルコンピュータ等の演算ユニット６とから構成されている。すなわち、上記の回折格子ユニット４が受発光プローブ３で受光された近赤外光を分光する分光手段として、上記のＩｎＧａＡｓアレイ型受光素子ユニット５が分光手段で分光されたスペクトル光を受光する受光手段として、上記の演算ユニット６が受光手段からの信号を用い生体成分の定量を行う演算手段及び皮膚性状に対応する検量式を複数記憶させる記憶手段として、それぞれ形成されている。そして、上記の演算ユニット６を使用することにより、皮膚性状を判断して分類する段階、分類した皮膚性状に対応する検量式を選択する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階がそれぞれ行なわれるものである。
【００２８】
また、皮膚に当接させ近赤外スペクトルを測定する受発光プローブ３は上記と同様に形成されるものであって、受発光プローブ３の受発光間隔Ｌ、すなわち、受発光プローブ３に設けた発光ファイバー７と受光ファイバー８の間隔Ｌは２ｍｍ以下に設定されている。
【００２９】
また、図１５に示す他の受発光プローブ３は、中心の受光ファイバー８の同心円状に受発光間距離（受発光間隔Ｌ１）３００μｍに複数本（６本）の発光ファイバー７ａと受発光間距離（受発光間隔Ｌ２）６５０μｍに複数本（１２本）の発光ファイバー７ｂを配置したものを用いた。受発光間距離３００μｍと６５０μｍからの発光は切り替え手段により切り替えられるようになっている。この切り替え手段が受発光間隔を選択できる受発光間隔選択手段として形成されるものであって、具体的には、光源１からの光をコリメートレンズにより平行光とした後、ハーフミラーにより３００μｍと６５０μｍの光ファイバー７ａ、７ｂへ１：２の割合で光量が分割されるように調整する。前記ハーフミラーで分割された光はそれぞれ集光ミラーや集光レンズで光ファイバー７ａ、７ｂの測定端と反対の端部に集光される。また、受発光間隔Ｌ１、Ｌ２の選択はこのハーフミラーと集光ミラーや集光レンズの間に配置した２つのメカニカルシャッターの開閉で行う。すなわち、３００μｍの受発光間隔Ｌ１を選択する時は３００μｍ側のメカニカルシャッターを開き、６５０μｍ側を閉じる。また、６５０μｍの受発光間隔Ｌ２を選択する時は６５０μｍ側のメカニカルシャッターを開き、３００μｍ側を閉じる。尚、この受発光プローブ３においても上記の受発光間距離（受発光間隔Ｌ１、Ｌ２）は３００μｍや６５０μｍに限定されるものではなく、上記と同様に２ｍｍ以下で皮膚性状（皮膚厚み）に応じて適宜に設定することができる。また、上記の受発光間隔Ｌは０．２ｍｍ以上に設定するのが好ましい。
【００３０】
一人の被験者（Ｐ１）から得られた皮膚（皮膚組織）の近赤外スペクトルを示すグラフを図３に、検量式の回帰係数を示すグラフを図４に示す。図３の近赤外スペクトルは図１の装置を用いて測定した。図４に示す検量式の回帰係数は経口でグルコース負荷を行い人為的に血糖値を変動させる実験を６回行い、実験で得られた皮膚の近赤外スペクトルと血液から測定した血糖値からなるデータを多変量解析することにより求めた。多変量解析はＰＬＳ回帰分析手法を用いている。また、解析に用いた近赤外スペクトルは微分処理あるいはＭＳＣ処理等の前処理をほどこさない生データを用いたが、これに限るものではない。被験者の性別は男、年齢は４１才、健常人である。皮膚の近赤外スペクトルは左手前腕内側部で行い、測定を行った部分の皮膚厚さは表皮層１２と真皮層１３の厚さの合計は約１．５ｍｍであった。皮膚厚さの測定はＣＯＲＴＥＸ社製（デンマーク）の超音波断層測定装置ＤｅｒｍａＳｃａｎ Ｃ Ｖｅｒ．３を用いて行った。
グルコース濃度、すなわち血糖値測定に用いられる検量式は
ｘｎ：波長ｎにおける吸光度
ｂ：定数
で表される。したがって、この検量式の回帰係数によりグルコース濃度の定量に重要な波長が推定できる。図４に検量線の回帰係数のグラフを示す。図４の回帰係数は１６００ｎｍ付近に特徴的な正のピークを有することがわかる。代表的な生体成分の吸光スペクトルを図５に示す。図５に示されるようにグルコース分子が１６００ｎｍ付近に他の生体成分に比較して特異な吸収を有することがわかる。このことから、上記実験で得られた検量式の回帰係数が生体皮膚組織中のグルコース分子の特異吸収を利用してグルコース濃度を定量していることがわかる。
【００３１】
上記の生体皮膚組織中グルコース濃度（血糖値）測定は上記従来例とほぼ同様であって、図６に示したような中心に配置した受光ファイバー８と、その受光ファイバー８を中心とする半径６５０μｍの円周上に等間隔で配置した１２本の発光ファイバー７からなる受発光プローブ３を皮膚面に９．８〜４９ｋＰａ（１００〜５００ｇｆ／ｃｍ2）、好ましくは２９．４ｋＰａ（３００ｇｆ／ｃｍ2）の一定加重で当接させ、近赤外スペクトルを測定することで行う。光源１からの光は円周上に配置された１２本の発光ファイバー７より照射され、皮膚組織を透過した光の一部が中心の受光ファイバー８で受光される。受光ファイバー８で受光される光の経路はＵ字型の経路Ａをとっていると推定でき、皮膚組織への到達深さは受光ファイバー８と発光ファイバー７間の距離で決まる。
【００３２】
人間の皮膚組織構造は、角質層を含む表皮層（表皮組織）１２、真皮層（真皮組織）１３、皮下組織層（皮下組織）１４の３層で構成される。真皮層１３には毛細血管、リンパ、神経組織が発達している。皮下組織層１４は主に脂肪組織で構成される。グルコースは水溶性のため脂肪組織中には存在しにくい。また、表皮層１２には毛細血管が発達していないため、血液で輸送されるグルコースが到達しにくい。そのため、血糖値の測定を皮膚の近赤外スペクトル測定を行うためには、真皮層１２の近赤外スペクトルを選択的に測定する必要がある。
【００３３】
実験に用いた受発光プローブ３は、受発光間隔Ｌを６５０μｍに設定したもので、表皮層１２と真皮層１３をあわせた皮膚厚さが１．５ｍｍである前記被験者から作成した検量式である。ある個体に対して作成された検量式は、その個体に適切であっても、皮膚厚さ等の皮膚性状に個体差があるため、そのまま他の個体に適用できるとは限らない。従来の他研究において個人毎の検量線が使われたのはそのためである。しかしながら、個体差があっても定量を目的とする生体成分（上記実験ではグルコース）は個体によらず共通であることから、その生体成分の変動（上記実験ではグルコースの変動）を捉えた検量式を用いれば、ある程度一般化された検量式を作成できることが期待できる。本発明に用いた検量式の一つは図４に示すような回帰係数形状を有する。この回帰係数を持つ検量式は、定量を目的とするグルコース分子の特異吸収波長である１６００ｎｍ付近の吸収を正のピークとして利用していることから、皮膚組織内のグルコース分子の変動を定量していることが期待できる。
【００３４】
次に、この検量式を用いて被験者（Ｐ１）の自身の血糖値変動を測定した結果を図７に示す。
【００３５】
同じ検量式を用いて他の４人の被験者（Ｐ２〜Ｐ５）の血糖値変動を測定する実験を合わせて実施した。被験者（Ｐ１）で作成した検量式を用いて、他の４人の被験者（Ｐ２〜Ｐ５）の血糖値変動を推定した場合の結果を図８に示す。また、この実験を行った被験者を含む被験者５名の属性を示す表を表１に示す。
【００３６】
【表１】
図７のグラフから、従来の技術と同様に、被験者（Ｐ１）から作成した検量式を用いて、同じ被験者（Ｐ１）の血糖値変動を推定する場合は相関係数０．９程度の良好な推定が可能であることがわかる。
【００３８】
次に、図８からわかるように、被験者（Ｐ１）で作成した検量式を用いて他の被験者の血糖値変動を推定した場合、被験者Ｐ２の血糖値変動については良好な測定が可能であるが、他の被験者Ｐ３、Ｐ４、Ｐ５についての血糖値測定は困難であることがわかった。図９に被験者（Ｐ１）で作成した検量式を用いて５名の被験者（Ｐ１〜Ｐ５）の血糖値変動を推定した結果のグラフを示す。この場合の相関係数は０．６９であった。
【００３９】
次に、被験者Ｐ２〜Ｐ５の被験者についても同様な６回のグルコース負荷実験を実施し、被験者毎に検量式を作成した。被験者（Ｐ１）で作成した検量式を含む６個の検量式を用いて、５名の被験者相互に血糖値変動の推定を行った。この結果を表２に示す。
【００４０】
【表２】
表２は作成した検量式に他の被験者毎のグルコース負荷実験で得られた近赤外スペクトルデータを代入し、血糖値変動を推定し、実測した血糖値と比較した結果をまとめたものである。この比較において相関係数（６回の当てはめの相関係数の平均値）が０．７以上の場合を○印、０．６以上０．７未満を△印、０．６未満を×印で示している。表２から分かるように、５名の被験者を用いたこの実験では、大きく２つの属性（Ａ群およびＢ群）に分かれ、その群内においてお互いに作成した検量式で血糖値の変動が測定できることが判明した。
【００４２】
表２と表１の結果を比較すると２群は皮膚厚さで層別できると考えられ、被験者の属性を分類する指標の一つとして皮膚厚さがあることが分かる。表１の皮膚厚さを求めた超音波断層写真を図１０に示す。皮膚厚さの値は、この断層写真より任意に選んだ３箇所の皮膚厚さ（表皮層と真皮層を合わせたもの）の平均値で算出した。
【００４３】
近赤外分光分析において精度よい検量式を作成するには、通常、多くのデータ蓄積をする必要がある。したがって、上記の血糖値測定においても検量式の作成には数日におよび被験者を拘束しデータ収集を行う必要があった。非侵襲的に、かつ、迅速に血糖値が測定できる近赤外線による定量手法の恩恵ははかり知れないが、そのための検量式作成に利用者の数日におよぶ拘束が必要とするのであれば、実用上の大きな問題と言わざるをえない。上記実験結果に示すように、我々は、利用者の個人毎で検量式を作成するのでなく、生体成分の測定対象である患者や被験者をある特性毎に分類し、分類した群に対応する検量式を前もって用意することで煩雑な検量式の作成プロセスを必要としない生体成分の定量方法および装置を提供できることがわかった。また、その分類は、皮膚厚さが大きな要因であることがわかった。また、予め複数の検量式を用意し、分類した群毎に適切な検量式を用いることで測定精度の向上がはかれることがわかった。
【００４４】
被験者の属性を分類する指標としては皮膚厚さ以外に、皮膚の角質水分量、皮膚組織の水分濃度、皮膚組織の密度、皮膚の色、皮膚の表面粗さ、性別、年齢、人種を使用することができる。
【００４５】
上記説明では定量分析を行う生体成分として主にグルコースあるいは血糖値を記述しているが、本発明はこれに限らず、組織水分量、中性脂肪、コレステロール、ＨｂＡ１ｃ（糖化ヘモグロビン）、フルクトサミン、アルブミン、グロブリン、尿酸等の生体成分の定量分析に利用できる。また、これらを代用特性とする皮膚の健康度や皮膚年齢、老化度、皮膚の張力（張り）等の測定に利用できる。
【００４６】
【実施例】
−参考例−
本参考例は超音波断層撮影装置を利用し、被験者の皮膚厚さを分類した例である。皮膚厚さの測定はＣＯＲＴＥＸ社製（デンマーク）の超音波断層測定装置ＤｅｒｍａＳｃａｎＣ Ｖｅｒ．３を用いて行った。皮膚厚さの測定結果を図１０に示す。
【００４７】
本参考例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。近赤外スペクトルから血糖値を推定するための演算手段では、表皮層と真皮層との合計の厚さが１．２ｍｍ以上の被験者の血糖値を推定する検量式と、１．２ｍｍ未満の被験者の血糖値を推定する検量式の２種類の検量式が準備されている。
【００４８】
表皮層と真皮層との合計の厚さが１．２ｍｍ以上の被験者の血糖値を推定する検量式は、
【００４９】
また、表皮層と真皮層との合計の厚さが１．２ｍｍ未満の被験者の血糖値を推定する検量式は、
【００５０】
尚、ａｎとｃｎは波長ｎにおける回帰係数、ｘｎは波長ｎにおける吸光度、ｂとｄは定数をそれぞれ示す。
【００５１】
皮膚の近赤外スペクトルによる血糖値の推定と、この測定値の測定精度を判断するための実測値の測定は以下の手順で行った。
（１）超音波断層撮影装置を利用し、被験者の皮膚厚さを測定する。
【００５２】
皮膚厚さの測定は、近赤外スペクトルを測定する同一部位である左手の前腕内側で行った。
（２）測定した皮膚厚さに基づいて、被験者の分類を行う。
（３）被験者の分類に応じた検量式を選択する。
（４）血糖値の実測を行う。
【００５３】
この段階は実際の測定では必要ないが、本参考例での測定精度を検証するために実施した。
（５）皮膚の近赤外スペクトルを測定する。
（６）測定した近赤外スペクトルを選択した検量式に代入し、血糖値を推定する。
【００５４】
上記手順は被験者が一人に固定される場合、２回目以降は皮膚厚さを測定する段階は省略し、分類された群の測定条件、検量式をそのまま用いてもよい。
【００５５】
本参考例での血糖定量の結果として、皮膚厚さ１．２ｍｍ以上と判断され、血糖値を推定された群の測定結果を図１１に、皮膚厚さ１．２ｍｍ未満と推定され、血糖値を推定された群の測定結果を図１２に示す。また、対照例として、同一データを被験者を分類することなく、一律に皮膚厚さ１．２ｍｍの血糖値推定に用いた検量式で血糖値推定を行った場合の結果を図９に示す。被験者を分類せず血糖値推定を行った場合の相関係数は０．６９で、被験者の分類を行った場合の相関係数は０．９３と０．８５であることから、本参考例による定量精度の向上が確認できた。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
【００５６】
本参考例では、皮膚厚さを実測する手段として超音波断層装置を用いたが、本発明はこれに限定するものではなく、断層像として計測しない超音波測定手法を用いても構わない。また、本発明は光コヒーレンストモグラフィを用いて皮膚厚さの実測を行うようにする。
【００５７】
超音波による皮膚厚さの測定は画像を用いて行う手法と、画像化せずに超音波の反射により距離を直接求める手法がある。従って、上記参考例のように断層画像を用いなくても超音波の反射信号を用いて皮膚厚さを求めることができる。
【００５８】
また、光コヒーレンストモグラフィーの代表的なものとして、低コヒーレンス干渉法と光周波数走査法がある。いずれの手法も数十μｍオーダーの分解能で生体の断層測定を可能としている。そして、近赤外レーザー光を前腕内側部に照射し、得られた皮膚断層画像から皮膚厚さを測定することができるものである。
−実施例１−
本実施例は測定した近赤外スペクトルより被験者の分類を行った例である。本実施例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。近赤外スペクトルの測定は、左手の前腕内側で行った。また、測定に用いた受発光プローブは従来の技術に示したものと同じ受発光ファイバー間が６５０μｍのものを用いた。
【００５９】
近赤外スペクトルを用いて行う被験者の分類は、主成分分析により行った。図１３に主成分分析により被験者の分類を行った結果を示す。図１３は、表皮層と真皮層との合計の厚さが１．２ｍｍ以上の被験者と、１．２ｍｍ未満の被験者が第２主成分（ＰＣ２）により明確に分類できることを示している（図１３の上部と下部）。
【００６０】
図１４に第２主成分のローディングプロットを示す。第２主成分のローディングプロットは波長１７２８ｎｍと１７６９ｎｍに特徴的なピークを有し、このピークは図５に示した中性脂肪のピークと一致する。したがって、第２主成分の寄与の大きさは皮下組織（脂肪組織）への光の到達度と解釈でき、皮膚厚さの代用特性として主成分分析の分類が利用できることを示している。すなわち、皮膚厚さの異なる個体差は皮下組織への光の到達量の違いとなり、最終的に近赤外スペクトルの形状差として表れる。この近赤外スペクトルの形状差を主成分分析等の多変量解析を用いて分類することで、被験者の分類が可能である。多変量解析による分類は主成分分析に限らず、クラスタ分析、判別分析等の定性分析手法を利用することができる。
【００６１】
血糖値を推定する検量式は、参考例と同様に２種類準備されている。皮膚の近赤外スペクトルによる血糖値の推定と、この測定値の測定精度を判断するための実測値の測定は以下の手順で行った。
（１）近赤外スペクトルの測定を行う。
（２）測定した近赤外スペクトルの主成分分析を行い、被験者の分類を行う。
（３）被験者の分類に応じた検量式を選択する。
（４）測定した近赤外スペクトルを選択した検量式に代入し、血糖値を推定する。
（５）血糖値の実測を行う。
【００６２】
この段階は実際の測定では必要ないが、本実施例での測定精度を検証するために実施した。
【００６３】
上記手順は被験者が一人に固定される場合、２回目以降は主成分分析する段階は省略し、分類された群の測定条件、検量式をそのまま用いてもよい。
【００６４】
本実施例での血糖定量の結果として、皮膚厚さ１．２ｍｍ以上と判断され、血糖値を推定された群の測定結果を図１１に、皮膚厚さ１．２ｍｍ未満と推定され、血糖値を推定された群の測定結果を図１２に示す。また、対照例として、同一データを被験者を分類することなく、一律に皮膚厚さ１．２ｍｍの血糖値推定に用いた検量式で血糖値推定を行った場合の結果を図９に示す。被験者を分離せず血糖値推定を行った場合の相関係数は０．６９で、被験者の分類を行った場合の相関係数は０．８５と０．９３であることから、本発明による定量精度の向上が確認できた。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
−実施例２−
本実施例は予備的に測定した近赤外スペクトルより被験者の分類を行った例である。
【００６５】
本実施例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。予備的な近赤外スペクトルの測定は、血糖値を推定するために近赤外スペクトルを測定する同一部位である左手の前腕内側で行った。また、測定に用いた受発光プローブは図１５に示すように、中心の受光ファイバー８の同心円状に受発光間距離３００μｍに６本の発光ファイバー７ａと受発光間距離６５０μｍに１２本の発光ファイバー７ｂを配置したものを用いた。受発光間距離３００μｍと６５０μｍからの発光は切り替え手段により切り替えられるようになっており、予備的な近赤外スペクトルの測定は受発光間距離６５０μｍを用いて行った。
【００６６】
予備的な近赤外スペクトルを用いて行う被験者の分類は、実施例１と同様に主成分分析により行った。本実施例では、主成分分析によって表皮層と真皮層との合計の厚さが１．２ｍｍ未満と判断した場合は、発光の切り替え手段によって受発光間隔３００μｍの近赤外スペクトルを測定する。これは、皮膚厚さが薄い群については真皮層を適切に透過する光路を設定した方がより血糖値の推定精度が向上するためで、受発光間距離を狭めることで皮膚の深さ方向への到達度を浅くする効果がある。
【００６７】
血糖値を推定する検量式は、他の実施例と同様に２種類準備されている。皮膚の近赤外スペクトルによる血糖値の推定と、この測定値の測定精度を判断するための実測値の測定は以下の手順で行った。
▲１▼ 予備的な皮膚の近赤外スペクトルの測定を受発光間隔６５０μｍで行う。
▲２▼ 測定した皮膚厚さに基づいて、被験者の分類を行う。
▲３▼ 被験者の分類に応じた検量式を選択する。
▲４▼ 血糖値の実測を行う。
【００６８】
この段階は実際の測定では必要ないが、本実施例での測定精度を検証するために実施した。
▲５▼ 皮膚厚さを１．２ｍｍ未満と判断した場合、受発光間距離３００μｍで皮膚の近赤外スペクトルを測定する。皮膚厚さを１．２ｍｍ以上と判断した場合、▲１▼で測定した皮膚の近赤外スペクトルをそのまま利用する。
▲６▼ 測定した近赤外スペクトルを選択した検量式に代入し、血糖値を推定する。
上記手順は被験者が一人に固定される場合、２回目以降は予備的な皮膚の近赤外スペクトルを測定する段階は省略し、分類された群の測定条件、検量式をそのまま用いてもよい。
【００６９】
本発明により従来の技術のような単一の測定条件でなく、個体差に応じて受発光ファイバー間隔を設定し、その分類に応じた検量式を準備することで定量精度の向上が可能である。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
−実施例３−
本実施例は測定した近赤外スペクトルより被験者の分類を行った例である。本実施例に用いた皮膚の近赤外スペクトルにより血糖値を推定するための装置は図１に示した装置とほぼ同様の装置を用いた。近赤外スペクトルの測定は、左手の前腕内側で行った。また、測定に用いた受発光プローブは従来の技術に示したものと同じ受発光ファイバー間が６５０μｍのものを用いた。
【００７０】
近赤外スペクトルを用いて行う被験者の分類は、近赤外スペクトルの１７２８ｎｍでの吸光度（Ａｂ１７２８）と１６７０ｎｍでの吸光度（Ａｂ１６７０）の比を用いて行った。すなわち、
Ａｂ１７２８／Ａｂ１６７０＞αの時は皮膚厚さ１．２ｍｍ未満
Ａｂ１７２８／Ａｂ１６７０＜α，Ａｂ１７２８／Ａｂ１６７０＝αの時は皮膚厚さ１．２ｍｍ以上
ここでαは定数
と判断し、皮膚性状の分類を行い、他の実施例と同様に、事前に準備した皮膚性状に対応する検量式に近赤外スペクトルを代入して血糖値の算出を行う。
【００７１】
本発明により、近赤外スペクトルによる皮膚性状の分類を脂肪成分由来の特異吸収波長を利用して行その分類に応じた検量式を準備することで定量精度の向上が可能である。また、個人毎に検量式作成することなく血糖値の推定が可能となった。
【００７２】
【発明の効果】
従来の技術では、個体差が大きい人間や生物の皮膚組織中の生体成分濃度を測定する場合に、通常、生体成分の定量を行う個体毎、あるいは測定部分毎に検量式を作成し定量分析を行うことが行われるが、本発明のように、被験者を皮膚性状に分類し、その分類毎に検量式を準備することで、個体差が大きい人間や生物の皮膚組織における生体成分濃度の定量、あるいはその成分濃度を代用特性とした特性値の定量を、個体毎、あるいは測定部分毎に検量式を作成することなくできるようになった。それにより、個体毎あるいは測定部分毎に検量線を作成する場合での、測定対象である人や生物を長時間拘束する必要がなくなり、測定対象の人や生物の負担を軽減することができた。さらに、測定者にも検量のための式を作成する等の複雑な作業の必要がなくなり、負担を軽減することができた。
【００７３】
また、個体差を考慮する際に、その個体の属性を判断し最適な測定条件や検量式を選択することで測定精度を向上させることもができた。
【図面の簡単な説明】
【図１】 本発明及び従来例を示す概略図である。
【図２】 同上の皮膚組織の断面模式図及び真皮のスペクトル測定を示す概念図である。
【図３】 同上の前腕内側部で測定した皮膚の近赤外スペクトルを示すグラフである。
【図４】 同上の血糖値定量を行う検量式の回帰係数を示すグラフである。
【図５】 同上の生体成分の吸収スペクトルを示すグラフである。
【図６】 同上の受発光プローブを示す正面図である。
【図７】 従来技術による血糖値測定結果を示すグラフである。
【図８】 本発明の一つの検量式での複数の被験者の血糖値測定を行った結果を示すグラフである。
【図９】 同上の一つの検量式での複数の被験者の血糖値測定を行った結果を示すグラフである。
【図１０】 同上の被験者の皮膚厚さを求めた超音波断層写真である。
【図１１】 同上の皮膚厚さ１．２ｍｍ以上の群に対して血糖値変動を求めた結果を示すグラフである。
【図１２】 同上の皮膚厚さ１．２ｍｍ未満の群に対して血糖値変動を求めた結果を示すグラフである。
【図１３】 同上の皮膚の近赤外スペクトルを主成分分析した結果を示すグラフである。
【図１４】 同上の主成分分析の第２主成分のローディングプロットである。
【図１５】 同上の実施例２に用いた受発光プローブを示す正面図である。
【符号の説明】
１ 光源
３ 受発光プローブ
６ 演算ユニット
Claims (6)
- 皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニットにより、目的とする生体成分濃度を定量する定量装置において、前記演算ユニットにより、皮膚厚さを指標として皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有し、前記皮膚厚さを指標として皮膚性状の分類を行うにあたって、皮膚の近赤外スペクトルの形状を実際の皮膚厚さの代用特性とすることを特徴とする生体成分の定量装置。
- 前記皮膚の近赤外スペクトルの形状を実際の皮膚厚さの代用特性として行う手法として、皮下脂肪の特異吸収波長である１７００〜１８００ｎｍの波長範囲から少なくとも一つの波長を選択し、その吸光度から皮膚性状の分類を行うことを特徴とする請求項１記載の生体成分の定量装置。
- 皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニットにより、目的とする生体成分濃度を定量する定量装置において、前記演算ユニットにより、皮膚厚さを指標として皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有し、前記皮膚厚さを指標として皮膚性状の分類を行うにあたって、光コヒーレンストモグラフィを利用することを特徴とする生体成分の定量装置。
- 皮膚の近赤外スペクトルを用いて生体成分の定量を行う演算手段と皮膚性状に対応する検量式を複数記憶した記憶手段とからなる演算ユニットにより、目的とする生体成分濃度を定量する定量装置において、前記演算ユニットにより、皮膚の近赤外スペクトルで皮膚性状を分類する段階、分類した皮膚性状に対応する近赤外スペクトルの測定条件を設定する段階、分類した皮膚性状に対応する検量式を選択する段階、選択された検量式に前記近赤外スペクトルを代入演算し、目的とする生体成分濃度を定量する段階を有することを特徴とする生体成分の定量装置。
- 受発光間隔を２ｍｍ以下に設定した受発光プローブと、前記受発光プローブに近赤外光を供給する光源と、前記受発光プローブからのスペクトル光を受光する受光手段と、前記近赤外光を分光する分光手段と、前記受光手段からの信号を用い生体成分の定量を行う演算手段、皮膚性状に対応する検量式を複数記憶させる記憶手段からなる生体成分の定量装置であり、生体成分を定量する際に、少なくとも、前記受発光プローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、前記受発光プローブを当接させた部分の皮膚性状を分類する段階、分類した皮膚性状に対応する検量式を選択する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階を経て目的とする生体成分濃度を定量することを特徴とする生体成分の定量装置。
- 受発光間隔を２ｍｍ以下に設定した受発光プロ−ブと、前記受発光間隔を選択できる受発光間隔選択手段と、前記受発光プローブに近赤外光を供給する光源と、前記受発光プローブからのスペクトル光を受光する受光手段と、前記近赤外光を分光する分光手段と、前記受光手段からの信号を用い生体成分の定量を行う演算手段、皮膚性状に対応する検量式を複数記憶させる記憶手段からなる生体成分の定量装置であり、生体成分を定量する際に、少なくとも前記受発光プローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、前記受発光プローブを当接させた部分の皮膚性状を分類する段階、分類した皮膚性状に対応する前記プローブの受発光間隔を選択する段階、分類した皮膚性状と受発光間隔に対応する検量式を選択する段階、設定した受発光間隔のプローブを皮膚表面に当接させ近赤外スペクトルを測定する段階、測定した前記近赤外スペクトルを選択された前記検量式に代入演算し、目的とする生体成分濃度を定量する段階を経て目的とする生体成分濃度を定量することを特徴とする生体成分の定量装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001350746A JP3931638B2 (ja) | 2001-11-15 | 2001-11-15 | 生体成分の定量装置 |
| EP20020803116 EP1452135B1 (en) | 2001-11-15 | 2002-11-14 | Method for calculating a biological component density of a subject |
| US10/478,276 US7333841B2 (en) | 2001-11-15 | 2002-11-14 | Method and device for calculating a biological component density of a subject |
| PCT/JP2002/011871 WO2003041582A1 (en) | 2001-11-15 | 2002-11-14 | Method of measuring biological component concentration and apparatus therefor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001350746A JP3931638B2 (ja) | 2001-11-15 | 2001-11-15 | 生体成分の定量装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003144421A JP2003144421A (ja) | 2003-05-20 |
| JP3931638B2 true JP3931638B2 (ja) | 2007-06-20 |
Family
ID=19163183
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001350746A Expired - Lifetime JP3931638B2 (ja) | 2001-11-15 | 2001-11-15 | 生体成分の定量装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US7333841B2 (ja) |
| EP (1) | EP1452135B1 (ja) |
| JP (1) | JP3931638B2 (ja) |
| WO (1) | WO2003041582A1 (ja) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH064222B2 (ja) | 1985-09-18 | 1994-01-19 | 株式会社岡本工作機械製作所 | 吸振器式砥石カバ− |
| WO2011013694A1 (ja) | 2009-07-28 | 2011-02-03 | パナソニック電工株式会社 | 血糖値推定装置 |
| JP2012085878A (ja) * | 2010-10-20 | 2012-05-10 | Seiko Epson Corp | 濃度定量装置及び濃度定量方法並びにプログラム |
| JP2012085877A (ja) * | 2010-10-20 | 2012-05-10 | Seiko Epson Corp | 濃度定量装置及び濃度定量方法並びにプログラム |
| JP2013138811A (ja) * | 2012-01-06 | 2013-07-18 | Seiko Epson Corp | 濃度定量装置、濃度定量方法及びプログラム |
| JP2013140126A (ja) * | 2012-01-06 | 2013-07-18 | Seiko Epson Corp | 濃度定量装置、濃度定量方法及びプログラム |
| US9265425B2 (en) | 2011-11-11 | 2016-02-23 | Hideo Ando | Measuring method of life activity, measuring device of life activity, transmission method of life activity detection signal, or service based on life activity information |
| US9456776B2 (en) | 2013-05-07 | 2016-10-04 | Hideo Ando | Detection method of life activity, measuring device of life activity, transmission method of life activity detection signal, or service based on life activity information |
| US9464983B2 (en) | 2010-07-12 | 2016-10-11 | Seiko Epson Corporation | Concentration determination apparatus, probe, concentration determination method, and program |
Families Citing this family (43)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7299080B2 (en) * | 1999-10-08 | 2007-11-20 | Sensys Medical, Inc. | Compact apparatus for noninvasive measurement of glucose through near-infrared spectroscopy |
| JP4533603B2 (ja) * | 2003-07-14 | 2010-09-01 | ポーラ化成工業株式会社 | 皮膚状態の測定方法 |
| JP4348467B2 (ja) | 2003-09-30 | 2009-10-21 | 学校法人日本大学 | 光を用いた生体の物質特性検査装置 |
| US7251517B2 (en) | 2004-06-30 | 2007-07-31 | Hitachi, Ltd. | Blood sugar level measuring apparatus |
| US7215983B2 (en) | 2004-06-30 | 2007-05-08 | Hitachi, Ltd. | Blood sugar level measuring apparatus |
| JP4586680B2 (ja) * | 2004-08-25 | 2010-11-24 | パナソニック電工株式会社 | 体内成分の定量分析用検量線の作成方法、および同検量線を用いた定量分析装置 |
| DE602005022388D1 (de) | 2004-08-25 | 2010-09-02 | Panasonic Elec Works Co Ltd | Quantitatives Analysegerät, das eine Kalibrierungskurve benutzt |
| US7856985B2 (en) | 2005-04-22 | 2010-12-28 | Cynosure, Inc. | Method of treatment body tissue using a non-uniform laser beam |
| KR20080050404A (ko) * | 2005-09-02 | 2008-06-05 | 포라 가세이 고교 가부시키가이샤 | 피부의 상태의 판정 방법 및 피부의 두께의 예측 방법 |
| US20070232872A1 (en) * | 2006-03-16 | 2007-10-04 | The Board Of Regents Of The University Of Texas System | Continuous noninvasive glucose monitoring in diabetic, non-diabetic, and critically ill patients with oct |
| US8485016B2 (en) | 2006-06-12 | 2013-07-16 | Trefford Simpson | Method of calibrating a biometric device |
| JP2008005988A (ja) * | 2006-06-28 | 2008-01-17 | Jasco Corp | 血糖値分析方法、及び血糖値分析装置 |
| US7586957B2 (en) | 2006-08-02 | 2009-09-08 | Cynosure, Inc | Picosecond laser apparatus and methods for its operation and use |
| JP5684443B2 (ja) * | 2007-06-06 | 2015-03-11 | 横河電機株式会社 | 生体成分測定装置 |
| JP2010005047A (ja) * | 2008-06-25 | 2010-01-14 | Panasonic Electric Works Co Ltd | 光学式生体成分測定用の検量モデルの作成方法及び非侵襲血糖値測定装置 |
| FR2940904B1 (fr) * | 2009-01-13 | 2012-08-31 | Urgo Laboratoires | Systeme de mesure de pression d'interface |
| JP4715930B2 (ja) * | 2009-02-03 | 2011-07-06 | パナソニック電工株式会社 | 光学式皮下脂肪厚測定装置 |
| JP2012519554A (ja) * | 2009-03-05 | 2012-08-30 | サイノシュア・インコーポレーテッド | 熱的外科手術安全装置および熱的外科手術方法 |
| WO2012158007A1 (es) * | 2011-05-18 | 2012-11-22 | Instituto Politecnico Nacional | Dispositivo láser para la medición de espesores de la epidermis |
| JP5838517B2 (ja) * | 2011-06-21 | 2016-01-06 | セイコーエプソン株式会社 | 濃度定量装置、濃度定量方法 |
| JP2013103094A (ja) * | 2011-11-16 | 2013-05-30 | Sony Corp | 測定装置、測定方法、プログラム及び記録媒体 |
| WO2013141734A1 (ru) * | 2012-03-20 | 2013-09-26 | Khokhoev Elbrus Marklenovich | Способ и устройство неинвазивного контроля уровня глюкозы в крови |
| WO2013158299A1 (en) | 2012-04-18 | 2013-10-24 | Cynosure, Inc. | Picosecond laser apparatus and methods for treating target tissues with same |
| US8843186B2 (en) | 2012-11-21 | 2014-09-23 | Folim G. Halaka | Non-invasive reagentless glucose determination |
| US9693694B2 (en) | 2012-11-21 | 2017-07-04 | Folim G. Halaka | Cancer cell detection using dielectrophoretic dynamic light scattering (DDLS) spectroscopy |
| JP5511033B1 (ja) | 2012-12-04 | 2014-06-04 | Ｎｅｃシステムテクノロジー株式会社 | 血糖値予測装置、測定装置、血糖値予測方法、及びプログラム |
| EP2973894A2 (en) | 2013-03-15 | 2016-01-20 | Cynosure, Inc. | Picosecond optical radiation systems and methods of use |
| JP6387610B2 (ja) * | 2013-12-27 | 2018-09-12 | ミツミ電機株式会社 | 生体情報測定装置 |
| JP5911113B2 (ja) * | 2014-02-05 | 2016-04-27 | Ｎｅｃソリューションイノベータ株式会社 | 測定装置 |
| US20160008556A1 (en) * | 2014-07-10 | 2016-01-14 | Elwha Llc | Needle insertion system and method for inserting a movable needle into a vertebrate subject |
| WO2016054079A1 (en) | 2014-09-29 | 2016-04-07 | Zyomed Corp. | Systems and methods for blood glucose and other analyte detection and measurement using collision computing |
| JP2016075522A (ja) * | 2014-10-03 | 2016-05-12 | セイコーエプソン株式会社 | 目的成分検量装置及び電子機器 |
| US9554738B1 (en) | 2016-03-30 | 2017-01-31 | Zyomed Corp. | Spectroscopic tomography systems and methods for noninvasive detection and measurement of analytes using collision computing |
| CN105866067A (zh) * | 2016-06-16 | 2016-08-17 | 江西出入境检验检疫局检验检疫综合技术中心 | 一种基于近红外光谱技术的茶粕中茶皂素快速测定方法 |
| KR102498120B1 (ko) | 2017-10-17 | 2023-02-09 | 삼성전자주식회사 | 생체정보 센서의 오차 보정 장치 및 방법과, 생체정보 추정 장치 및 방법 |
| KR102539144B1 (ko) | 2017-10-19 | 2023-06-01 | 삼성전자주식회사 | 개인화된 단위 스펙트럼 획득 장치 및 방법과, 생체 성분 추정 장치 및 방법 |
| CN112042066A (zh) | 2018-02-26 | 2020-12-04 | 赛诺秀股份有限公司 | 调q倾腔亚纳秒激光器 |
| KR102658237B1 (ko) * | 2018-10-02 | 2024-04-16 | 삼성전자주식회사 | 분석 물질의 혈중 농도 추정 장치 및 방법과, 모델 생성 장치 및 방법 |
| DK3809118T3 (da) * | 2019-10-17 | 2023-09-18 | Evonik Operations Gmbh | Fremgangsmåde til at forudsige en egenskabsværdi af et materiale ved anvendelse af hovedkomponentanalyse |
| CN115153533B (zh) * | 2020-02-26 | 2024-10-22 | 先阳科技有限公司 | 一种组织成分无创检测方法、装置、系统及可穿戴设备 |
| WO2022113608A1 (ja) * | 2020-11-24 | 2022-06-02 | サントリーホールディングス株式会社 | 健康支援装置、健康支援方法、および記録媒体 |
| JP7703112B2 (ja) * | 2022-08-09 | 2025-07-04 | バイオニクス株式会社 | バイタル情報測定装置、バイタル情報測定方法、生体管理システム、及び生体管理方法 |
| CN115963068B (zh) * | 2023-03-16 | 2023-05-05 | 北京心联光电科技有限公司 | 皮肤组织成分含量的测定方法及装置 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5460177A (en) * | 1993-05-07 | 1995-10-24 | Diasense, Inc. | Method for non-invasive measurement of concentration of analytes in blood using continuous spectrum radiation |
| JPH09215679A (ja) | 1996-02-09 | 1997-08-19 | Matsushita Electric Works Ltd | 体液成分濃度測定装置 |
| US5725480A (en) * | 1996-03-06 | 1998-03-10 | Abbott Laboratories | Non-invasive calibration and categorization of individuals for subsequent non-invasive detection of biological compounds |
| JP4212007B2 (ja) * | 1996-11-26 | 2009-01-21 | パナソニック電工株式会社 | 血液成分濃度の分析装置 |
| JP2000131322A (ja) | 1998-10-29 | 2000-05-12 | Matsushita Electric Works Ltd | グルコース濃度の定量方法及びその装置 |
| JP2000186998A (ja) * | 1998-12-22 | 2000-07-04 | Matsushita Electric Works Ltd | 生体スペクトル測定装置 |
| US6584340B1 (en) * | 1998-12-24 | 2003-06-24 | Matsushita Electric Industrial Co., Ltd. | Living body information measuring apparatus living body information measuring method body fat measuring apparatus body fat measuring method and program recording medium |
| AU754677B2 (en) * | 1999-01-22 | 2002-11-21 | Instrumentation Metrics, Inc. | System and method for noninvasive blood analyte measurements |
| US6587702B1 (en) * | 1999-01-22 | 2003-07-01 | Instrumentation Metrics, Inc | Classification and characterization of tissue through features related to adipose tissue |
| US6512937B2 (en) | 1999-07-22 | 2003-01-28 | Sensys Medical, Inc. | Multi-tier method of developing localized calibration models for non-invasive blood analyte prediction |
| JP2001212116A (ja) * | 2000-02-02 | 2001-08-07 | Matsushita Electric Ind Co Ltd | 血液成分測定装置 |
| US6738652B2 (en) * | 2000-06-15 | 2004-05-18 | Sensys Medical, Inc. | Classification and screening of test subjects according to optical thickness of skin |
-
2001
- 2001-11-15 JP JP2001350746A patent/JP3931638B2/ja not_active Expired - Lifetime
-
2002
- 2002-11-14 EP EP20020803116 patent/EP1452135B1/en not_active Expired - Lifetime
- 2002-11-14 US US10/478,276 patent/US7333841B2/en not_active Expired - Lifetime
- 2002-11-14 WO PCT/JP2002/011871 patent/WO2003041582A1/ja not_active Ceased
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH064222B2 (ja) | 1985-09-18 | 1994-01-19 | 株式会社岡本工作機械製作所 | 吸振器式砥石カバ− |
| WO2011013694A1 (ja) | 2009-07-28 | 2011-02-03 | パナソニック電工株式会社 | 血糖値推定装置 |
| US9464983B2 (en) | 2010-07-12 | 2016-10-11 | Seiko Epson Corporation | Concentration determination apparatus, probe, concentration determination method, and program |
| JP2012085878A (ja) * | 2010-10-20 | 2012-05-10 | Seiko Epson Corp | 濃度定量装置及び濃度定量方法並びにプログラム |
| JP2012085877A (ja) * | 2010-10-20 | 2012-05-10 | Seiko Epson Corp | 濃度定量装置及び濃度定量方法並びにプログラム |
| US9370323B2 (en) | 2010-10-20 | 2016-06-21 | Seiko Epson Corporation | Concentration determination apparatus, concentration determination method, and program |
| US9265425B2 (en) | 2011-11-11 | 2016-02-23 | Hideo Ando | Measuring method of life activity, measuring device of life activity, transmission method of life activity detection signal, or service based on life activity information |
| US9743837B2 (en) | 2011-11-11 | 2017-08-29 | Hideo Ando | Measuring method of life activity, measuring device of life activity, transmission method of life activity detection signal, or service based on life activity information |
| JP2013138811A (ja) * | 2012-01-06 | 2013-07-18 | Seiko Epson Corp | 濃度定量装置、濃度定量方法及びプログラム |
| JP2013140126A (ja) * | 2012-01-06 | 2013-07-18 | Seiko Epson Corp | 濃度定量装置、濃度定量方法及びプログラム |
| US9456776B2 (en) | 2013-05-07 | 2016-10-04 | Hideo Ando | Detection method of life activity, measuring device of life activity, transmission method of life activity detection signal, or service based on life activity information |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2003041582A1 (en) | 2003-05-22 |
| EP1452135B1 (en) | 2013-02-13 |
| US7333841B2 (en) | 2008-02-19 |
| EP1452135A1 (en) | 2004-09-01 |
| US20040142402A1 (en) | 2004-07-22 |
| JP2003144421A (ja) | 2003-05-20 |
| EP1452135A4 (en) | 2008-11-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3931638B2 (ja) | 生体成分の定量装置 | |
| JP4936203B2 (ja) | グルコース濃度の定量装置 | |
| US7343185B2 (en) | Measurement of body compounds | |
| KR100893432B1 (ko) | 조직 샘플 중의 표적 분석 물질 특성의 비침습적 측정방법 및 이를 위한 장치 | |
| US6675029B2 (en) | Apparatus and method for quantification of tissue hydration using diffuse reflectance spectroscopy | |
| US6668181B2 (en) | Method for quantification of stratum corneum hydration using diffuse reflectance spectroscopy | |
| US6587702B1 (en) | Classification and characterization of tissue through features related to adipose tissue | |
| US7725144B2 (en) | Determination of disease state using raman spectroscopy of tissue | |
| JP4472794B2 (ja) | グルコース濃度の定量装置 | |
| US20160354015A1 (en) | Non-invasive Measurement of Skin Thickness and Glucose Concentration with Raman Spectroscopy and Method of Calibration thereof | |
| JP2006126219A (ja) | 非侵襲性赤外分光法における多重スペクトル分析のための方法および装置 | |
| JP2002236097A (ja) | 非侵襲性近赤外分光法における多重スペクトル分析のための方法および装置 | |
| EP0630203A1 (en) | Non-invasive device and method for determining concentrations of various components of blood or tissue | |
| US20120078075A1 (en) | Determination of a measure of a glycation end-product or disease state using tissue fluorescence in combination with one or more other tests | |
| US8406839B2 (en) | Method and apparatus for determining blood analytes | |
| JPWO2003042180A1 (ja) | 生体成分濃度の測定方法及びその装置 | |
| US20010041829A1 (en) | Non-invasive method of determining skin thickness and characterizing layers of skin tissue in vivo | |
| Pleitez et al. | Infrared reflectometry of skin: Analysis of backscattered light from different skin layers | |
| Bodén et al. | Characterization of healthy skin using near infrared spectroscopy and skin impedance | |
| WO2019208561A1 (ja) | 血液成分の血中濃度測定方法、血中濃度測定装置およびプログラム | |
| JP2004321325A (ja) | 血糖値の定量方法 | |
| JPH10190A (ja) | 生体組織性状測定装置 | |
| JP4563075B2 (ja) | 血糖値測定装置 | |
| JP2000131322A (ja) | グルコース濃度の定量方法及びその装置 | |
| KR20240101725A (ko) | 캘리브레이션 방법 및 시스템 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041025 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061121 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070122 |
|
| 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: 20070220 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070305 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 3931638 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100323 Year of fee payment: 3 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100323 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: 20100323 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110323 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120323 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120323 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130323 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130323 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140323 Year of fee payment: 7 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| 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 |
|
| 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 |
|
| 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 |
|
| EXPY | Cancellation because of completion of term |