JP3604231B2 - グルコース濃度測定方法および装置 - Google Patents
グルコース濃度測定方法および装置 Download PDFInfo
- Publication number
- JP3604231B2 JP3604231B2 JP12179096A JP12179096A JP3604231B2 JP 3604231 B2 JP3604231 B2 JP 3604231B2 JP 12179096 A JP12179096 A JP 12179096A JP 12179096 A JP12179096 A JP 12179096A JP 3604231 B2 JP3604231 B2 JP 3604231B2
- Authority
- JP
- Japan
- Prior art keywords
- light
- glucose concentration
- backscattered
- intensity
- aqueous humor
- 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
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, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
-
- 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, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14532—Measuring 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
- G01B9/02002—Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies
- G01B9/02003—Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies using beat frequencies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
- G01B9/02002—Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies
- G01B9/02004—Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies using frequency scans
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
- G01B9/02007—Two or more frequencies or sources used for interferometric measurement
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
- G01B9/02012—Interferometers characterised by controlling or generating intrinsic radiation properties using temporal intensity variation
- G01B9/02014—Interferometers characterised by controlling or generating intrinsic radiation properties using temporal intensity variation by using pulsed light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02055—Reduction or prevention of errors; Testing; Calibration
- G01B9/02056—Passive reduction of errors
- G01B9/02057—Passive reduction of errors by using common path configuration, i.e. reference and object path almost entirely overlapping
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/0209—Low-coherence interferometers
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Ophthalmology & Optometry (AREA)
- Emergency Medicine (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Eye Examination Apparatus (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は生体におけるグルコース濃度の測定方法および装置に関し、詳細には眼球部の前眼房水内におけるグルコース濃度を非侵襲的に測定する方法および装置、並びにこの前眼房水内におけるグルコース濃度に基づいて血中グルコース濃度を非侵襲的に測定する方法および装置に関するものである。
【０００２】
【従来の技術】
血液中のグルコース濃度は個人差によりその平均レベルは異なるが、特に糖尿病疾患に対する投薬の要否を決定するための重要な指標値となっている。
【０００３】
ところでこの血中グルコース濃度は、食餌、肉体的活動、その他疾患の併発等によって極短時間の間に大きく変動する特性を有しており、急激な血中グルコース濃度の上昇によって緊急に投薬を要する場合も少なくない。
【０００４】
このため、このような疾患を有する患者についてはなるべく短い間隔で血中グルコース濃度をモニタすることが望まれているが、一方でこの血中グルコース濃度のモニタは通常、指先を穿切して実際に採血し、この血液を分析することにより血液中に含まれるグルコースの濃度（以下、単に血中グルコース濃度という）を測定することにより行なわれており、この穿切に伴う痛みのため患者に１日に何回も測定を強要するのは困難な状況にある。
【０００５】
そこで近年、このような欠点を有する上記侵襲型（侵入型）の測定に代えて、痛み等を伴わない非侵襲型（非侵入型）の測定方法が種々提案されている。
【０００６】
これらは主として、人の眼球部の角膜と水晶体との間にある前眼房を満たす眼房水中のグルコース濃度が、個人差はあるものの血中グルコース濃度と極めて高い相関関係を有していることに着目したものであり、眼房水中のグルコース濃度を非侵襲で測定するものである。
【０００７】
例えば特開昭５１−７５４９８号（米国特許第３，９５８，５６０ 号）によれば、眼房水中に入射した赤外線の旋光度を求めることにより、この旋光度と関連のあるグルコース濃度を得るものである。
【０００８】
また特表平６−５０３２４５ 号によればグルコースの誘導ラマン光を測定するものである。
【０００９】
さらに特開平６−２３７８９８ 号では、水晶体による反射光を光学的性質を測定するデバイスの発明が記載されている。その他、米国特許第５，４３３，１９７ 号においても眼房水中のグルコース濃度を測定する方法が記載されている。
【００１０】
【発明が解決しようとする課題】
しかしながら、上記特開平６−２３７８９８ 号に記載された発明は、角膜と眼房水との境界面での反射光を除去することができず、また角膜での吸収情報も一緒に検出されるため眼房水中のグルコース濃度決定の精度が低下する。
【００１１】
また、微小な吸光度変化を具体的にどのような技術的手段によって測定するかについての開示がないため実用上問題がある。
【００１２】
一方、米国特許第３，９５８，５６０ 号によれば、眼房水中のグルコース以外の多数の化合物が光学的に活性で偏光面の回転に関与する。また角膜は複屈折を示すため偏向面の回転を引き起こす。したがって旋光度により眼房水中のグルコース濃度を測定することは測定精度確保のうえで問題がある。
【００１３】
さらに特表平６−５０３２４５ 号によれば、グルコースの誘導ラマン光の測定のため大きな出力のポンプレーザ光を前眼房部に、視線光軸に対して垂直に導入するため実用的な測定系の構成が困難である。
【００１４】
本発明は上記事情に鑑みなされたものであって、眼房水中のグルコース濃度を非侵襲的に、かつ精度よく測定する装置および方法を提供することを目的とするものである。
【００１５】
また本発明はこのような非侵襲的に得られた眼房水中のグルコース濃度を用いて血液中のグルコース濃度を非侵襲的に測定する装置および方法を提供することをも目的とするものである。
【００１６】
【課題を解決するための手段】
本発明の第１のグルコース濃度測定方法は、
所定の光源から出射された低コヒーレンスな光を、互いに異なる２つの光路に沿って各別に進行する信号光と参照光とに分割し、
該信号光と参照光とで僅かな周波数差が生じるように両光のうち少なくとも一方を変調し、
所定の位置に予め配された眼球に前記信号光を照射し、
該眼球に照射された信号光の、角膜と前眼房との境界面による第１の後方散乱光と、前記参照光とを、該参照光の光路長を調整することにより干渉せしめ、
該干渉により得られた第１の干渉光の強度を測定し、
第１の干渉光の強度に基づいて前記第１の後方散乱光の強度を求め、
前記眼球に照射された信号光の、前眼房と水晶体との境界面による第２の後方散乱光と、前記参照光とを、該参照光の光路長を調整することにより干渉せしめ、
該干渉により得られた第２の干渉光の強度を測定し、
第２の干渉光の強度に基づいて前記第２の後方散乱光の強度を求め、
前記第１および第２の後方散乱光の強度に基づいて、前記前眼房内を満たす眼房水成分の光吸収特性を求め、
前記低コヒーレンスな光とは波長帯域の異なる他の複数の低コヒーレンスな光についても同様に前記光吸収特性を求め、
得られた複数の光吸収特性に基づいて前記眼房水成分中のグルコース濃度を求めることを特徴とするものである。
【００１７】
ここで上記低コヒーレンスな光としては、可干渉距離が数十μｍ程度と短い例えばＳＬＤ（Ｓｕｐｅｒ Ｌｕｍｉｎｅｓｃｅｎｔ Ｄｉｏｄｅ ）やＬＥＤ等が用いられる。なお実用上はより指向性の高いＳＬＤを用いるのが望ましい。
【００１８】
また、上記干渉光の強度を測定するとは、上記後方散乱光（信号光）と参照光との差周波数で強弱を繰り返すビート信号（干渉光）の強度を計測することを意味する。以下の発明においても同様である。
【００１９】
なお、上記本発明の第１のグルコース濃度測定方法において、上記各低コヒーレンスな光は、これら各低コヒーレンスな光の波長帯域よりも広い発光波長帯域の光の一部としてそれぞれ選択された光であってもよいし、互いに異なる複数の光源から各別に出射された光であってもよい。
【００２０】
また、複数の光源により構成した場合には、この波長帯域の異なる複数の光源から順次に低コヒーレンスな光を出射して、これらに対応した各干渉光を１つの光検出器を用いて検出すればよい。
【００２１】
本発明の第２のグルコース濃度測定方法は、
所定の光源から出射された、時間的に鋸歯状（例えば図６に示すよう）に周波数掃引されたコヒーレント光（例えばレーザー光）を互いに異なる２つの光路に沿って各別に進行する信号光と参照光とに分割し、
所定の位置に予め配された眼球に前記信号光を照射し、
該眼球の角膜と前眼房との境界面による前記信号光の第１の後方散乱光と、前記信号光（参照光と分割されてから眼球の角膜と前眼房との境界面に至るまでの信号光を意味する）および前記第１の後方散乱光（眼球の角膜と前眼房との境界面から参照光と干渉するまでのもの）と前記参照光（信号光と分割されてから第１の後方散乱光と干渉するまでのもの）との光路長差に基づく時間差をもって前記光源から出射した、前記第１の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを干渉せしめ、
該干渉により得られた第１の干渉光の強度を測定し、
該第１の干渉光の強度に基づいて前記第１の後方散乱光の強度を求め、
前記眼球の前眼房と水晶体との境界面による前記信号光の第２の後方散乱光と、前記信号光（参照光と分割されてから眼球の前眼房と水晶体との境界面に至るまでの信号光を意味する）および前記第２の後方散乱光（眼球の前眼房と水晶体との境界面から参照光と干渉するまでのもの）と前記参照光（信号光と分割されてから第２の後方散乱光と干渉するまでのもの）との光路長差に基づく時間差をもって前記光源から出射した、前記第２の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを干渉せしめ、
該干渉により得られた第２の干渉光の強度を測定し、
該第２の干渉光の強度に基づいて前記第２の後方散乱光の強度を求め、
前記第１および第２の後方散乱光の強度に基づいて、前記前眼房内を満たす眼房水成分の光吸収特性を求め、
前記コヒーレント光とは波長の異なる他の複数のコヒーレント光についても同様に前記光吸収特性を求め、
得られた複数の光吸収特性に基づいて前記眼房水成分中のグルコース濃度を求めることを特徴とするものである。
【００２２】
ここで上記第１の後方散乱光と干渉せしめられる参照光が「第１の後方散乱光とは周波数差を有する」のは、以下の理由による。
【００２３】
上記信号光が通過する光路長と第１の後方散乱光が通過する光路長との和と、参照光が通過する光路長とには差が設けられており（参照光の光路長の方が短い場合だけでなく、参照光の光路長の方を長く設定することも勿論可能である。）、この光路長差によって、例えば参照光の光路長の方が短い場合には、波面整合（干渉）せしめられる位置に到達するのは参照光の方が第１の後方散乱光よりも早い。
【００２４】
すなわち第１の後方散乱光がその位置に到達したときには、この第１の後方散乱光を生じせしめた信号光と分割された参照光は既にこの位置を通過しており、この信号光が光源から出射された時点よりも遅い時刻に出射されたコヒーレント光の一部である参照光が到達する。
【００２５】
この遅い時刻に出射されたコヒーレント光は、時間的に周波数掃引されているため、第１の後方散乱光とは僅かに周波数差を有するものとなっている。
【００２６】
なお、上記第２の後方散乱光と干渉せしめられる参照光は、第１の後方散乱光と干渉せしめられる参照光とも異なる周波数を有している。これは、第１の後方散乱光が角膜と前眼房との境界による散乱光であるのに対し、第２の後方散乱光は、その境界よりもさらに眼球の深部である前眼房と水晶体との境界による散乱光であるため、両者の光路長差（前眼房の厚さの２倍）分に応じた時間差が生じ、この結果、第２の後方散乱光と干渉せしめられる参照光は第１の後方散乱光と干渉せしめられる参照光よりも遅い時刻に光源から出射されたコヒーレント光によるものとなる。
【００２７】
なお、上記本発明の第２のグルコース濃度測定方法において、上記各コヒーレント光は、単一の光源からそれぞれ選択的に出射された光であってもよいし、または互いに異なる複数の光源から各別に出射された光であってもよい。
【００２８】
本発明の第３のグルコース濃度測定方法は、
所定の光源から出射された超短パルス光を眼球に照射し、
該眼球の角膜と前眼房との境界面による前記超短パルス光の第１の後方散乱光の強度および前眼房と水晶体との境界面による第２の後方散乱光の強度を各別に測定し、
該第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求め、
前記超短パルス光とは波長の異なる他の複数の超短パルス光についても同様に前記光吸収特性を求め、
該超短パルス光の波長ごとに得られた複数の光吸収特性に基づいて前記眼房水成分中のグルコース濃度を求めることを特徴とするものである。
【００２９】
ここで上記超短パルス光とは、少なくとも上記第１の後方散乱光の強度と、前眼房と水晶体との境界面による第２の後方散乱光の強度とを、時間的に分離して各別に測定し得る程度の非常に短時間（例えばフェムト秒〜ピコ秒単位程度）だけ発光するパルス状の光、例えばモードロックＴｉ：サファイアレーザーなど等を意味する。このような超短パルス光を用いることにより、第１の後方散乱光に対して前眼房を往復する距離に対応した時間だけ遅れた第２の後方散乱光を、例えばストリークカメラ等の時間分解可能の光検出器を用いて、第１の後方散乱光から分離して検出することができる。
【００３０】
なおこの方法においても、各超短パルス光は、単一の光源からそれぞれ選択的に出射された光であってもよいし、または、互いに異なる複数の光源から各別に出射された光であってもよい。
【００３１】
本発明の第４のグルコース濃度測定方法は、上述の各グルコース濃度測定方法のうちいずれかにより求められた眼房水成分中のグルコース濃度ごとに、対応する血中グルコース濃度を各別に侵襲的に求めて、その方法により求められた前記眼房水成分中のグルコース濃度と侵襲的に得た血中グルコース濃度との相関を予め求め、
この相関を求めた以後においては、上記各グルコース濃度測定方法のうちいずれかにより求めた眼房水成分中のグルコース濃度と前記相関とに基づいて、血中グルコース濃度を非侵襲的に求めることを特徴とするものである。
【００３２】
本発明の第１のグルコース濃度測定装置は、上記本発明の第１のグルコース濃度測定方法を実施するための装置であって、
互いに発光波長帯域の異なる複数の低コヒーレンスな光を出射する光源装置と、
該光源装置から出射された低コヒーレンスな光を互いに異なる２つの光路に沿って各別に進行する参照光と眼球に入射される信号光とに分割する光路分割手段と、
該信号光と参照光とで僅かな周波数差が生じるように両光のうち少なくとも一方を変調する、該少なくとも一方の光路上に設けられた変調手段と、
前記参照光が進行する光路の長さを調整する光路長調整手段と、
前記眼球の角膜と前眼房との境界面による信号光の第１の後方散乱と前記参照光、および、前記眼球の前眼房と水晶体との境界面による前記信号光の第２の後方散乱光と前記参照光とを、それぞれ波面整合させる波面整合手段と、
前記参照光と前記第１の後方散乱光との波面整合による第１の干渉光および前記参照光と前記第２の後方散乱光との波面整合による第２の干渉光の強度を光電的に検出する光検出器と、
該各干渉光の強度に基づいて前記第１および第２の後方散乱光の強度を求めるヘテロダイン演算手段と、
前記第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段と、
前記複数の低コヒーレンスな光ごとに得られた複数の光吸収特性に基づいて、前記眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段とを備えたことを特徴とするものである。
【００３３】
ここで上記光源装置は、前記各コヒーレンスな光の波長帯域よりも広い発光波長帯域の低コヒーレンスな光を出射する単一の光源と、
該発光波長域の広い低コヒーレンスな光から、前記各低コヒーレンスな光を波長選択する波長選択手段とにより構成されているものであってもよいし、前記各低コヒーレンスな光を各別に出射する複数の光源により構成されているものであってもよい。
【００３４】
本発明の第２のグルコース濃度測定装置は、上記本発明の第２のグルコース濃度測定方法を実施するための装置であって、
図６に示すように、時間的に鋸歯状に周波数掃引された、互いに波長の異なる複数のコヒーレント光を出射する光源装置と、
該光源装置から出射された周波数掃引されたコヒーレント光を互いに異なる２つの光路に沿って各別に進行する参照光と眼球に入射される信号光とに分割する光路分割手段と、
前記眼球の角膜と前眼房との境界面による前記信号光の第１の後方散乱光と、前記信号光および前記第１の後方散乱光と前記参照光との光路長差に基づく時間差をもって前記光源から出射した、前記第１の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを、および、前記眼球の前眼房と水晶体との境界面による前記信号光の第２の後方散乱光と、前記信号光および前記第２の後方散乱光と前記参照光との光路長差に基づく時間差をもって前記光源から出射した、前記第２の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを、それぞれ波面整合させる波面整合手段と、
前記第１の後方散乱光と該第１の後方散乱光に対して僅かな周波数差を有する参照光との波面整合により得られた第１の干渉光、および前記第２の後方散乱光と該第２の後方散乱光に対して僅かな周波数差を有する参照光との波面整合により得られた第２の干渉光の強度を光電的に検出する光検出器と、
該各干渉光の強度に基づいて前記第１および第２の後方散乱光の強度を求めるヘテロダイン演算手段と、
前記第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段と、
前記複数のコヒーレント光ごとに得られた複数の光吸収特性に基づいて、前記眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段とを備えたことを特徴とするものである。
【００３５】
このグルコース濃度測定装置においても、光源装置は、前記複数のコヒーレント光のうちいずれか１つを選択的に出射し得る単一の光源と、該光源に対して、前記複数のコヒーレント光のうち選択的にいずれか１つを出射させる制御を施す制御手段とにより構成されているものであってもよいし、または、前記各コヒーレント光を各別に出射する複数の光源により構成されているものであってもよい。
【００３６】
本発明の第３のグルコース濃度測定装置は、上記本発明の第３のグルコース濃度測定方法を実施するための装置であって、
互いに波長の異なる複数の超短パルス光を出射する光源装置と、
該超短パルス光を眼球に入射せしめ、該眼球の角膜と前眼房との境界面による該超短パルス光の第１の後方散乱光の強度、および前眼房と水晶体との境界面による該超短パルス光の第２の後方散乱光の強度を時系列的に各別に求める光時間領域後方散乱測定手段と、
前記第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段と、
前記複数の超短パルス光ごとに得られた複数の光吸収特性に基づいて、前記眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段とを備えたことを特徴とするものである。
【００３７】
また上記光源装置は、前記複数の超短パルス光のうちいずれか１つを選択的に出射し得る単一の光源と、該光源に対して、前記複数の超短パルス光のうち選択的にいずれか１つを出射させる制御を施す制御手段とにより構成されているものであってもよいし、または、前記各超短パルス光を各別に出射する複数の光源により構成されているものであってもよい。
【００３８】
本発明の第４のグルコース濃度測定装置は、上記各グルコース濃度測定装置のうちいずれかにより求められた眼房水成分中のグルコース濃度と予め例えば侵襲的に得た血中グルコース濃度との相関を示すテーブルを備え、
上記各グルコース濃度測定装置のうちいずれかにより求められた眼房水成分中のグルコース濃度および前記テーブルに基づいて、前記血中グルコース濃度を非侵襲的に求めることを特徴とするものである。
【００３９】
【発明の効果】
本発明の第１および第２のグルコース濃度測定方法並びに本発明の第１および第２のグルコース濃度測定装置によれば、光ヘテロダイン後方散乱測定法を用いることにより、まず角膜と前眼房の境界面による入射光の微弱な第１の後方散乱光ＩＲ２および前眼房と水晶体の境界面による入射光の微弱な第２の後方散乱光ＩＲ３を各別に精度よく検出することができる。
【００４０】
ここで図２に示すように、眼球部への入射光の強度をＩ０ 、空気／角膜間の反射率をＲ１ 、角膜／前眼房間の反射率をＲ２ 、前眼房／水晶体間の反射率をＲ３ 、角膜による入射光の片光路の光学的吸収率をα１ 、前眼房（眼房水）による入射光の片光路の光学的吸収率をα２ とすると、第１の後方散乱光の強度ＩＲ２、第２の後方散乱光の強度ＩＲ３は下記式（１），（２）のように表すことができる。
【００４１】
この精度よく検出された各後方散乱光ＩＲ２，ＩＲ３の比（ＩＲ３／ＩＲ２）を計算すると、
となる。ＩＲ３，ＩＲ２，Ｒ３ ，Ｒ２ は既知であるから眼房水による入射光の片光路の光学的吸収率α２ を求めることができる。
【００４２】
ここで眼房水にはグルコース以外にもＮａＣｌをはじめとして複数の成分が含まれているため、上記の作用を複数の波長の入射光について繰り返し行ない、多変量解析を含む公知な近赤外分光分析法を適用することにより、眼房水中のグルコースだけの濃度を特定することができる。
【００４３】
このようにして非侵襲的に得られた眼房水中のグルコース濃度と従来の侵襲的に得られた血中グルコース濃度との相関関係を例えば変換テーブル等として各患者ごとに予め準備しておくことにより、その後は血中グルコース濃度を測定することなく眼房水中のグルコース濃度を測定すればよく、患者に苦痛を与えずに繰り返し精度のよい測定を行なうことが可能となる。
【００４４】
本発明の第３のグルコース濃度測定方法および本発明の第３のグルコース濃度測定装置によれば、第１の後方散乱光の強度および第２の後方散乱光の強度を図８に示すように時間的に分離して各別に測定することができ、この図８に示した波形を時間について積分することによりその光量、すなわち上記式（１）に示すＩＲ２および式（２）に示すＩＲ３を求めることができる。
【００４５】
したがって式（３）により、眼房水による入射光の片光路の光学的吸収率α２ を求めることができる。
【００４６】
そして本願の上記発明の場合と同様に、複数の波長の入射光について繰り返し行ない、多変量解析を含む公知な近赤外分光分析法により眼房水中のグルコース濃度を特定することができる。
【００４７】
【発明の実施の形態】
以下、本発明の具体的な実施の形態について図面を用いて説明する。
【００４８】
図１は本発明の第１のグルコース濃度測定方法を実施するための基本的な装置構成を示す図である。
【００４９】
図示のグルコース濃度測定装置は、複数の波長帯域の、可干渉距離が数十μｍ程度の低コヒーレンスな光を出射する光源装置１０と、この光源装置１０から出射された低コヒーレンスな光を眼球１００ に照射せしめ、眼球１００ の角膜１１０ と前眼房１２０ との境界面による入射した光の第１の後方散乱光の強度、および前眼房１２０ と水晶体１３０ との境界面による入射した光の第２の後方散乱光の強度、を光ヘテロダイン後方散乱測定法により求める光ヘテロダイン測定手段２０と、これら第１および第２の後方散乱光の強度に基づいて、前眼房１２０ を満たす眼房水成分の入射光の波長に対する光吸収率を求める光吸収特性分析手段３０と、光源装置１０から出射する光の波長ごとに得られた複数の光吸収率に基づいて、多変量解析を含む公知の近赤外分光分析法により眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段４０を備えた構成である。
【００５０】
ここで図１に示した基本的な構成をより具体化した構成を図３に示す。
【００５１】
図示のグルコース濃度測定装置では、図１に示した光源装置１０として、中心波長がλ１であって波長帯域の広い低コヒーレンスな光を出射するＳＬＤ１１ａ ，中心波長がλ２であって波長帯域の広い低コヒーレンスな光を出射するＳＬＤ１１ｂ ，…，中心波長がλｎであって波長帯域の広い低コヒーレンスな光を出射するＳＬＤ１１ｘ という複数の光源を備え、
光ヘテロダイン測定手段２０として、光源装置１０から出射された低コヒーレンスな光を互いに異なる２つの光路に沿って各別に進行する信号光と参照光とに分割し、また２つの光を波面整合する４つのビームスプリッタ２１ａ ，２１ｂ ，２１ｃ ，２１ｄ と、信号光と参照光とで僅かな周波数差Δｆが生じるように両光を各別に変調する光音響変調器２５ａ ，２５ｂ と、参照光が進行する光路の長さを調整する、光軸方向に移動可能の反射ミラー２７と、眼球１００ に導光された信号光の後方散乱光と参照光とが波面整合して得られた干渉光の強度を測定する光検出器２８と、この干渉光の強度に基づいて光ヘテロダイン後方散乱測定法により後方散乱光の強度を求める光ヘテロダイン演算手段２９とを備えた構成である。
【００５２】
ここで各ＳＬＤ１１ａ ，１１ｂ ，…，１１ｘ は順次に光を出射するものである。なお光音響変調器は必ずしも本実施形態のように、両光（信号光および参照光）をそれぞれ変調するものである必要はなく、少なくとも一方を変調して両光に僅かな周波数差を生じせしめればよい。
【００５３】
次に本実施形態のグルコース濃度測定装置の作用について説明する。
【００５４】
まずＳＬＤ１１ａ から中心波長λ１の低コヒーレンスな光が出射される。このとき他の光源は光を出射しない。出射された光は光路上の最初のビームスプリッタ２１ａ により、信号光と参照光とに分割される。これら信号光および参照光は各光路上にそれぞれ設けられた光音響変調器２５ａ ，２５ｂ により変調され、両光には僅かな周波数差Δｆが生じる。
【００５５】
次に、変調された両光のうち信号光はビームスプリッタ２１ｂ を透過して眼球１００ に入射する。
【００５６】
眼球１００ に入射した光は大きく４つの光に分割される。すなわち、
▲１▼空気／角膜間の境界面で反射する光、
▲２▼角膜／前眼房間の境界面で反射する光、
▲３▼前眼房／水晶体間の境界面で反射する光、
▲４▼水晶体を透過する光、である。
【００５７】
これらのうち後方散乱光として、▲１▼〜▲３▼の光が入射光である信号光とは反対の方向に眼球１００ から出射する。
【００５８】
眼球１００ から出射したこれらの各後方散乱光はビームスプリッタ２１ｂ により反射され次のビームスプリッタ２１ｄ を透過する。
【００５９】
一方、変調された後の参照光は、ビームスプリッタ２１ｃ により反射されてミラー２７に入射し、このミラー２７で反射されて再度ビームスプリッタ２１ｃ に入射し、このビームスプリッタ２１ｃ を透過し、次のビームスプリッタ２１ｄ により反射される。
【００６０】
このとき、ミラー２７をその光軸方向（矢印Ｘ方向）に移動すると、参照光がビームスプリッタ２１ｄ に到達するまでに通過する光路長が変化する。
【００６１】
また、参照光、眼球に入射する信号光、および眼球で反射された後の後方散乱光は、可干渉距離の短い低コヒーレンスな光であるため、前述の３つの後方散乱光のうち、信号光として通過した距離および後方散乱光としてビームスプリッタ２１ｄ に到達するまでの通過距離の和が、参照光のビームスプリッタ２１ｄ に到達するまでの通過距離に略等しい後方散乱光だけが参照光と干渉し、この干渉する両光の差周波数（Δｆ）で強弱を繰り返すビート信号が発生する。
【００６２】
すなわち、
（１）▲１▼の後方散乱光と参照光とが干渉するときは、
ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｂ 間の距離とビームスプリッタ２１ｂ 〜空気／角膜間の境界面間の往復距離とビームスプリッタ２１ｂ 〜ビームスプリッタ２１ｄ 間の距離の和が、ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｃ 間の距離とビームスプリッタ２１ｃ 〜ミラー２７間の境界面間の往復距離とビームスプリッタ２１ｃ 〜ビームスプリッタ２１ｄ 間の距離の和とが略等しく、
（２）▲２▼の後方散乱光が参照光と干渉するときは、
ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｂ 間の距離とビームスプリッタ２１ｂ 〜角膜／前眼房間の境界面間の往復距離とビームスプリッタ２１ｂ 〜ビームスプリッタ２１ｄ 間の距離の和が、ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｃ 間の距離とビームスプリッタ２１ｃ 〜ミラー２７間の境界面間の往復距離とビームスプリッタ２１ｃ 〜ビームスプリッタ２１ｄ 間の距離の和とが略等しく、
（３）▲３▼の後方散乱光が参照光と干渉するときは、
ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｂ 間の距離とビームスプリッタ２１ｂ 〜前眼房／水晶体間の境界面間の往復距離とビームスプリッタ２１ｂ 〜ビームスプリッタ２１ｄ 間の距離の和が、ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｃ 間の距離とビームスプリッタ２１ｃ 〜ミラー２７間の境界面間の往復距離とビームスプリッタ２１ｃ 〜ビームスプリッタ２１ｄ 間の距離の和とが略等しい。
【００６３】
したがって、まずビームスプリッタ２１ｃ に近い位置に配置したミラー２７を、信号光の空気／角膜１１０ 間の境界面による後方散乱光と参照光との干渉によるビート信号が光検出器２８で認められる位置まで、序々にビームスプリッタ２１ｃ から遠ざける方向に移動する。
【００６４】
ここでさらにミラー２７をビームスプリッタ２１ｃ から遠ざける方向に移動するとビート信号は発生しなくなるが、ミラー２７の移動を続けると再度ビート信号が光検出器２８により検出される。この２度目のビート信号は、信号光の角膜１１０ ／前眼房１２０ 間の境界面による後方散乱光（以下、第１の後方散乱光という）と参照光とが干渉することにより発生したものであり、最初のビート信号の検出から角膜１１０ の厚さ分だけミラー２７を移動した時に発生する。
【００６５】
この第１の後方散乱光と参照光とが干渉することにより発生したビート信号を、第１の干渉光としてその強度ＢＲ２を光検出器２８により検出する。
【００６６】
さらにミラー２７をビームスプリッタ２１ｃ から遠ざける方向に移動すると、またビート信号は発生しなくなるが、ミラー２７の移動を続けると３度目のビート信号が光検出器２８により検出される。この３度目のビート信号は、信号光の前眼房１２０ ／水晶体１３０ 間の境界面による後方散乱光（以下、第２の後方散乱光という）と参照光とが干渉することにより発生したものであり、２度目のビート信号の検出から前眼房１２０ の厚さ分だけミラー２７を移動した時に発生する。
【００６７】
この第２の後方散乱光と参照光とが干渉することにより発生したビート信号を第２の干渉光としてその強度ＢＲ３を光検出器２８により検出する。
【００６８】
このようにして得られた上記２つのビート信号の強度ＢＲ２，ＢＲ３は電気信号に変換されてヘテロダイン演算手段２９に入力され、ヘテロダイン演算手段２９はこれらの電気信号に基づいて第１の後方散乱光の強度ＩＲ２および第２の後方散乱光の強度ＩＲ３を算出する。
【００６９】
すなわち各ビート信号の強度ＢＲ２，ＢＲ３はそれぞれ第１の後方散乱光の強度ＩＲ２と参照光の強度との積の平方根、第２の後方散乱光の強度ＩＲ３と参照光の強度との積の平方根であり、参照光の強度については予め測定した値を用いればよく、この結果、各後方散乱光の強度ＩＲ２およびＩＲ３を求めることができる。
【００７０】
ここで図２を参照して、光吸収特性分析手段３０による作用を説明する。
【００７１】
まず眼球１００ への入射光である信号光の強度をＩ０ 、空気／角膜１１０ 間の反射率をＲ１ 、角膜１１０ ／前眼房１２０ 間の反射率をＲ２ 、前眼房１２０ ／水晶体１３０ 間の反射率をＲ３ 、角膜１１０ による入射光の片光路の光学的吸収率をα１ 、前眼房（眼房水）１２０ による入射光の片光路の光学的吸収率をα２ とすると、第１の後方散乱光の強度ＩＲ２、第２の後方散乱光の強度ＩＲ３は下記式（１），（２）のように表すことができる。
【００７２】
図２および式（２）からも解されるように、第２の後方散乱光の強度ＩＲ３は、眼房水１２０ による吸収率α２ と角膜１１０ による吸収率α１ との双方に依存している。
【００７３】
ここで、光吸収特性分析手段３０は各後方散乱光の強度ＩＲ２，ＩＲ３の比（ＩＲ３／ＩＲ２）を下記式（３）にしたがって計算する。
【００７４】
ＩＲ３，ＩＲ２は上記ヘテロダイン演算手段２９により精度よく求められたものであり、Ｒ３ ，Ｒ２ は既知であるから、眼房水のみによる光学的吸収率α２ が求められる。
【００７５】
ここで眼房水にはグルコース以外にもＮａＣｌをはじめとして複数の成分が含まれているため、上記１つの波長λ１の入射光によって得られた眼房水１２０ の光学的吸収率α２ （この波長λ１の入射光による眼房水１２０ の光学的吸収率α２ を以下、α２ （１） と表記する）のみによっては、眼房水中のグルコース濃度を特定することができない。
【００７６】
そこで次に、２番目のＳＬＤ１１ｂ から中心波長λ２の低コヒーレンスな光が出射される。このとき最初のＳＬＤ１１ａ を含めて他のＳＬＤは光を出射しない。
【００７７】
ＳＬＤ１１ｂ から出射された中心波長λ２の低コヒーレンスな光も、最初の中心波長λ１の光と同様に光路上の最初のビームスプリッタ２１ａ により信号光と参照光とに分割され、信号光は眼球１００ に入射し、参照光はミラー２７に反射され、ビームスプリッタ２１ｄ において得られた干渉光に基づいて、光検出器２８、ヘテロダイン演算手段２９および光吸収特性分析手段３０による前述と同様の作用により、第２の中心波長λ２の入射光による眼房水１２０ の光学的吸収率α２ （２） が求められる。
【００７８】
以下、各ＳＬＤ１１ｃ ，…，ＳＬＤ１１ｘ についても同様の操作を繰り返すことにより、複数の波長λ１，λ２，…，λｎの各入射光による眼房水１２０ の光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） を得ることができる。
【００７９】
このようにして得られた入射光の波長ごとの光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） はグルコース濃度算出手段４０に入力され、グルコース濃度算出手段４０はこれらの複数の光学的吸収率に基づいて、多変量解析を含む公知な近赤外分光分析法により、複数の成分を含む眼房水１２０ からグルコースの濃度を特定する。
【００８０】
このように本実施形態のグルコース濃度測定装置によれば、眼房水中のグルコース濃度を非侵襲的に精度よく測定することができる。
【００８１】
そしてこのようにして得られた眼房水中のグルコース濃度は、予め各患者ごとに上記装置で求めた眼房水中のグルコース濃度と従来の侵襲的に得た血中グルコース濃度との相関関係を例えば変換テーブル等として、眼房水／血中グルコース濃度変換手段５０のメモリ５２に記憶せしめておき、上記グルコース濃度測定装置からある患者の眼房水中グルコース濃度が変換手段５０に入力されたときに、眼房水／血中グルコース濃度変換手段５０の変換演算手段５１がこのメモリ５２に記憶せしめられた変換テーブルを参照して、その患者の眼房水中グルコース濃度に応じた血中グルコース濃度を出力するようにすれば、テーブル作成後においては、眼房水中のグルコース濃度を上記装置により非侵襲的に計測してその値を求め、得られた眼房水中グルコース濃度とテーブルとに基づいて血中グルコース濃度を算出することができ、血中グルコース濃度の非侵襲的測定を実現することができる。
【００８２】
後述する以下の実施形態においても眼房水／血中グルコース濃度変換手段５０を設ければ上記と同様に効果を得ることができる。
【００８３】
なお、本実施形態のグルコース濃度測定装置は、光源装置１０として、出射する光の波長帯域が互いに異なる複数の光源１１ａ ，１１ｂ ，…，１１ｘ を備えたものについて説明したが、本発明のグルコース濃度測定方法，装置は必ずしもこのように複数の光源を備えた構成である必要はなく、例えば図４（Ａ）に示すように、上述の各低コヒーレンスな光の波長帯域の全体を含むような発光波長域のより広い光を出射する１つの光源１１と、この光源１１から出射された発光波長域の広い光から上述の各低コヒーレンスな光（所定の波長帯域のもの）を選択する波長選択手段１２とからなる構成としてもよい。
【００８４】
また本実施形態では、非侵襲的に得られた眼房水中のグルコース濃度を用いて非侵襲的に血液中のグルコース濃度を得るものであり、本発明の第３のグルコース濃度測定方法，装置の実施例形態でもある。
【００８５】
図５は、本発明の第２のグルコース濃度測定方法を実施するための具体的な一形態を示す図である。
【００８６】
図示のグルコース濃度測定装置の基本的な構成は図１に示したものと同様であるが、図１に示した構成の光源装置１０を、波長λ１のレーザー光を出射する光源１１ａ ′，波長λ２のレーザー光を出射する光源１１ｂ ′，…，波長λｎのレーザー光を出射する光源１１ｘ ′と、これらのレーザー光の入力を受けて図６に示すように各光の周波数を一定の帯域で時間的に鋸歯状に掃引する周波数掃引手段１３とを備えた光源装置１０′としたものである。なお、本実施形態においてはコヒーレント光の一形態としてレーザー光を適用したものについて説明するが、本発明の第２のグルコース濃度測定方法（装置）においてはコヒーレント光であれば、特にレーザー光に限るものではない。
【００８７】
さらに図１に示した構成の光ヘテロダイン測定手段２０を、ビームスプリッタ２３、光路分割手段としてのおよび波面整合手段としての反射板２６と、光検出器２８と、干渉光の強度に基づいて後方散乱光の強度を求めるヘテロダイン演算手段２９′とを備えた構成としたものである。
【００８８】
周波数掃引手段１３により周波数掃引されて光源装置１０′から出射したレーザー光はビームスプリッタ２３を透過して眼球部１００ に入射せしめられ、一方、反射板２６は、眼球１００ に照射されたレーザー光の各反射面で反射された後方散乱光と、この後方散乱光よりも時間的に遅れて光源装置１０′から出射されたレーザー光であって後方散乱光に対して僅かな周波数差を有しその反射面２６ａ において反射した光（参照光）とを波面整合して干渉光を生じせしめる作用を果たすものである。
【００８９】
またヘテロダイン演算手段２９′は、反射板２６により波面整合された後方散乱光とこの後方散乱光より遅れて光源装置１０′から出射され、後方散乱光と同時に反射板２６に到達して反射したレーザー光（参照光）との干渉光の、光検出器２８により検出された強度に基づいて、光ヘテロダイン後方散乱測定法により後方散乱光の強度を求める作用をなすものである。
【００９０】
なお、本実施形態のグルコース濃度測定装置においても、図１および図３に示した光吸収特性分析手段３０およびグルコース濃度算出手段４０を備えているが、その作用は図１または図３に示したものと同様であるため図示を省略している。
【００９１】
次に本実施形態のグルコース濃度測定装置の作用について説明する。
【００９２】
まず光源１１ａ ′から波長λ１のレーザー光が出射される。この光は光源装置１０′の周波数掃引手段１３によりその周波数が掃引されて光源装置１０′から出射される。ここで光源装置１０′から時刻ｔ０ において出射された光の周波数をｆ０ とする（図６参照）。
【００９３】
光源装置１０′から出射された周波数ｆ０ の光は、その光路上に配されたビームスプリッタ２３を透過して反射板２６に入射する。この反射板２６はその反射面２６ａ が半透鏡に形成されているため、入射した光の一部を反射し、残りの透過を許容する。この反射した光は参照光として、透過した光は信号光として作用する。
【００９４】
したがって、周波数ｆ０ の光もその大半が反射面２６ａ を透過し、一部が反射面２６ａ で反射される。
【００９５】
ここでまず反射面２６ａ を透過した周波数ｆ０ の光は、眼球１００ に入射し、前述と同様に、
▲１▼空気／角膜間の境界面で反射する光、
▲２▼角膜／前眼房間の境界面で反射する光、
▲３▼前眼房／水晶体間の境界面で反射する光、
▲４▼水晶体を透過する光、に分割される。
【００９６】
そしてこれらのうち後方散乱光として、▲１▼〜▲３▼の光が入射方向とは反対の方向に眼球１００ から出射し、反射板２６に到達する。
【００９７】
ここで、反射面２６ａ を透過した周波数ｆ０ の信号光が眼球１００ に入射し、後方散乱光として再度、反射板２６の反射面２６ａ に到達するまでの光路長について考察すると、
（１）上記▲１▼の後方散乱光の場合は、反射板２６の反射面２６ａ 〜空気／角膜１１０ 間の境界面までの往復距離、
（２）上記▲２▼の後方散乱光の場合は、反射板２６の反射面２６ａ 〜空気／角膜１１０ 間の境界面までの距離と、空気／角膜１１０ 間の境界面〜角膜１１０ ／前眼房１２０ 間の境界面までの距離との和の往復分、
（３）上記▲３▼の後方散乱光の場合は、反射板２６の反射面２６ａ 〜空気／角膜１１０ 間の境界面までの距離と、空気／角膜１１０ 間の境界面〜角膜１１０ ／前眼房１２０ 間の境界面までの距離と、角膜１１０ ／前眼房１２０ 間の境界面〜前眼房１２０ ／水晶体１３０ 間の境界面までの距離との和の往復分、
となる。
【００９８】
したがって、これらの各後方散乱光の光路長（眼球１００ に入射するまでの光路長を含む）は互いに異なるものであり、この結果、各後方散乱光が反射板２６の反射面２６ａ に到達するには、この光路長差に応じた時間差が生じる。
【００９９】
一方、光源装置１０′からは時間的に周波数掃引されたレーザー光が連続的に出射されており、これらの光も連続的に反射板２６に到達する。
【０１００】
この結果、例えば上記▲１▼の後方散乱光（周波数ｆ０ ）は、反射板２６の反射面２６ａ において、時刻ｔ１ （図６参照）に光源装置１０′から出射された周波数ｆ１ の光と出会って波面整合される。したがって、この波面整合による干渉光のビート信号は両光の差周波数（ｆ１ −ｆ０ ）で強弱を繰り返す。
【０１０１】
また、上記▲２▼の後方散乱光（周波数ｆ０ ）は、反射板２６の反射面２６ａ において、時刻ｔ２ に光源装置１０′から出射された周波数ｆ２ の光と出会って波面整合される。したがって、この波面整合による干渉光のビート信号の周波数は両光の差周波数である（ｆ２ −ｆ０ ）となる。
【０１０２】
さらに、上記▲３▼の後方散乱光（周波数ｆ０ ）は、反射板２６の反射面２６ａ において、時刻ｔ３ に光源装置１０′から出射された周波数ｆ３ の光と出会って波面整合される。したがって、この波面整合による干渉光のビート信号の周波数は両光の差周波数である（ｆ３ −ｆ０ ）となる。
【０１０３】
光源装置１０′からは周波数が刻々と変化するレーザー光が連続的に出射されるが、上記各光路長差は不変であるため、上記各後方散乱光との干渉によるビート信号の周波数は常に一定に保たれる。
【０１０４】
そしてこのように反射板２６の反射面２６ａ において各後方散乱光ごとにそれぞれ干渉した光はビームスプリッタ２３により反射されて光検出器２８に入射する。
【０１０５】
光検出器２８は、これら強弱の周波数が互いに異なる３つのビート信号が合成された光の強度を検出し、この強度を電気信号に変換してヘテロダイン演算手段２９′に入力する。
【０１０６】
ヘテロダイン演算手段２９′は、入力された電気信号に基づいて、各後方散乱光ごとのビート信号を分離し、この分離した各ビート信号の強度に基づいて各後方散乱光の強度を求める。
【０１０７】
なお、各後方散乱光ごとのビート信号の分離は、角膜１１０ および前眼房１２０ の光軸方向の厚さが既知であれば、その厚さに応じた光の通過時間の差に基づいて周波数の掃引速度の関係から、各後方散乱光によるビート信号の周波数差を容易に求めることができ、または角膜１１０ および前眼房１２０ の厚さが未知であっても、得られた合成ビート信号について公知の周波数分析を行なうことにより、各後方散乱光によるビート信号の周波数を容易に求めることができる。
【０１０８】
このようにヘテロダイン演算手段２９′により求められた各後方散乱光の強度は光吸収特性分析手段３０に入力される。
【０１０９】
そして、光源装置１０′の２番目以下の光源１１ｂ ′，…，１１ｘ ′についても上述と同様の操作を繰り返すことにより、複数の波長の各入射光による眼房水１２０ の光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） を得ることができ、得られた入射光の波長ごとの光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） はグルコース濃度算出手段４０に入力され、グルコース濃度算出手段４０はこれらの複数の光学的吸収率に基づいて、多変量解析を含む公知な近赤外分光分析法により、複数の成分を含む眼房水１２０ からグルコースの濃度を特定する。
【０１１０】
このように本実施形態のグルコース濃度測定装置によれば、眼房水中のグルコース濃度を非侵襲的に精度よく測定することができる。
【０１１１】
なお、本実施形態のグルコース濃度測定装置においても光源装置１０′として、各コヒーレント光を各別に出射する複数の光源１１ａ ′，１１ｂ ′，…，１１ｘ ′からなる構成について説明したが、必ずしもこのように各別の光源であるものに限らず、図４（Ｂ）に示すように、これらの各コヒーレント光を選択的に出射し得る単一の光源１１′と、光源１１′に対して、上記複数のコヒーレント光のうち選択的にいずれか１つを出射させる制御を施す制御手段１４とからなるものであってもよい。
【０１１２】
図７は本発明の第３のグルコース濃度測定方法を実施するための装置構成を示す図である。
【０１１３】
図示のグルコース濃度測定装置は、波長λ１の超短パルス光を出射するＴｉ：サファイアレーザー光源（以下、単に光源という）１１ａ ″，波長λ２の超短パルス光を出射する光源１１ｂ ″，…，波長λｎの超短パルス光を出射する１１ｘ ″という複数の超短パルス光を出射する光源装置１０″と、これらの超短パルス光を眼球１１０ に入射せしめ、眼球１００ の角膜１１０ と前眼房１２０ との境界面による超短パルス光の第１の後方散乱光の強度および前眼房１２０ と水晶体１３０ との境界面による超短パルス光の第２の後方散乱光の強度を時系列的に各別に求める、例えばストリークカメラ等の光検出器を含む光時間領域後方散乱測定手段６０と、第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段３０と、複数の超短パルス光ごとに得られた複数の光吸収特性に基づいて、眼房水１２０ 成分中のグルコース濃度を求めるグルコース濃度算出手段４０とを備えた構成である。
【０１１４】
ここで上記超短パルス光は、例えばフェムト秒〜ピコ秒という極めて短い時間だけ発光するパルス状の光を意味する。
【０１１５】
また光時間領域後方散乱測定手段６０は具体的には、光源装置１０″から出射された超短パルス光を、眼球１００ に入射する信号光と入射光による眼球１００ の各後方散乱光を検出するタイミングをとるための参照光とに分割するビームスプリッタ６１と、この参照光を導光するミラー６２と、入射光による眼球１００ の各後方散乱光を図８に示すように時間的に分離して各別に測定することができるストリークカメラ等の光検出器６３と、このように時間的に分離して検出された各後方散乱光の強度を時間について積分することによりその光量を求める眼房水透過率演算手段６４とを備えた構成である。
【０１１６】
なお参照光は、上述した後方散乱光の検出のタイミングをとるだけでなく、眼球１００ へ入射する信号光の光量をモニタするという意義をも有して光検出器６３に入力される。
【０１１７】
次に本実施形態のグルコース濃度測定装置の作用について説明する。
【０１１８】
まず光源１１ａ ″から波長λ１の超短パルス光が出射される。この光は光時間領域後方散乱測定手段６０に入力され、まずビームスプリッタ６１により信号光と参照光とに分割される。
【０１１９】
このうち信号光は眼球１００ に入射する。眼球１００ に入射した信号光は既に各実施形態の説明において既述したように、４つの光に大きく分割される。そしてそれら４つの光のうち、
▲１▼空気／角膜間の境界面で反射する光、
▲２▼角膜／前眼房間の境界面で反射する光、
▲３▼前眼房／水晶体間の境界面で反射する光、がそれぞれ後方散乱光として眼球１００ から出射する。
【０１２０】
このとき、各後方散乱光は信号光の入射方向についての反射面の深さが互いに相違するため、眼球１００ から出射するタイミングが異なり、上記▲１▼の後方散乱光、▲２▼の後方散乱光、▲３▼の後方散乱光、の順に眼球１００ から出射する。そして各後方散乱光は入射した信号光による反射光であるため、いずれも超短パルス光であり、この超短パルス光の発光時間は上記反射面の深さの差を通過する時間差に対して短いため、各後方散乱光は完全に分離されて眼球１００ から順次に出射する。
【０１２１】
このように順次に眼球１００ から出射された各後方散乱光は、ビームスプリッタ６１で反射されて順次に時間分解可能の光検出器６３に入力する。
【０１２２】
この光検出器６３には、各後方散乱光が入射するのに先だって、ビームスプリッタ６１で反射され、各ミラー６２により導光された参照光が入力される。この参照光が光検出器６３による検出タイミングのトリガーとして作用し、この参照光自体の強度が検出されるとともに、各後方散乱光の強度も各別に検出される。
【０１２３】
光検出器６３により時間分解されて時系列的に各別に検出された参照光の強度および各後方散乱光の強度の様子を図８に示す。
【０１２４】
光検出器６３はこのように検出した各光の強度を電気信号に変換して眼房水透過率演算手段６４に入力する。
【０１２５】
眼房水透過率演算手段６４は入力された各光の強度に対応した電気信号を時間について積分することにより各光の光量を求める。求められた各光の光量は光吸収特性分析手段３０に入力され、光吸収特性分析手段３０はこの各光の光量に応じて前記各実施形態の場合と同様の作用により、眼房水１２０ の光学的吸収率を求める。
【０１２６】
以上の作用と同様の操作を、上記光源装置１０″から出射する超短パルス光の波長を変えて、すなわち互いに異なる単一波長のパルス光を出射する複数の光源１１ａ ″，１１ｂ ″，…，１１ｘ ″を順次に代えることによって行ない、各入射光の波長に応じた眼房水１２０ の複数の光学的吸収率を求める。
【０１２７】
以上の操作により得られた入射光の波長ごとの光学的吸収率は、グルコース濃度算出手段４０に入力され、グルコース濃度算出手段４０はこれらの複数の光学的吸収率に基づいて、多変量解析を含む公知な近赤外分光分析法により、複数の成分を含む眼房水１２０ からグルコースの濃度を特定する。
【０１２８】
このように本実施形態のグルコース濃度測定装置によれば、眼房水中のグルコース濃度を非侵襲的に精度よく測定することができる。
【０１２９】
なお、本実施形態のグルコース濃度測定装置においても、光源装置１０″として、出射する超短パルス光の波長が互いに異なる複数の光源１１ａ ″，１１ｂ ″，…，１１ｘ ″からなる構成について説明したが、必ずしも複数の光源を備えた構成である必要はなく、前述の第２の実施形態の光源装置１０′の場合と同様に、これらの各超短パルス光を選択的に出射し得る単一の光源と、この光源に対して、上記複数の超短パルス光のうち選択的にいずれか１つを出射させる制御を施す制御手段とからなるものであってもよい。
【図面の簡単な説明】
【図１】本発明の第１のグルコース濃度測定方法を実施するための基本的な装置構成を示す図
【図２】眼球への入射光と後方散乱光との関係を示す図
【図３】本発明の第１のグルコース濃度測定方法を実施するためのより具体的な装置構成を示す図
【図４】（Ａ）図３に示した実施形態における光源装置の他の実施形態を示す図、（Ｂ）図５に示した実施形態における光源装置の他の実施形態を示す図
【図５】本発明の第２のグルコース濃度測定方法を実施するための具体的な実施形態の構成を示す図
【図６】周波数掃引の様子を示すグラフ
【図７】本発明の第３のグルコース濃度測定方法を実施するための具体的な実施形態の構成を示す図
【図８】時間分解可能の光検出器により検出された各光の強度の概念を示すグラフ
【符号の説明】
１０ 光源装置
１１ａ ，…，１１ｘ 光源
２５ａ ，２５ｂ 光音響変調器
２７ ミラー
２８ 光検出器
２９ ヘテロダイン演算手段
３０ 光吸収特性分析手段
４０ グルコース濃度算出手段
５０ 変換手段
１００ 眼球
１１０ 角膜
１２０ 前眼房（眼房水）
１３０ 水晶体
【発明の属する技術分野】
本発明は生体におけるグルコース濃度の測定方法および装置に関し、詳細には眼球部の前眼房水内におけるグルコース濃度を非侵襲的に測定する方法および装置、並びにこの前眼房水内におけるグルコース濃度に基づいて血中グルコース濃度を非侵襲的に測定する方法および装置に関するものである。
【０００２】
【従来の技術】
血液中のグルコース濃度は個人差によりその平均レベルは異なるが、特に糖尿病疾患に対する投薬の要否を決定するための重要な指標値となっている。
【０００３】
ところでこの血中グルコース濃度は、食餌、肉体的活動、その他疾患の併発等によって極短時間の間に大きく変動する特性を有しており、急激な血中グルコース濃度の上昇によって緊急に投薬を要する場合も少なくない。
【０００４】
このため、このような疾患を有する患者についてはなるべく短い間隔で血中グルコース濃度をモニタすることが望まれているが、一方でこの血中グルコース濃度のモニタは通常、指先を穿切して実際に採血し、この血液を分析することにより血液中に含まれるグルコースの濃度（以下、単に血中グルコース濃度という）を測定することにより行なわれており、この穿切に伴う痛みのため患者に１日に何回も測定を強要するのは困難な状況にある。
【０００５】
そこで近年、このような欠点を有する上記侵襲型（侵入型）の測定に代えて、痛み等を伴わない非侵襲型（非侵入型）の測定方法が種々提案されている。
【０００６】
これらは主として、人の眼球部の角膜と水晶体との間にある前眼房を満たす眼房水中のグルコース濃度が、個人差はあるものの血中グルコース濃度と極めて高い相関関係を有していることに着目したものであり、眼房水中のグルコース濃度を非侵襲で測定するものである。
【０００７】
例えば特開昭５１−７５４９８号（米国特許第３，９５８，５６０ 号）によれば、眼房水中に入射した赤外線の旋光度を求めることにより、この旋光度と関連のあるグルコース濃度を得るものである。
【０００８】
また特表平６−５０３２４５ 号によればグルコースの誘導ラマン光を測定するものである。
【０００９】
さらに特開平６−２３７８９８ 号では、水晶体による反射光を光学的性質を測定するデバイスの発明が記載されている。その他、米国特許第５，４３３，１９７ 号においても眼房水中のグルコース濃度を測定する方法が記載されている。
【００１０】
【発明が解決しようとする課題】
しかしながら、上記特開平６−２３７８９８ 号に記載された発明は、角膜と眼房水との境界面での反射光を除去することができず、また角膜での吸収情報も一緒に検出されるため眼房水中のグルコース濃度決定の精度が低下する。
【００１１】
また、微小な吸光度変化を具体的にどのような技術的手段によって測定するかについての開示がないため実用上問題がある。
【００１２】
一方、米国特許第３，９５８，５６０ 号によれば、眼房水中のグルコース以外の多数の化合物が光学的に活性で偏光面の回転に関与する。また角膜は複屈折を示すため偏向面の回転を引き起こす。したがって旋光度により眼房水中のグルコース濃度を測定することは測定精度確保のうえで問題がある。
【００１３】
さらに特表平６−５０３２４５ 号によれば、グルコースの誘導ラマン光の測定のため大きな出力のポンプレーザ光を前眼房部に、視線光軸に対して垂直に導入するため実用的な測定系の構成が困難である。
【００１４】
本発明は上記事情に鑑みなされたものであって、眼房水中のグルコース濃度を非侵襲的に、かつ精度よく測定する装置および方法を提供することを目的とするものである。
【００１５】
また本発明はこのような非侵襲的に得られた眼房水中のグルコース濃度を用いて血液中のグルコース濃度を非侵襲的に測定する装置および方法を提供することをも目的とするものである。
【００１６】
【課題を解決するための手段】
本発明の第１のグルコース濃度測定方法は、
所定の光源から出射された低コヒーレンスな光を、互いに異なる２つの光路に沿って各別に進行する信号光と参照光とに分割し、
該信号光と参照光とで僅かな周波数差が生じるように両光のうち少なくとも一方を変調し、
所定の位置に予め配された眼球に前記信号光を照射し、
該眼球に照射された信号光の、角膜と前眼房との境界面による第１の後方散乱光と、前記参照光とを、該参照光の光路長を調整することにより干渉せしめ、
該干渉により得られた第１の干渉光の強度を測定し、
第１の干渉光の強度に基づいて前記第１の後方散乱光の強度を求め、
前記眼球に照射された信号光の、前眼房と水晶体との境界面による第２の後方散乱光と、前記参照光とを、該参照光の光路長を調整することにより干渉せしめ、
該干渉により得られた第２の干渉光の強度を測定し、
第２の干渉光の強度に基づいて前記第２の後方散乱光の強度を求め、
前記第１および第２の後方散乱光の強度に基づいて、前記前眼房内を満たす眼房水成分の光吸収特性を求め、
前記低コヒーレンスな光とは波長帯域の異なる他の複数の低コヒーレンスな光についても同様に前記光吸収特性を求め、
得られた複数の光吸収特性に基づいて前記眼房水成分中のグルコース濃度を求めることを特徴とするものである。
【００１７】
ここで上記低コヒーレンスな光としては、可干渉距離が数十μｍ程度と短い例えばＳＬＤ（Ｓｕｐｅｒ Ｌｕｍｉｎｅｓｃｅｎｔ Ｄｉｏｄｅ ）やＬＥＤ等が用いられる。なお実用上はより指向性の高いＳＬＤを用いるのが望ましい。
【００１８】
また、上記干渉光の強度を測定するとは、上記後方散乱光（信号光）と参照光との差周波数で強弱を繰り返すビート信号（干渉光）の強度を計測することを意味する。以下の発明においても同様である。
【００１９】
なお、上記本発明の第１のグルコース濃度測定方法において、上記各低コヒーレンスな光は、これら各低コヒーレンスな光の波長帯域よりも広い発光波長帯域の光の一部としてそれぞれ選択された光であってもよいし、互いに異なる複数の光源から各別に出射された光であってもよい。
【００２０】
また、複数の光源により構成した場合には、この波長帯域の異なる複数の光源から順次に低コヒーレンスな光を出射して、これらに対応した各干渉光を１つの光検出器を用いて検出すればよい。
【００２１】
本発明の第２のグルコース濃度測定方法は、
所定の光源から出射された、時間的に鋸歯状（例えば図６に示すよう）に周波数掃引されたコヒーレント光（例えばレーザー光）を互いに異なる２つの光路に沿って各別に進行する信号光と参照光とに分割し、
所定の位置に予め配された眼球に前記信号光を照射し、
該眼球の角膜と前眼房との境界面による前記信号光の第１の後方散乱光と、前記信号光（参照光と分割されてから眼球の角膜と前眼房との境界面に至るまでの信号光を意味する）および前記第１の後方散乱光（眼球の角膜と前眼房との境界面から参照光と干渉するまでのもの）と前記参照光（信号光と分割されてから第１の後方散乱光と干渉するまでのもの）との光路長差に基づく時間差をもって前記光源から出射した、前記第１の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを干渉せしめ、
該干渉により得られた第１の干渉光の強度を測定し、
該第１の干渉光の強度に基づいて前記第１の後方散乱光の強度を求め、
前記眼球の前眼房と水晶体との境界面による前記信号光の第２の後方散乱光と、前記信号光（参照光と分割されてから眼球の前眼房と水晶体との境界面に至るまでの信号光を意味する）および前記第２の後方散乱光（眼球の前眼房と水晶体との境界面から参照光と干渉するまでのもの）と前記参照光（信号光と分割されてから第２の後方散乱光と干渉するまでのもの）との光路長差に基づく時間差をもって前記光源から出射した、前記第２の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを干渉せしめ、
該干渉により得られた第２の干渉光の強度を測定し、
該第２の干渉光の強度に基づいて前記第２の後方散乱光の強度を求め、
前記第１および第２の後方散乱光の強度に基づいて、前記前眼房内を満たす眼房水成分の光吸収特性を求め、
前記コヒーレント光とは波長の異なる他の複数のコヒーレント光についても同様に前記光吸収特性を求め、
得られた複数の光吸収特性に基づいて前記眼房水成分中のグルコース濃度を求めることを特徴とするものである。
【００２２】
ここで上記第１の後方散乱光と干渉せしめられる参照光が「第１の後方散乱光とは周波数差を有する」のは、以下の理由による。
【００２３】
上記信号光が通過する光路長と第１の後方散乱光が通過する光路長との和と、参照光が通過する光路長とには差が設けられており（参照光の光路長の方が短い場合だけでなく、参照光の光路長の方を長く設定することも勿論可能である。）、この光路長差によって、例えば参照光の光路長の方が短い場合には、波面整合（干渉）せしめられる位置に到達するのは参照光の方が第１の後方散乱光よりも早い。
【００２４】
すなわち第１の後方散乱光がその位置に到達したときには、この第１の後方散乱光を生じせしめた信号光と分割された参照光は既にこの位置を通過しており、この信号光が光源から出射された時点よりも遅い時刻に出射されたコヒーレント光の一部である参照光が到達する。
【００２５】
この遅い時刻に出射されたコヒーレント光は、時間的に周波数掃引されているため、第１の後方散乱光とは僅かに周波数差を有するものとなっている。
【００２６】
なお、上記第２の後方散乱光と干渉せしめられる参照光は、第１の後方散乱光と干渉せしめられる参照光とも異なる周波数を有している。これは、第１の後方散乱光が角膜と前眼房との境界による散乱光であるのに対し、第２の後方散乱光は、その境界よりもさらに眼球の深部である前眼房と水晶体との境界による散乱光であるため、両者の光路長差（前眼房の厚さの２倍）分に応じた時間差が生じ、この結果、第２の後方散乱光と干渉せしめられる参照光は第１の後方散乱光と干渉せしめられる参照光よりも遅い時刻に光源から出射されたコヒーレント光によるものとなる。
【００２７】
なお、上記本発明の第２のグルコース濃度測定方法において、上記各コヒーレント光は、単一の光源からそれぞれ選択的に出射された光であってもよいし、または互いに異なる複数の光源から各別に出射された光であってもよい。
【００２８】
本発明の第３のグルコース濃度測定方法は、
所定の光源から出射された超短パルス光を眼球に照射し、
該眼球の角膜と前眼房との境界面による前記超短パルス光の第１の後方散乱光の強度および前眼房と水晶体との境界面による第２の後方散乱光の強度を各別に測定し、
該第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求め、
前記超短パルス光とは波長の異なる他の複数の超短パルス光についても同様に前記光吸収特性を求め、
該超短パルス光の波長ごとに得られた複数の光吸収特性に基づいて前記眼房水成分中のグルコース濃度を求めることを特徴とするものである。
【００２９】
ここで上記超短パルス光とは、少なくとも上記第１の後方散乱光の強度と、前眼房と水晶体との境界面による第２の後方散乱光の強度とを、時間的に分離して各別に測定し得る程度の非常に短時間（例えばフェムト秒〜ピコ秒単位程度）だけ発光するパルス状の光、例えばモードロックＴｉ：サファイアレーザーなど等を意味する。このような超短パルス光を用いることにより、第１の後方散乱光に対して前眼房を往復する距離に対応した時間だけ遅れた第２の後方散乱光を、例えばストリークカメラ等の時間分解可能の光検出器を用いて、第１の後方散乱光から分離して検出することができる。
【００３０】
なおこの方法においても、各超短パルス光は、単一の光源からそれぞれ選択的に出射された光であってもよいし、または、互いに異なる複数の光源から各別に出射された光であってもよい。
【００３１】
本発明の第４のグルコース濃度測定方法は、上述の各グルコース濃度測定方法のうちいずれかにより求められた眼房水成分中のグルコース濃度ごとに、対応する血中グルコース濃度を各別に侵襲的に求めて、その方法により求められた前記眼房水成分中のグルコース濃度と侵襲的に得た血中グルコース濃度との相関を予め求め、
この相関を求めた以後においては、上記各グルコース濃度測定方法のうちいずれかにより求めた眼房水成分中のグルコース濃度と前記相関とに基づいて、血中グルコース濃度を非侵襲的に求めることを特徴とするものである。
【００３２】
本発明の第１のグルコース濃度測定装置は、上記本発明の第１のグルコース濃度測定方法を実施するための装置であって、
互いに発光波長帯域の異なる複数の低コヒーレンスな光を出射する光源装置と、
該光源装置から出射された低コヒーレンスな光を互いに異なる２つの光路に沿って各別に進行する参照光と眼球に入射される信号光とに分割する光路分割手段と、
該信号光と参照光とで僅かな周波数差が生じるように両光のうち少なくとも一方を変調する、該少なくとも一方の光路上に設けられた変調手段と、
前記参照光が進行する光路の長さを調整する光路長調整手段と、
前記眼球の角膜と前眼房との境界面による信号光の第１の後方散乱と前記参照光、および、前記眼球の前眼房と水晶体との境界面による前記信号光の第２の後方散乱光と前記参照光とを、それぞれ波面整合させる波面整合手段と、
前記参照光と前記第１の後方散乱光との波面整合による第１の干渉光および前記参照光と前記第２の後方散乱光との波面整合による第２の干渉光の強度を光電的に検出する光検出器と、
該各干渉光の強度に基づいて前記第１および第２の後方散乱光の強度を求めるヘテロダイン演算手段と、
前記第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段と、
前記複数の低コヒーレンスな光ごとに得られた複数の光吸収特性に基づいて、前記眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段とを備えたことを特徴とするものである。
【００３３】
ここで上記光源装置は、前記各コヒーレンスな光の波長帯域よりも広い発光波長帯域の低コヒーレンスな光を出射する単一の光源と、
該発光波長域の広い低コヒーレンスな光から、前記各低コヒーレンスな光を波長選択する波長選択手段とにより構成されているものであってもよいし、前記各低コヒーレンスな光を各別に出射する複数の光源により構成されているものであってもよい。
【００３４】
本発明の第２のグルコース濃度測定装置は、上記本発明の第２のグルコース濃度測定方法を実施するための装置であって、
図６に示すように、時間的に鋸歯状に周波数掃引された、互いに波長の異なる複数のコヒーレント光を出射する光源装置と、
該光源装置から出射された周波数掃引されたコヒーレント光を互いに異なる２つの光路に沿って各別に進行する参照光と眼球に入射される信号光とに分割する光路分割手段と、
前記眼球の角膜と前眼房との境界面による前記信号光の第１の後方散乱光と、前記信号光および前記第１の後方散乱光と前記参照光との光路長差に基づく時間差をもって前記光源から出射した、前記第１の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを、および、前記眼球の前眼房と水晶体との境界面による前記信号光の第２の後方散乱光と、前記信号光および前記第２の後方散乱光と前記参照光との光路長差に基づく時間差をもって前記光源から出射した、前記第２の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを、それぞれ波面整合させる波面整合手段と、
前記第１の後方散乱光と該第１の後方散乱光に対して僅かな周波数差を有する参照光との波面整合により得られた第１の干渉光、および前記第２の後方散乱光と該第２の後方散乱光に対して僅かな周波数差を有する参照光との波面整合により得られた第２の干渉光の強度を光電的に検出する光検出器と、
該各干渉光の強度に基づいて前記第１および第２の後方散乱光の強度を求めるヘテロダイン演算手段と、
前記第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段と、
前記複数のコヒーレント光ごとに得られた複数の光吸収特性に基づいて、前記眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段とを備えたことを特徴とするものである。
【００３５】
このグルコース濃度測定装置においても、光源装置は、前記複数のコヒーレント光のうちいずれか１つを選択的に出射し得る単一の光源と、該光源に対して、前記複数のコヒーレント光のうち選択的にいずれか１つを出射させる制御を施す制御手段とにより構成されているものであってもよいし、または、前記各コヒーレント光を各別に出射する複数の光源により構成されているものであってもよい。
【００３６】
本発明の第３のグルコース濃度測定装置は、上記本発明の第３のグルコース濃度測定方法を実施するための装置であって、
互いに波長の異なる複数の超短パルス光を出射する光源装置と、
該超短パルス光を眼球に入射せしめ、該眼球の角膜と前眼房との境界面による該超短パルス光の第１の後方散乱光の強度、および前眼房と水晶体との境界面による該超短パルス光の第２の後方散乱光の強度を時系列的に各別に求める光時間領域後方散乱測定手段と、
前記第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段と、
前記複数の超短パルス光ごとに得られた複数の光吸収特性に基づいて、前記眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段とを備えたことを特徴とするものである。
【００３７】
また上記光源装置は、前記複数の超短パルス光のうちいずれか１つを選択的に出射し得る単一の光源と、該光源に対して、前記複数の超短パルス光のうち選択的にいずれか１つを出射させる制御を施す制御手段とにより構成されているものであってもよいし、または、前記各超短パルス光を各別に出射する複数の光源により構成されているものであってもよい。
【００３８】
本発明の第４のグルコース濃度測定装置は、上記各グルコース濃度測定装置のうちいずれかにより求められた眼房水成分中のグルコース濃度と予め例えば侵襲的に得た血中グルコース濃度との相関を示すテーブルを備え、
上記各グルコース濃度測定装置のうちいずれかにより求められた眼房水成分中のグルコース濃度および前記テーブルに基づいて、前記血中グルコース濃度を非侵襲的に求めることを特徴とするものである。
【００３９】
【発明の効果】
本発明の第１および第２のグルコース濃度測定方法並びに本発明の第１および第２のグルコース濃度測定装置によれば、光ヘテロダイン後方散乱測定法を用いることにより、まず角膜と前眼房の境界面による入射光の微弱な第１の後方散乱光ＩＲ２および前眼房と水晶体の境界面による入射光の微弱な第２の後方散乱光ＩＲ３を各別に精度よく検出することができる。
【００４０】
ここで図２に示すように、眼球部への入射光の強度をＩ０ 、空気／角膜間の反射率をＲ１ 、角膜／前眼房間の反射率をＲ２ 、前眼房／水晶体間の反射率をＲ３ 、角膜による入射光の片光路の光学的吸収率をα１ 、前眼房（眼房水）による入射光の片光路の光学的吸収率をα２ とすると、第１の後方散乱光の強度ＩＲ２、第２の後方散乱光の強度ＩＲ３は下記式（１），（２）のように表すことができる。
【００４１】
【００４２】
ここで眼房水にはグルコース以外にもＮａＣｌをはじめとして複数の成分が含まれているため、上記の作用を複数の波長の入射光について繰り返し行ない、多変量解析を含む公知な近赤外分光分析法を適用することにより、眼房水中のグルコースだけの濃度を特定することができる。
【００４３】
このようにして非侵襲的に得られた眼房水中のグルコース濃度と従来の侵襲的に得られた血中グルコース濃度との相関関係を例えば変換テーブル等として各患者ごとに予め準備しておくことにより、その後は血中グルコース濃度を測定することなく眼房水中のグルコース濃度を測定すればよく、患者に苦痛を与えずに繰り返し精度のよい測定を行なうことが可能となる。
【００４４】
本発明の第３のグルコース濃度測定方法および本発明の第３のグルコース濃度測定装置によれば、第１の後方散乱光の強度および第２の後方散乱光の強度を図８に示すように時間的に分離して各別に測定することができ、この図８に示した波形を時間について積分することによりその光量、すなわち上記式（１）に示すＩＲ２および式（２）に示すＩＲ３を求めることができる。
【００４５】
したがって式（３）により、眼房水による入射光の片光路の光学的吸収率α２ を求めることができる。
【００４６】
そして本願の上記発明の場合と同様に、複数の波長の入射光について繰り返し行ない、多変量解析を含む公知な近赤外分光分析法により眼房水中のグルコース濃度を特定することができる。
【００４７】
【発明の実施の形態】
以下、本発明の具体的な実施の形態について図面を用いて説明する。
【００４８】
図１は本発明の第１のグルコース濃度測定方法を実施するための基本的な装置構成を示す図である。
【００４９】
図示のグルコース濃度測定装置は、複数の波長帯域の、可干渉距離が数十μｍ程度の低コヒーレンスな光を出射する光源装置１０と、この光源装置１０から出射された低コヒーレンスな光を眼球１００ に照射せしめ、眼球１００ の角膜１１０ と前眼房１２０ との境界面による入射した光の第１の後方散乱光の強度、および前眼房１２０ と水晶体１３０ との境界面による入射した光の第２の後方散乱光の強度、を光ヘテロダイン後方散乱測定法により求める光ヘテロダイン測定手段２０と、これら第１および第２の後方散乱光の強度に基づいて、前眼房１２０ を満たす眼房水成分の入射光の波長に対する光吸収率を求める光吸収特性分析手段３０と、光源装置１０から出射する光の波長ごとに得られた複数の光吸収率に基づいて、多変量解析を含む公知の近赤外分光分析法により眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段４０を備えた構成である。
【００５０】
ここで図１に示した基本的な構成をより具体化した構成を図３に示す。
【００５１】
図示のグルコース濃度測定装置では、図１に示した光源装置１０として、中心波長がλ１であって波長帯域の広い低コヒーレンスな光を出射するＳＬＤ１１ａ ，中心波長がλ２であって波長帯域の広い低コヒーレンスな光を出射するＳＬＤ１１ｂ ，…，中心波長がλｎであって波長帯域の広い低コヒーレンスな光を出射するＳＬＤ１１ｘ という複数の光源を備え、
光ヘテロダイン測定手段２０として、光源装置１０から出射された低コヒーレンスな光を互いに異なる２つの光路に沿って各別に進行する信号光と参照光とに分割し、また２つの光を波面整合する４つのビームスプリッタ２１ａ ，２１ｂ ，２１ｃ ，２１ｄ と、信号光と参照光とで僅かな周波数差Δｆが生じるように両光を各別に変調する光音響変調器２５ａ ，２５ｂ と、参照光が進行する光路の長さを調整する、光軸方向に移動可能の反射ミラー２７と、眼球１００ に導光された信号光の後方散乱光と参照光とが波面整合して得られた干渉光の強度を測定する光検出器２８と、この干渉光の強度に基づいて光ヘテロダイン後方散乱測定法により後方散乱光の強度を求める光ヘテロダイン演算手段２９とを備えた構成である。
【００５２】
ここで各ＳＬＤ１１ａ ，１１ｂ ，…，１１ｘ は順次に光を出射するものである。なお光音響変調器は必ずしも本実施形態のように、両光（信号光および参照光）をそれぞれ変調するものである必要はなく、少なくとも一方を変調して両光に僅かな周波数差を生じせしめればよい。
【００５３】
次に本実施形態のグルコース濃度測定装置の作用について説明する。
【００５４】
まずＳＬＤ１１ａ から中心波長λ１の低コヒーレンスな光が出射される。このとき他の光源は光を出射しない。出射された光は光路上の最初のビームスプリッタ２１ａ により、信号光と参照光とに分割される。これら信号光および参照光は各光路上にそれぞれ設けられた光音響変調器２５ａ ，２５ｂ により変調され、両光には僅かな周波数差Δｆが生じる。
【００５５】
次に、変調された両光のうち信号光はビームスプリッタ２１ｂ を透過して眼球１００ に入射する。
【００５６】
眼球１００ に入射した光は大きく４つの光に分割される。すなわち、
▲１▼空気／角膜間の境界面で反射する光、
▲２▼角膜／前眼房間の境界面で反射する光、
▲３▼前眼房／水晶体間の境界面で反射する光、
▲４▼水晶体を透過する光、である。
【００５７】
これらのうち後方散乱光として、▲１▼〜▲３▼の光が入射光である信号光とは反対の方向に眼球１００ から出射する。
【００５８】
眼球１００ から出射したこれらの各後方散乱光はビームスプリッタ２１ｂ により反射され次のビームスプリッタ２１ｄ を透過する。
【００５９】
一方、変調された後の参照光は、ビームスプリッタ２１ｃ により反射されてミラー２７に入射し、このミラー２７で反射されて再度ビームスプリッタ２１ｃ に入射し、このビームスプリッタ２１ｃ を透過し、次のビームスプリッタ２１ｄ により反射される。
【００６０】
このとき、ミラー２７をその光軸方向（矢印Ｘ方向）に移動すると、参照光がビームスプリッタ２１ｄ に到達するまでに通過する光路長が変化する。
【００６１】
また、参照光、眼球に入射する信号光、および眼球で反射された後の後方散乱光は、可干渉距離の短い低コヒーレンスな光であるため、前述の３つの後方散乱光のうち、信号光として通過した距離および後方散乱光としてビームスプリッタ２１ｄ に到達するまでの通過距離の和が、参照光のビームスプリッタ２１ｄ に到達するまでの通過距離に略等しい後方散乱光だけが参照光と干渉し、この干渉する両光の差周波数（Δｆ）で強弱を繰り返すビート信号が発生する。
【００６２】
すなわち、
（１）▲１▼の後方散乱光と参照光とが干渉するときは、
ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｂ 間の距離とビームスプリッタ２１ｂ 〜空気／角膜間の境界面間の往復距離とビームスプリッタ２１ｂ 〜ビームスプリッタ２１ｄ 間の距離の和が、ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｃ 間の距離とビームスプリッタ２１ｃ 〜ミラー２７間の境界面間の往復距離とビームスプリッタ２１ｃ 〜ビームスプリッタ２１ｄ 間の距離の和とが略等しく、
（２）▲２▼の後方散乱光が参照光と干渉するときは、
ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｂ 間の距離とビームスプリッタ２１ｂ 〜角膜／前眼房間の境界面間の往復距離とビームスプリッタ２１ｂ 〜ビームスプリッタ２１ｄ 間の距離の和が、ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｃ 間の距離とビームスプリッタ２１ｃ 〜ミラー２７間の境界面間の往復距離とビームスプリッタ２１ｃ 〜ビームスプリッタ２１ｄ 間の距離の和とが略等しく、
（３）▲３▼の後方散乱光が参照光と干渉するときは、
ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｂ 間の距離とビームスプリッタ２１ｂ 〜前眼房／水晶体間の境界面間の往復距離とビームスプリッタ２１ｂ 〜ビームスプリッタ２１ｄ 間の距離の和が、ビームスプリッタ２１ａ 〜ビームスプリッタ２１ｃ 間の距離とビームスプリッタ２１ｃ 〜ミラー２７間の境界面間の往復距離とビームスプリッタ２１ｃ 〜ビームスプリッタ２１ｄ 間の距離の和とが略等しい。
【００６３】
したがって、まずビームスプリッタ２１ｃ に近い位置に配置したミラー２７を、信号光の空気／角膜１１０ 間の境界面による後方散乱光と参照光との干渉によるビート信号が光検出器２８で認められる位置まで、序々にビームスプリッタ２１ｃ から遠ざける方向に移動する。
【００６４】
ここでさらにミラー２７をビームスプリッタ２１ｃ から遠ざける方向に移動するとビート信号は発生しなくなるが、ミラー２７の移動を続けると再度ビート信号が光検出器２８により検出される。この２度目のビート信号は、信号光の角膜１１０ ／前眼房１２０ 間の境界面による後方散乱光（以下、第１の後方散乱光という）と参照光とが干渉することにより発生したものであり、最初のビート信号の検出から角膜１１０ の厚さ分だけミラー２７を移動した時に発生する。
【００６５】
この第１の後方散乱光と参照光とが干渉することにより発生したビート信号を、第１の干渉光としてその強度ＢＲ２を光検出器２８により検出する。
【００６６】
さらにミラー２７をビームスプリッタ２１ｃ から遠ざける方向に移動すると、またビート信号は発生しなくなるが、ミラー２７の移動を続けると３度目のビート信号が光検出器２８により検出される。この３度目のビート信号は、信号光の前眼房１２０ ／水晶体１３０ 間の境界面による後方散乱光（以下、第２の後方散乱光という）と参照光とが干渉することにより発生したものであり、２度目のビート信号の検出から前眼房１２０ の厚さ分だけミラー２７を移動した時に発生する。
【００６７】
この第２の後方散乱光と参照光とが干渉することにより発生したビート信号を第２の干渉光としてその強度ＢＲ３を光検出器２８により検出する。
【００６８】
このようにして得られた上記２つのビート信号の強度ＢＲ２，ＢＲ３は電気信号に変換されてヘテロダイン演算手段２９に入力され、ヘテロダイン演算手段２９はこれらの電気信号に基づいて第１の後方散乱光の強度ＩＲ２および第２の後方散乱光の強度ＩＲ３を算出する。
【００６９】
すなわち各ビート信号の強度ＢＲ２，ＢＲ３はそれぞれ第１の後方散乱光の強度ＩＲ２と参照光の強度との積の平方根、第２の後方散乱光の強度ＩＲ３と参照光の強度との積の平方根であり、参照光の強度については予め測定した値を用いればよく、この結果、各後方散乱光の強度ＩＲ２およびＩＲ３を求めることができる。
【００７０】
ここで図２を参照して、光吸収特性分析手段３０による作用を説明する。
【００７１】
まず眼球１００ への入射光である信号光の強度をＩ０ 、空気／角膜１１０ 間の反射率をＲ１ 、角膜１１０ ／前眼房１２０ 間の反射率をＲ２ 、前眼房１２０ ／水晶体１３０ 間の反射率をＲ３ 、角膜１１０ による入射光の片光路の光学的吸収率をα１ 、前眼房（眼房水）１２０ による入射光の片光路の光学的吸収率をα２ とすると、第１の後方散乱光の強度ＩＲ２、第２の後方散乱光の強度ＩＲ３は下記式（１），（２）のように表すことができる。
【００７２】
【００７３】
ここで、光吸収特性分析手段３０は各後方散乱光の強度ＩＲ２，ＩＲ３の比（ＩＲ３／ＩＲ２）を下記式（３）にしたがって計算する。
【００７４】
【００７５】
ここで眼房水にはグルコース以外にもＮａＣｌをはじめとして複数の成分が含まれているため、上記１つの波長λ１の入射光によって得られた眼房水１２０ の光学的吸収率α２ （この波長λ１の入射光による眼房水１２０ の光学的吸収率α２ を以下、α２ （１） と表記する）のみによっては、眼房水中のグルコース濃度を特定することができない。
【００７６】
そこで次に、２番目のＳＬＤ１１ｂ から中心波長λ２の低コヒーレンスな光が出射される。このとき最初のＳＬＤ１１ａ を含めて他のＳＬＤは光を出射しない。
【００７７】
ＳＬＤ１１ｂ から出射された中心波長λ２の低コヒーレンスな光も、最初の中心波長λ１の光と同様に光路上の最初のビームスプリッタ２１ａ により信号光と参照光とに分割され、信号光は眼球１００ に入射し、参照光はミラー２７に反射され、ビームスプリッタ２１ｄ において得られた干渉光に基づいて、光検出器２８、ヘテロダイン演算手段２９および光吸収特性分析手段３０による前述と同様の作用により、第２の中心波長λ２の入射光による眼房水１２０ の光学的吸収率α２ （２） が求められる。
【００７８】
以下、各ＳＬＤ１１ｃ ，…，ＳＬＤ１１ｘ についても同様の操作を繰り返すことにより、複数の波長λ１，λ２，…，λｎの各入射光による眼房水１２０ の光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） を得ることができる。
【００７９】
このようにして得られた入射光の波長ごとの光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） はグルコース濃度算出手段４０に入力され、グルコース濃度算出手段４０はこれらの複数の光学的吸収率に基づいて、多変量解析を含む公知な近赤外分光分析法により、複数の成分を含む眼房水１２０ からグルコースの濃度を特定する。
【００８０】
このように本実施形態のグルコース濃度測定装置によれば、眼房水中のグルコース濃度を非侵襲的に精度よく測定することができる。
【００８１】
そしてこのようにして得られた眼房水中のグルコース濃度は、予め各患者ごとに上記装置で求めた眼房水中のグルコース濃度と従来の侵襲的に得た血中グルコース濃度との相関関係を例えば変換テーブル等として、眼房水／血中グルコース濃度変換手段５０のメモリ５２に記憶せしめておき、上記グルコース濃度測定装置からある患者の眼房水中グルコース濃度が変換手段５０に入力されたときに、眼房水／血中グルコース濃度変換手段５０の変換演算手段５１がこのメモリ５２に記憶せしめられた変換テーブルを参照して、その患者の眼房水中グルコース濃度に応じた血中グルコース濃度を出力するようにすれば、テーブル作成後においては、眼房水中のグルコース濃度を上記装置により非侵襲的に計測してその値を求め、得られた眼房水中グルコース濃度とテーブルとに基づいて血中グルコース濃度を算出することができ、血中グルコース濃度の非侵襲的測定を実現することができる。
【００８２】
後述する以下の実施形態においても眼房水／血中グルコース濃度変換手段５０を設ければ上記と同様に効果を得ることができる。
【００８３】
なお、本実施形態のグルコース濃度測定装置は、光源装置１０として、出射する光の波長帯域が互いに異なる複数の光源１１ａ ，１１ｂ ，…，１１ｘ を備えたものについて説明したが、本発明のグルコース濃度測定方法，装置は必ずしもこのように複数の光源を備えた構成である必要はなく、例えば図４（Ａ）に示すように、上述の各低コヒーレンスな光の波長帯域の全体を含むような発光波長域のより広い光を出射する１つの光源１１と、この光源１１から出射された発光波長域の広い光から上述の各低コヒーレンスな光（所定の波長帯域のもの）を選択する波長選択手段１２とからなる構成としてもよい。
【００８４】
また本実施形態では、非侵襲的に得られた眼房水中のグルコース濃度を用いて非侵襲的に血液中のグルコース濃度を得るものであり、本発明の第３のグルコース濃度測定方法，装置の実施例形態でもある。
【００８５】
図５は、本発明の第２のグルコース濃度測定方法を実施するための具体的な一形態を示す図である。
【００８６】
図示のグルコース濃度測定装置の基本的な構成は図１に示したものと同様であるが、図１に示した構成の光源装置１０を、波長λ１のレーザー光を出射する光源１１ａ ′，波長λ２のレーザー光を出射する光源１１ｂ ′，…，波長λｎのレーザー光を出射する光源１１ｘ ′と、これらのレーザー光の入力を受けて図６に示すように各光の周波数を一定の帯域で時間的に鋸歯状に掃引する周波数掃引手段１３とを備えた光源装置１０′としたものである。なお、本実施形態においてはコヒーレント光の一形態としてレーザー光を適用したものについて説明するが、本発明の第２のグルコース濃度測定方法（装置）においてはコヒーレント光であれば、特にレーザー光に限るものではない。
【００８７】
さらに図１に示した構成の光ヘテロダイン測定手段２０を、ビームスプリッタ２３、光路分割手段としてのおよび波面整合手段としての反射板２６と、光検出器２８と、干渉光の強度に基づいて後方散乱光の強度を求めるヘテロダイン演算手段２９′とを備えた構成としたものである。
【００８８】
周波数掃引手段１３により周波数掃引されて光源装置１０′から出射したレーザー光はビームスプリッタ２３を透過して眼球部１００ に入射せしめられ、一方、反射板２６は、眼球１００ に照射されたレーザー光の各反射面で反射された後方散乱光と、この後方散乱光よりも時間的に遅れて光源装置１０′から出射されたレーザー光であって後方散乱光に対して僅かな周波数差を有しその反射面２６ａ において反射した光（参照光）とを波面整合して干渉光を生じせしめる作用を果たすものである。
【００８９】
またヘテロダイン演算手段２９′は、反射板２６により波面整合された後方散乱光とこの後方散乱光より遅れて光源装置１０′から出射され、後方散乱光と同時に反射板２６に到達して反射したレーザー光（参照光）との干渉光の、光検出器２８により検出された強度に基づいて、光ヘテロダイン後方散乱測定法により後方散乱光の強度を求める作用をなすものである。
【００９０】
なお、本実施形態のグルコース濃度測定装置においても、図１および図３に示した光吸収特性分析手段３０およびグルコース濃度算出手段４０を備えているが、その作用は図１または図３に示したものと同様であるため図示を省略している。
【００９１】
次に本実施形態のグルコース濃度測定装置の作用について説明する。
【００９２】
まず光源１１ａ ′から波長λ１のレーザー光が出射される。この光は光源装置１０′の周波数掃引手段１３によりその周波数が掃引されて光源装置１０′から出射される。ここで光源装置１０′から時刻ｔ０ において出射された光の周波数をｆ０ とする（図６参照）。
【００９３】
光源装置１０′から出射された周波数ｆ０ の光は、その光路上に配されたビームスプリッタ２３を透過して反射板２６に入射する。この反射板２６はその反射面２６ａ が半透鏡に形成されているため、入射した光の一部を反射し、残りの透過を許容する。この反射した光は参照光として、透過した光は信号光として作用する。
【００９４】
したがって、周波数ｆ０ の光もその大半が反射面２６ａ を透過し、一部が反射面２６ａ で反射される。
【００９５】
ここでまず反射面２６ａ を透過した周波数ｆ０ の光は、眼球１００ に入射し、前述と同様に、
▲１▼空気／角膜間の境界面で反射する光、
▲２▼角膜／前眼房間の境界面で反射する光、
▲３▼前眼房／水晶体間の境界面で反射する光、
▲４▼水晶体を透過する光、に分割される。
【００９６】
そしてこれらのうち後方散乱光として、▲１▼〜▲３▼の光が入射方向とは反対の方向に眼球１００ から出射し、反射板２６に到達する。
【００９７】
ここで、反射面２６ａ を透過した周波数ｆ０ の信号光が眼球１００ に入射し、後方散乱光として再度、反射板２６の反射面２６ａ に到達するまでの光路長について考察すると、
（１）上記▲１▼の後方散乱光の場合は、反射板２６の反射面２６ａ 〜空気／角膜１１０ 間の境界面までの往復距離、
（２）上記▲２▼の後方散乱光の場合は、反射板２６の反射面２６ａ 〜空気／角膜１１０ 間の境界面までの距離と、空気／角膜１１０ 間の境界面〜角膜１１０ ／前眼房１２０ 間の境界面までの距離との和の往復分、
（３）上記▲３▼の後方散乱光の場合は、反射板２６の反射面２６ａ 〜空気／角膜１１０ 間の境界面までの距離と、空気／角膜１１０ 間の境界面〜角膜１１０ ／前眼房１２０ 間の境界面までの距離と、角膜１１０ ／前眼房１２０ 間の境界面〜前眼房１２０ ／水晶体１３０ 間の境界面までの距離との和の往復分、
となる。
【００９８】
したがって、これらの各後方散乱光の光路長（眼球１００ に入射するまでの光路長を含む）は互いに異なるものであり、この結果、各後方散乱光が反射板２６の反射面２６ａ に到達するには、この光路長差に応じた時間差が生じる。
【００９９】
一方、光源装置１０′からは時間的に周波数掃引されたレーザー光が連続的に出射されており、これらの光も連続的に反射板２６に到達する。
【０１００】
この結果、例えば上記▲１▼の後方散乱光（周波数ｆ０ ）は、反射板２６の反射面２６ａ において、時刻ｔ１ （図６参照）に光源装置１０′から出射された周波数ｆ１ の光と出会って波面整合される。したがって、この波面整合による干渉光のビート信号は両光の差周波数（ｆ１ −ｆ０ ）で強弱を繰り返す。
【０１０１】
また、上記▲２▼の後方散乱光（周波数ｆ０ ）は、反射板２６の反射面２６ａ において、時刻ｔ２ に光源装置１０′から出射された周波数ｆ２ の光と出会って波面整合される。したがって、この波面整合による干渉光のビート信号の周波数は両光の差周波数である（ｆ２ −ｆ０ ）となる。
【０１０２】
さらに、上記▲３▼の後方散乱光（周波数ｆ０ ）は、反射板２６の反射面２６ａ において、時刻ｔ３ に光源装置１０′から出射された周波数ｆ３ の光と出会って波面整合される。したがって、この波面整合による干渉光のビート信号の周波数は両光の差周波数である（ｆ３ −ｆ０ ）となる。
【０１０３】
光源装置１０′からは周波数が刻々と変化するレーザー光が連続的に出射されるが、上記各光路長差は不変であるため、上記各後方散乱光との干渉によるビート信号の周波数は常に一定に保たれる。
【０１０４】
そしてこのように反射板２６の反射面２６ａ において各後方散乱光ごとにそれぞれ干渉した光はビームスプリッタ２３により反射されて光検出器２８に入射する。
【０１０５】
光検出器２８は、これら強弱の周波数が互いに異なる３つのビート信号が合成された光の強度を検出し、この強度を電気信号に変換してヘテロダイン演算手段２９′に入力する。
【０１０６】
ヘテロダイン演算手段２９′は、入力された電気信号に基づいて、各後方散乱光ごとのビート信号を分離し、この分離した各ビート信号の強度に基づいて各後方散乱光の強度を求める。
【０１０７】
なお、各後方散乱光ごとのビート信号の分離は、角膜１１０ および前眼房１２０ の光軸方向の厚さが既知であれば、その厚さに応じた光の通過時間の差に基づいて周波数の掃引速度の関係から、各後方散乱光によるビート信号の周波数差を容易に求めることができ、または角膜１１０ および前眼房１２０ の厚さが未知であっても、得られた合成ビート信号について公知の周波数分析を行なうことにより、各後方散乱光によるビート信号の周波数を容易に求めることができる。
【０１０８】
このようにヘテロダイン演算手段２９′により求められた各後方散乱光の強度は光吸収特性分析手段３０に入力される。
【０１０９】
そして、光源装置１０′の２番目以下の光源１１ｂ ′，…，１１ｘ ′についても上述と同様の操作を繰り返すことにより、複数の波長の各入射光による眼房水１２０ の光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） を得ることができ、得られた入射光の波長ごとの光学的吸収率α２ （１） ，α２ （２） ，…，α２ （ｎ） はグルコース濃度算出手段４０に入力され、グルコース濃度算出手段４０はこれらの複数の光学的吸収率に基づいて、多変量解析を含む公知な近赤外分光分析法により、複数の成分を含む眼房水１２０ からグルコースの濃度を特定する。
【０１１０】
このように本実施形態のグルコース濃度測定装置によれば、眼房水中のグルコース濃度を非侵襲的に精度よく測定することができる。
【０１１１】
なお、本実施形態のグルコース濃度測定装置においても光源装置１０′として、各コヒーレント光を各別に出射する複数の光源１１ａ ′，１１ｂ ′，…，１１ｘ ′からなる構成について説明したが、必ずしもこのように各別の光源であるものに限らず、図４（Ｂ）に示すように、これらの各コヒーレント光を選択的に出射し得る単一の光源１１′と、光源１１′に対して、上記複数のコヒーレント光のうち選択的にいずれか１つを出射させる制御を施す制御手段１４とからなるものであってもよい。
【０１１２】
図７は本発明の第３のグルコース濃度測定方法を実施するための装置構成を示す図である。
【０１１３】
図示のグルコース濃度測定装置は、波長λ１の超短パルス光を出射するＴｉ：サファイアレーザー光源（以下、単に光源という）１１ａ ″，波長λ２の超短パルス光を出射する光源１１ｂ ″，…，波長λｎの超短パルス光を出射する１１ｘ ″という複数の超短パルス光を出射する光源装置１０″と、これらの超短パルス光を眼球１１０ に入射せしめ、眼球１００ の角膜１１０ と前眼房１２０ との境界面による超短パルス光の第１の後方散乱光の強度および前眼房１２０ と水晶体１３０ との境界面による超短パルス光の第２の後方散乱光の強度を時系列的に各別に求める、例えばストリークカメラ等の光検出器を含む光時間領域後方散乱測定手段６０と、第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段３０と、複数の超短パルス光ごとに得られた複数の光吸収特性に基づいて、眼房水１２０ 成分中のグルコース濃度を求めるグルコース濃度算出手段４０とを備えた構成である。
【０１１４】
ここで上記超短パルス光は、例えばフェムト秒〜ピコ秒という極めて短い時間だけ発光するパルス状の光を意味する。
【０１１５】
また光時間領域後方散乱測定手段６０は具体的には、光源装置１０″から出射された超短パルス光を、眼球１００ に入射する信号光と入射光による眼球１００ の各後方散乱光を検出するタイミングをとるための参照光とに分割するビームスプリッタ６１と、この参照光を導光するミラー６２と、入射光による眼球１００ の各後方散乱光を図８に示すように時間的に分離して各別に測定することができるストリークカメラ等の光検出器６３と、このように時間的に分離して検出された各後方散乱光の強度を時間について積分することによりその光量を求める眼房水透過率演算手段６４とを備えた構成である。
【０１１６】
なお参照光は、上述した後方散乱光の検出のタイミングをとるだけでなく、眼球１００ へ入射する信号光の光量をモニタするという意義をも有して光検出器６３に入力される。
【０１１７】
次に本実施形態のグルコース濃度測定装置の作用について説明する。
【０１１８】
まず光源１１ａ ″から波長λ１の超短パルス光が出射される。この光は光時間領域後方散乱測定手段６０に入力され、まずビームスプリッタ６１により信号光と参照光とに分割される。
【０１１９】
このうち信号光は眼球１００ に入射する。眼球１００ に入射した信号光は既に各実施形態の説明において既述したように、４つの光に大きく分割される。そしてそれら４つの光のうち、
▲１▼空気／角膜間の境界面で反射する光、
▲２▼角膜／前眼房間の境界面で反射する光、
▲３▼前眼房／水晶体間の境界面で反射する光、がそれぞれ後方散乱光として眼球１００ から出射する。
【０１２０】
このとき、各後方散乱光は信号光の入射方向についての反射面の深さが互いに相違するため、眼球１００ から出射するタイミングが異なり、上記▲１▼の後方散乱光、▲２▼の後方散乱光、▲３▼の後方散乱光、の順に眼球１００ から出射する。そして各後方散乱光は入射した信号光による反射光であるため、いずれも超短パルス光であり、この超短パルス光の発光時間は上記反射面の深さの差を通過する時間差に対して短いため、各後方散乱光は完全に分離されて眼球１００ から順次に出射する。
【０１２１】
このように順次に眼球１００ から出射された各後方散乱光は、ビームスプリッタ６１で反射されて順次に時間分解可能の光検出器６３に入力する。
【０１２２】
この光検出器６３には、各後方散乱光が入射するのに先だって、ビームスプリッタ６１で反射され、各ミラー６２により導光された参照光が入力される。この参照光が光検出器６３による検出タイミングのトリガーとして作用し、この参照光自体の強度が検出されるとともに、各後方散乱光の強度も各別に検出される。
【０１２３】
光検出器６３により時間分解されて時系列的に各別に検出された参照光の強度および各後方散乱光の強度の様子を図８に示す。
【０１２４】
光検出器６３はこのように検出した各光の強度を電気信号に変換して眼房水透過率演算手段６４に入力する。
【０１２５】
眼房水透過率演算手段６４は入力された各光の強度に対応した電気信号を時間について積分することにより各光の光量を求める。求められた各光の光量は光吸収特性分析手段３０に入力され、光吸収特性分析手段３０はこの各光の光量に応じて前記各実施形態の場合と同様の作用により、眼房水１２０ の光学的吸収率を求める。
【０１２６】
以上の作用と同様の操作を、上記光源装置１０″から出射する超短パルス光の波長を変えて、すなわち互いに異なる単一波長のパルス光を出射する複数の光源１１ａ ″，１１ｂ ″，…，１１ｘ ″を順次に代えることによって行ない、各入射光の波長に応じた眼房水１２０ の複数の光学的吸収率を求める。
【０１２７】
以上の操作により得られた入射光の波長ごとの光学的吸収率は、グルコース濃度算出手段４０に入力され、グルコース濃度算出手段４０はこれらの複数の光学的吸収率に基づいて、多変量解析を含む公知な近赤外分光分析法により、複数の成分を含む眼房水１２０ からグルコースの濃度を特定する。
【０１２８】
このように本実施形態のグルコース濃度測定装置によれば、眼房水中のグルコース濃度を非侵襲的に精度よく測定することができる。
【０１２９】
なお、本実施形態のグルコース濃度測定装置においても、光源装置１０″として、出射する超短パルス光の波長が互いに異なる複数の光源１１ａ ″，１１ｂ ″，…，１１ｘ ″からなる構成について説明したが、必ずしも複数の光源を備えた構成である必要はなく、前述の第２の実施形態の光源装置１０′の場合と同様に、これらの各超短パルス光を選択的に出射し得る単一の光源と、この光源に対して、上記複数の超短パルス光のうち選択的にいずれか１つを出射させる制御を施す制御手段とからなるものであってもよい。
【図面の簡単な説明】
【図１】本発明の第１のグルコース濃度測定方法を実施するための基本的な装置構成を示す図
【図２】眼球への入射光と後方散乱光との関係を示す図
【図３】本発明の第１のグルコース濃度測定方法を実施するためのより具体的な装置構成を示す図
【図４】（Ａ）図３に示した実施形態における光源装置の他の実施形態を示す図、（Ｂ）図５に示した実施形態における光源装置の他の実施形態を示す図
【図５】本発明の第２のグルコース濃度測定方法を実施するための具体的な実施形態の構成を示す図
【図６】周波数掃引の様子を示すグラフ
【図７】本発明の第３のグルコース濃度測定方法を実施するための具体的な実施形態の構成を示す図
【図８】時間分解可能の光検出器により検出された各光の強度の概念を示すグラフ
【符号の説明】
１０ 光源装置
１１ａ ，…，１１ｘ 光源
２５ａ ，２５ｂ 光音響変調器
２７ ミラー
２８ 光検出器
２９ ヘテロダイン演算手段
３０ 光吸収特性分析手段
４０ グルコース濃度算出手段
５０ 変換手段
１００ 眼球
１１０ 角膜
１２０ 前眼房（眼房水）
１３０ 水晶体
Claims (4)
- 時間的に鋸歯状に周波数掃引された、互いに波長の異なる複数のコヒーレント光を順次出射する光源装置と、
該光源装置から出射された周波数掃引されたコヒーレント光を互いに異なる２つの光路に沿って各別に進行する参照光と眼球に入射される信号光とに分割する光路分割手段と、
前記眼球の角膜と前眼房との境界面による前記信号光の第１の後方散乱光と、前記信号光および前記第１の後方散乱光と前記参照光との光路長差に基づく時間差をもって前記光源から出射した、前記第１の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを、および、前記眼球の前眼房と水晶体との境界面による前記信号光の第２の後方散乱光と、前記信号光および前記第２の後方散乱光と前記参照光との光路長差に基づく時間差をもって前記光源から出射した、前記第２の後方散乱光とは周波数差を有するコヒーレント光による前記参照光とを、それぞれ波面整合させる波面整合手段と、
前記第１の後方散乱光と該第１の後方散乱光に対して僅かな周波数差を有する参照光との波面整合により得られた第１の干渉光、および前記第２の後方散乱光と該第２の後方散乱光に対して僅かな周波数差を有する参照光との波面整合により得られた第２の干渉光の強度を光電的に検出する光検出器と、
該各干渉光の強度に基づいて前記第１および第２の後方散乱光の強度を求めるヘテロダイン演算手段と、
前記第１および第２の後方散乱光の強度に基づいて前記前眼房を満たす眼房水成分の光吸収特性を求める光吸収特性分析手段と、
前記複数のコヒーレント光ごとに得られた複数の光吸収特性に基づいて、前記眼房水成分中のグルコース濃度を求めるグルコース濃度算出手段とを備えたことを特徴とするグルコース濃度測定装置。 - 前記光源装置が、前記複数のコヒーレント光のうちいずれか１つを選択的に出射し得る単一の光源と、該光源に対して、前記複数のコヒーレント光のうち選択的にいずれか１つを出射させる制御を施す制御手段とにより構成されていることを特徴とする請求項１記載のグルコース濃度測定装置。
- 前記光源装置が、前記各コヒーレント光を各別に出射する複数の光源により構成されていることを特徴とする請求項１記載のグルコース濃度測定装置。
- 請求項１から３のうちいずれか１項に記載のグルコース濃度測定装置により求められた眼房水成分中のグルコース濃度と予め求めた血中グルコース濃度との相関を示すテーブルを備え、
請求項１から３のうちいずれか１項に記載のグルコース濃度測定装置により求められた眼房水成分中のグルコース濃度および前記テーブルに基づいて、前記血中グルコース濃度を非侵襲的に求めることを特徴とするグルコース濃度測定装置。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12179096A JP3604231B2 (ja) | 1996-05-16 | 1996-05-16 | グルコース濃度測定方法および装置 |
| EP97107946A EP0807812B1 (en) | 1996-05-16 | 1997-05-15 | Glucose concentration measuring method and apparatus |
| DE69712255T DE69712255T2 (de) | 1996-05-16 | 1997-05-15 | Verfahren und Vorrichtung zur Messung der Glukosekonzentration |
| US08/857,205 US5835215A (en) | 1996-05-16 | 1997-05-15 | Glucose concentration measuring method and apparatus with short coherence source and heterodyne interferometer |
| US09/138,563 US5969815A (en) | 1996-05-16 | 1998-08-24 | Glucose concentration measuring method and apparatus with a coherent source and heterodyne interferometer |
| US09/138,124 US5961449A (en) | 1996-05-16 | 1998-08-24 | Glucose concentration measuring method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12179096A JP3604231B2 (ja) | 1996-05-16 | 1996-05-16 | グルコース濃度測定方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09299333A JPH09299333A (ja) | 1997-11-25 |
| JP3604231B2 true JP3604231B2 (ja) | 2004-12-22 |
Family
ID=14819971
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP12179096A Expired - Lifetime JP3604231B2 (ja) | 1996-05-16 | 1996-05-16 | グルコース濃度測定方法および装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (3) | US5835215A (ja) |
| EP (1) | EP0807812B1 (ja) |
| JP (1) | JP3604231B2 (ja) |
| DE (1) | DE69712255T2 (ja) |
Families Citing this family (88)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6240306B1 (en) | 1995-08-09 | 2001-05-29 | Rio Grande Medical Technologies, Inc. | Method and apparatus for non-invasive blood analyte measurement with fluid compartment equilibration |
| US6405069B1 (en) | 1996-01-31 | 2002-06-11 | Board Of Regents, The University Of Texas System | Time-resolved optoacoustic method and system for noninvasive monitoring of glucose |
| JP3604231B2 (ja) * | 1996-05-16 | 2004-12-22 | 富士写真フイルム株式会社 | グルコース濃度測定方法および装置 |
| US20010048738A1 (en) * | 1997-04-03 | 2001-12-06 | Sbc Technology Resourses, Inc. | Profile management system including user interface for accessing and maintaining profile data of user subscribed telephony services |
| US7890158B2 (en) | 2001-06-05 | 2011-02-15 | Lumidigm, Inc. | Apparatus and method of biometric determination using specialized optical spectroscopy systems |
| JP3543923B2 (ja) * | 1997-12-25 | 2004-07-21 | 富士写真フイルム株式会社 | グルコース濃度測定装置 |
| US6137585A (en) * | 1998-05-15 | 2000-10-24 | Laser Diagnostic Technologies, Inc. | Method and apparatus for recording three-dimensional distribution of light backscattering potential in transparent and semi-transparent structures |
| US6151127A (en) * | 1998-05-28 | 2000-11-21 | The General Hospital Corporation | Confocal microscopy |
| DE19830541B9 (de) * | 1998-07-08 | 2005-04-21 | Schrader, Bernhard, Prof. Dr.-Ing. | Verfahren und Vorrichtung zur automatischen Messung der Konzentration von Komponenten des Augen-Vorderkammerwassers sowie Verwendung derselben |
| JP2000037355A (ja) | 1998-07-24 | 2000-02-08 | Fuji Photo Film Co Ltd | グルコース濃度測定方法および装置 |
| WO2000016692A1 (en) * | 1998-09-18 | 2000-03-30 | Proniewicz Walter K | Noninvasive measurement of blood sugar based on optoelectronic observations of the eye |
| US6853854B1 (en) | 1998-09-18 | 2005-02-08 | Q Step Technologies, Llc | Noninvasive measurement system |
| AU6783900A (en) * | 1999-08-17 | 2001-03-13 | Board Of Regents, The University Of Texas System | Methods for noninvasive analyte sensing |
| US6725073B1 (en) | 1999-08-17 | 2004-04-20 | Board Of Regents, The University Of Texas System | Methods for noninvasive analyte sensing |
| AU3213500A (en) * | 1999-09-17 | 2001-04-24 | Walter K. Proniewicz | Noninvasive measurement of blood sugar level |
| US6816605B2 (en) | 1999-10-08 | 2004-11-09 | Lumidigm, Inc. | Methods and systems for biometric identification of individuals using linear optical spectroscopy |
| DE19955268C2 (de) * | 1999-11-17 | 2001-09-06 | Isis Optronics Gmbh | Verfahren und Vorrichtung zur interferometrischen Untersuchung streuender Objekte |
| US6574490B2 (en) | 2001-04-11 | 2003-06-03 | Rio Grande Medical Technologies, Inc. | System for non-invasive measurement of glucose in humans |
| US6704588B2 (en) | 2001-06-16 | 2004-03-09 | Rafat R. Ansari | Method and apparatus for the non-invasive measurement of blood glucose levels in humans |
| US6650915B2 (en) | 2001-09-13 | 2003-11-18 | Fovioptics, Inc. | Non-invasive measurement of blood analytes using photodynamics |
| US20040027582A1 (en) * | 2002-08-09 | 2004-02-12 | Lev Dulman | Method and apparatus for determining sample composition with an interferometer |
| US20040027583A1 (en) * | 2002-08-09 | 2004-02-12 | Lev Dulman | Pre-established reference scale for interferometric topological metrology |
| US6999181B2 (en) * | 2002-08-09 | 2006-02-14 | Angstrovision, Inc. | Advanced signal processing technique for translating fringe line disturbances into sample height at a particular position above an interferometer's sample stage |
| US7136168B2 (en) * | 2002-08-09 | 2006-11-14 | Angstrovision, Inc. | Interferometric topological metrology with pre-established reference scale |
| US6895264B2 (en) | 2002-08-26 | 2005-05-17 | Fovioptics Inc. | Non-invasive psychophysical measurement of glucose using photodynamics |
| DE10260256B9 (de) * | 2002-12-20 | 2007-03-01 | Carl Zeiss | Interferometersystem und Meß-/Bearbeitungswerkzeug |
| WO2004060154A1 (en) * | 2003-01-07 | 2004-07-22 | Intelligent Photonics Control Corp. | Non-invasive blood monitor |
| DE10302849A1 (de) * | 2003-01-23 | 2004-07-29 | Carl Zeiss Meditec Ag | Verfahren und Anordnung zur Messung der Dispersion in transparenten Medien |
| KR20060002923A (ko) | 2003-04-04 | 2006-01-09 | 루미다임 인크. | 다중 스펙트럼 생체인식 센서 |
| US7460696B2 (en) | 2004-06-01 | 2008-12-02 | Lumidigm, Inc. | Multispectral imaging biometrics |
| US7668350B2 (en) | 2003-04-04 | 2010-02-23 | Lumidigm, Inc. | Comparative texture analysis of tissue for biometric spoof detection |
| US7751594B2 (en) | 2003-04-04 | 2010-07-06 | Lumidigm, Inc. | White-light spectral biometric sensors |
| US6968222B2 (en) | 2003-05-02 | 2005-11-22 | Oculir, Inc. | Methods and device for non-invasive analyte measurement |
| US6958039B2 (en) | 2003-05-02 | 2005-10-25 | Oculir, Inc. | Method and instruments for non-invasive analyte measurement |
| US6975892B2 (en) | 2003-10-21 | 2005-12-13 | Oculir, Inc. | Methods for non-invasive analyte measurement from the conjunctiva |
| US7307734B2 (en) * | 2003-08-14 | 2007-12-11 | University Of Central Florida | Interferometric sensor for characterizing materials |
| GB2407378B (en) | 2003-10-24 | 2006-09-06 | Lein Applied Diagnostics Ltd | Ocular property measuring apparatus and method therefor |
| GB2409033C (en) * | 2003-12-12 | 2006-05-24 | Lein Applied Diagnostics Ltd | Extended focal region measuring apparatus and method |
| JP4505852B2 (ja) * | 2004-04-13 | 2010-07-21 | 学校法人早稲田大学 | 眼底分光像撮影装置 |
| US20080009688A1 (en) * | 2004-04-14 | 2008-01-10 | Oculir, Inc. | Methods for non-invasive analyte measurement |
| US7184148B2 (en) | 2004-05-14 | 2007-02-27 | Medeikon Corporation | Low coherence interferometry utilizing phase |
| US7327463B2 (en) | 2004-05-14 | 2008-02-05 | Medrikon Corporation | Low coherence interferometry utilizing magnitude |
| US7190464B2 (en) * | 2004-05-14 | 2007-03-13 | Medeikon Corporation | Low coherence interferometry for detecting and characterizing plaques |
| US7474408B2 (en) | 2004-05-14 | 2009-01-06 | Medeikon Corporation | Low coherence interferometry utilizing phase |
| US20050254059A1 (en) * | 2004-05-14 | 2005-11-17 | Alphonse Gerard A | Low coherence interferometric system for optical metrology |
| US7242480B2 (en) * | 2004-05-14 | 2007-07-10 | Medeikon Corporation | Low coherence interferometry for detecting and characterizing plaques |
| US8229185B2 (en) | 2004-06-01 | 2012-07-24 | Lumidigm, Inc. | Hygienic biometric sensors |
| US8787630B2 (en) | 2004-08-11 | 2014-07-22 | Lumidigm, Inc. | Multispectral barcode imaging |
| US7751862B2 (en) * | 2004-08-19 | 2010-07-06 | Fp Technology | Frequency resolved imaging system |
| JP4563130B2 (ja) * | 2004-10-04 | 2010-10-13 | 株式会社トプコン | 光画像計測装置 |
| US7248907B2 (en) * | 2004-10-23 | 2007-07-24 | Hogan Josh N | Correlation of concurrent non-invasively acquired signals |
| US7167736B2 (en) * | 2004-11-04 | 2007-01-23 | Q Step Technologies, Inc. | Non-invasive measurement system and method for measuring the concentration of an optically-active substance |
| GB2422660C (en) * | 2005-01-27 | 2018-04-25 | H Icheck Ltd | Improved device for monitoring body functions |
| US7801338B2 (en) | 2005-04-27 | 2010-09-21 | Lumidigm, Inc. | Multispectral biometric sensors |
| US7627357B2 (en) * | 2005-06-30 | 2009-12-01 | General Electric Company | System and method for non-invasive glucose monitoring |
| DE102005034332A1 (de) * | 2005-07-22 | 2007-01-25 | Carl Zeiss Meditec Ag | Einrichtung und Verfahren zur Beobachtung, Dokumentation und/oder Diagnose des Augenhintergrundes |
| US7488930B2 (en) | 2006-06-02 | 2009-02-10 | Medeikon Corporation | Multi-channel low coherence interferometer |
| US8355545B2 (en) | 2007-04-10 | 2013-01-15 | Lumidigm, Inc. | Biometric detection using spatial, temporal, and/or spectral techniques |
| WO2008100329A2 (en) | 2006-07-19 | 2008-08-21 | Lumidigm, Inc. | Multibiometric multispectral imager |
| US8175346B2 (en) | 2006-07-19 | 2012-05-08 | Lumidigm, Inc. | Whole-hand multispectral biometric imaging |
| US7995808B2 (en) | 2006-07-19 | 2011-08-09 | Lumidigm, Inc. | Contactless multispectral biometric capture |
| US7804984B2 (en) | 2006-07-31 | 2010-09-28 | Lumidigm, Inc. | Spatial-spectral fingerprint spoof detection |
| US7801339B2 (en) | 2006-07-31 | 2010-09-21 | Lumidigm, Inc. | Biometrics with spatiospectral spoof detection |
| DE102006048849B4 (de) * | 2006-10-16 | 2018-05-09 | Carl Zeiss Ag | Verfahren und Vorrichtung zum Messen der Konzentration zumindest eines Stoffes aus einer Gruppe von n sich in einem Probenmaterial befindenden und die Dispersion des Probenmaterials beeinflussenden Stoffen |
| JP2008104628A (ja) * | 2006-10-25 | 2008-05-08 | Tokyo Institute Of Technology | 眼球の結膜強膜撮像装置 |
| DE102006054556A1 (de) * | 2006-11-20 | 2008-05-21 | Zimmer Medizinsysteme Gmbh | Vorrichtung und Verfahren zum nicht-invasiven, optischen Erfassen von chemischen und physikalischen Blutwerten und Körperinhaltsstoffen |
| CN100512752C (zh) * | 2007-01-12 | 2009-07-15 | 天津市先石光学技术有限公司 | 光程可选择的无创人体成分检测方法及检测装置 |
| WO2008134135A2 (en) | 2007-03-21 | 2008-11-06 | Lumidigm, Inc. | Biometrics based on locally consistent features |
| GB2451443B (en) | 2007-07-30 | 2012-12-26 | Lein Applied Diagnostics Ltd | Optical measurement apparatus and method therefor |
| GB2451441B (en) | 2007-07-30 | 2012-07-11 | Lein Applied Diagnostics Ltd | Optical alignment apparatus and method thereof |
| GB2457302B (en) | 2008-02-11 | 2013-04-10 | Lein Applied Diagnostics Ltd | Measurement apparatus and method therefor |
| DE102008013821B4 (de) * | 2008-03-10 | 2010-11-18 | Westphal, Peter, Dr. | Verfahren und Vorrichtung zur Messung gelöster Stoffe im menschlichen oder tierischen Augen-Kammerwasser |
| GB2467340B (en) | 2009-01-30 | 2013-11-13 | Lein Applied Diagnostics Ltd | Signal sample trigger apparatus, data acquisition system and method of sampling an analogue signal |
| WO2011028620A1 (en) | 2009-08-26 | 2011-03-10 | Lumidigm, Inc. | Multiplexed biometric imaging and dual-imager biometric sensor |
| JP4973751B2 (ja) * | 2009-11-26 | 2012-07-11 | 横河電機株式会社 | 生体成分測定装置 |
| US8570149B2 (en) | 2010-03-16 | 2013-10-29 | Lumidigm, Inc. | Biometric imaging using an optical adaptive interface |
| JP5506504B2 (ja) * | 2010-03-31 | 2014-05-28 | キヤノン株式会社 | 撮像装置及び撮像方法 |
| US8639306B2 (en) | 2011-02-07 | 2014-01-28 | Brien Holden Vision Diagnostics, Inc. | Noninvasive eye-property monitoring, including aqueous-humor glucose monitoring as an indication of blood glucose level |
| EP2680741A4 (en) | 2011-03-02 | 2014-08-27 | Brien Holden Vision Diagnostics Inc | RETINOGRAPH SYSTEM AND METHODOLOGY |
| US9662004B2 (en) | 2011-04-29 | 2017-05-30 | Taiwan Biophotonic Corporation | Apparatus for non-invasive glucose monitoring |
| EP2517624A3 (en) | 2011-04-29 | 2013-03-20 | Industrial Technology Research Institute | Apparatus and method for non-invasive blood glucose monitoring and method for analysing biological molecule |
| US9724022B2 (en) | 2011-04-29 | 2017-08-08 | Taiwan Biophotonic Corporation | Apparatus for non-invasive glucose monitoring |
| TWI495864B (zh) * | 2012-04-27 | 2015-08-11 | Taiwan Biophotonic Corp | 非侵入式血糖監測裝置與方法以及生化分子的分析方法 |
| US9335265B2 (en) * | 2014-02-07 | 2016-05-10 | Agilent Technologies, Inc. | Spectrographic system utilizing a chirped, pulsed optical source |
| US9459201B2 (en) | 2014-09-29 | 2016-10-04 | Zyomed Corp. | Systems and methods for noninvasive blood glucose and other analyte detection and measurement using collision computing |
| KR102335739B1 (ko) | 2014-12-19 | 2021-12-06 | 삼성전자주식회사 | 비 침습적 혈당 측정 방법 및 이를 위한 장치 |
| 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 |
| ES2774983B2 (es) | 2019-01-22 | 2021-06-10 | Univ Sevilla | Dispositivo portable y metodo para la estimacion no invasiva del nivel de glucosa en sangre |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3958560A (en) * | 1974-11-25 | 1976-05-25 | Wayne Front March | Non-invasive automatic glucose sensor system |
| US5025785A (en) * | 1982-09-10 | 1991-06-25 | Weiss Jeffrey N | Diabetes detection method |
| US4834111A (en) * | 1987-01-12 | 1989-05-30 | The Trustees Of Columbia University In The City Of New York | Heterodyne interferometer |
| CA2013768C (en) * | 1990-04-03 | 1995-12-05 | Joco Djurdjevic | Gas injection system |
| US5209231A (en) * | 1990-11-02 | 1993-05-11 | University Of Connecticut | Optical glucose sensor apparatus and method |
| US5243983A (en) * | 1990-12-14 | 1993-09-14 | Georgia Tech Research Corporation | Non-invasive blood glucose measurement system and method using stimulated raman spectroscopy |
| WO1992019930A1 (en) * | 1991-04-29 | 1992-11-12 | Massachusetts Institute Of Technology | Method and apparatus for optical imaging and measurement |
| US5433197A (en) * | 1992-09-04 | 1995-07-18 | Stark; Edward W. | Non-invasive glucose measurement method and apparatus |
| DE4243142A1 (de) * | 1992-12-19 | 1994-06-23 | Boehringer Mannheim Gmbh | Vorrichtung zur in-vivo-Bestimmung einer optischen Eigenschaft des Kammerwassers des Auges |
| US5533509A (en) * | 1993-08-12 | 1996-07-09 | Kurashiki Boseki Kabushiki Kaisha | Method and apparatus for non-invasive measurement of blood sugar level |
| JPH10500338A (ja) * | 1994-05-19 | 1998-01-13 | ベーリンガー・マンハイム・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング | バイオロジカル・サンプル中の分析物を決定するための方法および装置 |
| US5553617A (en) * | 1995-01-20 | 1996-09-10 | Hughes Aircraft Company | Noninvasive method and apparatus for determining body chemistry |
| JP3604231B2 (ja) * | 1996-05-16 | 2004-12-22 | 富士写真フイルム株式会社 | グルコース濃度測定方法および装置 |
| TW342447B (en) * | 1996-11-11 | 1998-10-11 | Cherng Jou | Noninvasive polarized common path optical heterodyne glucose monitoring system |
-
1996
- 1996-05-16 JP JP12179096A patent/JP3604231B2/ja not_active Expired - Lifetime
-
1997
- 1997-05-15 DE DE69712255T patent/DE69712255T2/de not_active Expired - Lifetime
- 1997-05-15 EP EP97107946A patent/EP0807812B1/en not_active Expired - Lifetime
- 1997-05-15 US US08/857,205 patent/US5835215A/en not_active Expired - Lifetime
-
1998
- 1998-08-24 US US09/138,563 patent/US5969815A/en not_active Expired - Lifetime
- 1998-08-24 US US09/138,124 patent/US5961449A/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| EP0807812B1 (en) | 2002-05-02 |
| DE69712255T2 (de) | 2002-08-29 |
| US5969815A (en) | 1999-10-19 |
| US5961449A (en) | 1999-10-05 |
| EP0807812A1 (en) | 1997-11-19 |
| JPH09299333A (ja) | 1997-11-25 |
| US5835215A (en) | 1998-11-10 |
| DE69712255D1 (de) | 2002-06-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3604231B2 (ja) | グルコース濃度測定方法および装置 | |
| JP3543923B2 (ja) | グルコース濃度測定装置 | |
| JP2000037355A (ja) | グルコース濃度測定方法および装置 | |
| JP5215664B2 (ja) | 光学コヒーレンス撮像のシステムおよび方法 | |
| CA2700133C (en) | Three diode optical bridge system | |
| US8632181B2 (en) | SS OCT interferometry for measuring a sample | |
| US20180070830A1 (en) | Systems and methods for time-resolved diffuse correlation spectroscopy | |
| US20100016732A1 (en) | Apparatus and method for neural-signal capture to drive neuroprostheses or control bodily function | |
| US8728062B2 (en) | Device and method for controlling a laser therapy of the eye | |
| US20140268163A1 (en) | Methods and Apparatus Related to Multi Wavelength Photothermal Optical Coherence Tomography | |
| EP2259713B1 (en) | Apparatus and method using light retro-reflected from a retina to non- invasively measure the blood concentration of a substance | |
| RU2007147839A (ru) | Датчик глюкозы | |
| WO2006009572A2 (en) | A non-invasive analysis system | |
| JP3731784B2 (ja) | グルコース濃度測定方法および装置 | |
| US20130163003A1 (en) | Apparatus and method for optical swept-source coherence tomography | |
| JPH10111238A (ja) | 散乱体の吸収情報の計測方法及び装置 | |
| WO2010009452A1 (en) | Method and apparatus for neural-signal capture to drive neuroprostheses or control bodily function | |
| US20240023846A1 (en) | Transabdominal fetal oximetry based on frequency-modulated continuous-wave near-infrared spectroscopy | |
| JP4477568B2 (ja) | 成分濃度測定装置及び成分濃度測定装置制御方法 | |
| US20220007975A1 (en) | Portable Device and Method for Non-Invasive Blood Glucose Level Estimation | |
| JP2003202287A (ja) | 散乱吸収体測定方法及び装置 | |
| US10779755B2 (en) | Non-invasive blood sugar measurement method and device using optical reflectometry | |
| JP2017140159A (ja) | 赤外分光法を用いた非侵襲血糖値測定方法 | |
| JPH11244243A (ja) | グルコース濃度測定方法および装置 | |
| JP2018175481A (ja) | 眼球の光計測装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040301 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040309 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040422 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20040518 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040720 |
|
| A911 | Transfer of reconsideration by examiner before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20040831 |
|
| 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: 20040928 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040928 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071008 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: 20071008 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081008 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091008 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101008 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111008 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121008 Year of fee payment: 8 |