US 7,184,148 B2Grant

Low coherence interferometry utilizing phase

Issue Date:2007-02-27

•25 Claims

•19 Drawing Sheets

Abstract

A method for determining a characteristic of an analyte in a biological sample, the method comprising: directing broadband light by means of a sensing light path at the biological sample, at a target depth defined by the sensing light path and a reference light path; and receiving the broadband light reflected from the biological sample by means of the sensing light path. The method also includes: directing the broadband light by means of the reference light path at a reflecting device; receiving the broadband light reflected from the reflecting device by means of the reference light path; and interfering the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device. The method further includes: detecting the broadband light resulting from interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device; and modulating an effective light path length of at least one of the reference light path and the sensing light path to enhance interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device. The method continues with: determining a magnitude of change of the effective light path length introduced by the modulating when interference is enhanced; determining a variation in an index of refraction from a ratio of the magnitude of change of the effective light path length and the target depth; and determining the characteristic of the analyte in the biological sample from the variation in the index of refraction.

Metadata

Assignee

Medeikon Corporation

Inventor

Gerard A. Alphonse

Application Information

Application Number:US 10/845,849

Filing Date:2004-05-14

Priority Date:2004-05-14

Art Unit:2877

Classifications

IPC:

G01B9/02

Patent Drawings (19 sheets)

Description

Background

[0001] The invention concerns a method for low coherence interferometry of a biological sample using phase. The term “biological sample” denotes a body fluid or tissue of an organism. Biological samples are generally optically heterogeneous, that is, they contain a plurality of scattering centers scattering irradiated light. In the case of biological tissue, especially skin tissue, the cell walls and other intra-tissue components form the scattering centers.

[0002] Generally, for the qualitative and quantitative analysis in such biological samples, reagents or systems of reagents are used that chemically react with the particular component(s) to be determined. The reaction results in a physically detectable change in the solution of reaction, for instance a change in its color, which can be measured as a measurement quantity. By calibrating with standard samples of known concentration, a correlation is determined between the values of the measurement quantity measured at different concentrations and the particular concentration. These procedures allow accurate and sensitive analyses, but on the other hand they require removing a liquid sample, especially a blood sample, from the body for the analysis (“invasive analysis”).

[0003] The American Diabetes Association (ADA) estimates that diabetes afflicts nearly 17 million people in the United States. Diabetes can lead to severe complications over time, including heart failure, kidney failure, blindness, and loss of limb due to poor peripheral circulation. According to ADA, complications arising from diabetes cost the U.S. health care system in excess of $132 Billion.

[0004] Diabetes complications are largely due to years of poor blood glucose control. The Diabetes Care and Complications Trial (DCCT) carried out by the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) demonstrated that more frequent monitoring of blood glucose and insulin levels can prevent many of the long-term complications of diabetes.

[0005] Monitoring of blood glucose concentration is key to managing the therapy of diabetes patients. Monitoring results are used to adjust nutrition, medication, and exercise in order to achieve the best possible glucose control, reducing the complications and mortality associated with diabetes. At present, the most widely used method for monitoring of blood glucose by diabetes patients involves chemical analysis of blood samples taken by puncturing the finger or forearm. This method is painful, requires relatively complex operations, is inconvenient due to disruption of daily life, and may become difficult to perform in the long term due to calluses on the fingers and poor circulation. As a result, the average diabetic patient tests his/her blood glucose levels less than twice a day versus the recommended four or more times per day. Non-invasive blood glucose monitoring techniques with accuracies equal to or better than the current chemical glucose methods are therefore needed.

[0006] Accordingly, a number of procedures and apparatus have been suggested to determine glucose in blood, tissue and other biological samples in vivo and in a non-invasive manner. Existing non-invasive procedures for glucose determination include nuclear magnetic resonance (NMR), electron spin resonance (ESR) and infrared spectroscopy. However, none of these procedures have achieved practical significance. Large and costly equipment is required, which are wholly unsuitable for routine analysis or even for patient self-checking (home monitoring).

[0007] One of the most promising approaches for non-invasive glucose monitoring is based on optical techniques. Optical glucose monitoring techniques are particularly attractive in that they are relatively fast, use non-ionizing radiation, and generally do not require consumable reagents. Several optical glucose monitoring techniques have been proposed so far, with varying degrees of success. Several of these techniques are discussed herein as background, however, once again, none of these techniques has attained significant commercial success relative to invasive techniques.

[0008] One approach is Near-Infrared (NIR)/Mid-Infrared (MIR) spectroscopy. In infrared spectroscopy, radiation from external light sources is transmitted through or reflected by a body part. Spectroscopic techniques are used to analyze the amount of radiation absorbed at each wavelength by the body part constituents and to compare the absorption data to known data for glucose. Practical implementation of a glucose sensor based on these principles is very difficult and several wavelengths are required. Infrared (IR) spectra are sensitive to physical and chemical factors such as temperature, pH, and scattering. Furthermore, spectroscopy is affected by skin pigmentation, use of medications that absorb various IR wavelengths, alterations in blood levels of hemoglobin or other proteins that absorb IR, changes in body temperature, and alterations in the state of hydration or nutrition. In addition, the NIR spectrum of glucose is very similar to that of other sugars, including fructose, which is often used by diabetics. Therefore, the signal (i.e. the change in the absorption spectrum as a function of glucose concentration) is very small compared to noise and to interference resulting especially from the water spectral absorption and other strongly absorbing components.

[0009] Another approach is Raman Spectroscopy. With Raman spectroscopy, Raman spectra are observed when incident radiation is inelastically scattered. The loss or gain of photon energy are independent of the excitation frequency and provide specific information about the chemical structure of the sample. The Raman signal is very weak, requiring long data acquisition time, making the device sensitive to light source fluctuations. Measurements are subject to high background noise because of tissue autofluorescence. Scatter and reabsorption in biological tissues make detection of Raman frequency shifts due to physiological concentrations difficult.

[0010] Another spectroscopic approach is based on photoacoustics. In photoacoustic spectroscopy, a laser beam pulse is used to rapidly heat the tissue and generate an acoustic pressure wave that can be measured by a microphone or other transducer. The acoustic signal is analyzed to infer blood glucose concentration. Measurements are affected by chemical interferences from biological molecules as well as physical interference from temperature and pressure changes. Current instruments are complex and sensitive to environmental conditions.

[0011] Another optical approach considered of glucose monitoring is based on employing polarimetry. Glucose concentration changes the polarization of light fields. The eye's aqueous humor has been suggested as the medium for this technique as skin is not a feasible site due to its high light scattering properties. However, polarization measurements are affected by optical rotation due to cornea, and by other optically active substances. Other interfering factors include saccadic motion and corneal birefringence. In addition, there is a significant lag between blood glucose changes and glucose changes in intra-ocular fluids, of up to 30 minutes.

[0012] Yet, another approach employed for glucose monitoring is based on light scattering. Changes in glucose levels induce changes in light scattering properties, generally, of the skin. U.S. Pat. No. 6,226,089 to Hakamata discloses detecting the intensities of backscattering light generated by predetermined interfaces of an eyeball when a laser beam emitted from a semiconductor laser is projected onto the eyeball in a predetermined position. The absorbance or refractive index of the aqueous humor in the anterior chamber of the eyeball is determined on the basis of the intensities of the backscattering light, and the glucose concentration in the aqueous humor is determined on the basis of the absorbance or refractive index in the aqueous humor. Light scattering effects are evident in the near-infrared range, where water absorption is much weaker than at larger wavelengths (medium- and far-infrared). However, techniques that rely on the backscattered light from the aqueous humor of the eye are affected by optical rotation due to cornea, and by other optically active substances. Other interfering factors include saccadic motion and corneal birefringence. Finally, it should be appreciated that there is often a significant time lag, (e.g., up to 30 minutes) between blood glucose changes and glucose changes of the intra-ocular fluids.

[0013] Low-Coherence Interferometry (LCI) is one technique for analyzing skin light scattering properties. Low Coherence Interferometry (LCI) is an optical technique that allows for accurate, analysis of the scattering properties of heterogeneous optical media such as biological tissue. In LCI, light from a broad bandwidth light source is first split into sample and reference light beams which are both retro-reflected, from a targeted region of the sample and from a reference mirror, respectively, and are subsequently recombined to generate an interference signal. Constructive interference between the sample and reference beams occurs only if the optical path difference between them is less than the coherence length of the source.

[0014] U.S. Pat. No. 5,710,630 to Essenpreis et al. describes a glucose measuring apparatus for the analytical determination of the glucose concentration in a biological sample and comprising a light source to generate the measuring light, light irradiation means comprising a light aperture by means of which the measuring light is irradiated into the biological sample through a boundary surface thereof, a primary-side measuring light path from the light source to the boundary surface, light receiving means for the measuring light emerging from a sample boundary surface following interaction with said sample, and a secondary-side sample light path linking the boundary surface where the measuring light emerges from the sample with a photodetector. The apparatus being characterized in that the light source and the photodetector are connected by a reference light path of defined optical length and in that an optic coupler is inserted into the secondary-side measurement light path which combines the secondary-side measuring light path with the reference light path in such manner that they impinge on the photodetector at the same location thereby generating an interference signal. A glucose concentration is determined utilizing the optical path length of the secondary-side measuring light path inside the sample derived from the interference signal.

Brief Summary

[0015] The above-discussed and other drawbacks and deficiencies of the prior art are overcome or alleviated by the measurement system and methodology disclosed herein. Disclosed herein in an exemplary embodiment is a method for determining a characteristic of an analyte in a biological sample, the method comprising: directing broadband light by means of a sensing light path at the biological sample, at a target depth defined by the sensing light path and a reference light path; and receiving the broadband light reflected from the biological sample by means of the sensing light path. The method also includes: directing the broadband light by means of the reference light path at a reflecting device; receiving the broadband light reflected from the reflecting device by means of the reference light path; and interfering the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device. The method further includes: detecting the broadband light resulting from interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device; and modulating an effective light path length of at least one of the reference light path and the sensing light path to enhance interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device. The method continues with: determining a magnitude of change of the effective light path length introduced by the modulating when interference is enhanced; determining a variation in an index of refraction from a ratio of the magnitude of change of the effective light path length and the target depth; and determining the characteristic of the analyte in the biological sample from the variation in the index of refraction.

[0016] Also disclosed herein in an exemplary embodiment is a system for determining a characteristic of an analyte in a biological sample, the system comprising: a broadband light source for providing a broadband light; a sensing light path receptive to the broadband light from the broadband light source, the sensing light path configured to direct the broadband light at the biological sample and to receive the broadband light reflected from the biological sample; and a reflecting device. The system also includes: a reference light path receptive to the broadband light from the broadband light source. The reference light path is configured to direct the broadband light at the reflecting device and to receive the broadband light reflected from the reflecting device, the reference light path coupled with the sensing light path to facilitate interference of the broadband light reflected from the biological sample and the broadband light reflected from the fixed reflecting device. The reference light path and the sensing light path cooperating to define a target depth. The system further includes: a detector receptive to the broadband light resulting from interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device; and a modulator associated with at least one of the reference light path and the sensing light path; the modulator, modulating an effective light path length of the at least one of the reference light path and the sensing light path to enhance interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device; and a processing means configured to determine a magnitude of change of the effective light path length introduced by the modulator when interference is enhanced, a variation in an index of refraction from a ratio of the magnitude of change of the effective light path length and the target depth, and the characteristic of the analyte in the biological sample from the variation in the index of refraction.

[0017] Further disclosed herein in another exemplary embodiment is a method for determining a characteristic of an analyte in a biological sample, the method comprising: directing broadband light by means of a sensing light path at the biological sample, at a target depth defined by the sensing light path and a reference light path; receiving the broadband light reflected from the biological sample by means of the sensing light path; directing the broadband light by means of the reference light path at a reflecting device; and receiving the broadband light reflected from the reflecting device by means of the reference light path. The method also includes: interfering the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device; detecting the broadband light resulting from interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device; and modulating an effective light path length of at least one of the reference light path and the sensing light path in accordance with a modulator function having a phase component. The method further includes: determining a magnitude of change of the effective light path length introduced by the modulating when the phase component of the modulator function is generally null; determining a variation in an index of refraction from the magnitude of change of the effective light path length and the target depth; and determining the characteristic of the analyte in the biological sample from the variation in the index of refraction.

[0018] Further disclosed herein in yet another exemplary embodiment is a system for determining a characteristic of an analyte in a biological sample, the system comprising: a broadband light source for providing a broadband light; and a sensing light path receptive to the broadband light from the broadband light source. The sensing light path is configured to direct the broadband light at the biological sample and to receive the broadband light reflected from the biological sample. The system also includes a reflecting device and a reference light path receptive to the broadband light from the broadband light source. The reference light path is configured to direct the broadband light at the reflecting device and to receive the broadband light reflected from the reflecting device. The reference light path is coupled with the sensing light path to facilitate interference of the broadband light reflected from the biological sample and the broadband light reflected from the fixed reflecting device. The reference light path and the sensing light path cooperating to define a target depth. The system further includes: a detector receptive to the broadband light resulting from interference of the broadband light reflected from the biological sample and the broadband light reflected from the reflecting device; and a modulator associated with at least one of the reference light path and the sensing light path. The modulator, modulating an effective light path length of at least one of the reference light path and the sensing light path in accordance with a modulator function having a phase component. The system further includes a processing means configured to determine: a magnitude of change of the effective light path length introduced by the modulator when the phase component of the modulator function is generally null; a variation in an index of refraction from the magnitude of change of the effective light path length and the target depth; and the characteristic of the analyte in the biological sample from the variation in the index of refraction.

Brief Description of the Drawings

[0019] These and other features and advantages of the present invention may be best understood by reading the accompanying detailed description of the exemplary embodiments while referring to the accompanying figures wherein like elements are numbered alike in the several figures in which:

[0020] FIG. 1 is a basic all-fiber low-coherence interferometer (LCI);

[0021] FIG. 2 depicts a plot of the envelope function G(Δl) and of the interference signal G(Δl)cos φ_s;

[0022] FIG. 3 depicts a range of unambiguous measurement for a periodic interference signal;

[0023] FIG. 4A depicts a minimum configuration interferometer system in accordance with an exemplary embodiment of the invention;

[0024] FIG. 4B depicts a configuration of an interferometer system in accordance with an exemplary embodiment of the invention;

[0025] FIG. 5 depicts a plot of the interference signal for a single ramp;

[0026] FIG. 6 depicts an a simplified block depicting a detection scheme employing ramp modulation;

[0027] FIG. 7A depicts an interference signal corresponding to Δn=0 at a 1 mm depth;

[0028] FIG. 7B depicts an interference signal corresponding to Δn=0.0003 at a 1 mm depth;

[0029] FIG. 8A shows a balanced interferometric signal resultant from sinusoidal modulation depicting a heavy portion oscillating at the frequency f over the peak;

[0030] FIG. 8B shows an unbalanced interferometric signal, and a shifted modulation pattern;

[0031] FIG. 9 depicts an implementation of a sinusoidal modulation or homodyne detection scheme;

[0032] FIG. 10 depicts an adaptation of the interferometer system of FIGS. 4A and 4B with a calibration strip;

[0033] FIG. 11A depicts an interface for extension modules in accordance with another exemplary embodiment of the invention;

[0034] FIG. 11B depicts an interface for extension in accordance with another exemplary embodiment of the invention;

[0035] FIG. 11C depicts another interface for extension in accordance with yet another exemplary embodiment of the invention;

[0036] FIG. 12 depicts an adaptation of the interferometer system of FIGS. 4A and 4B for ranging measurements in accordance with another exemplary embodiment;

[0037] FIG. 13 depicts another adaptation of the interferometer system of FIGS. 4A and 4B for ranging measurements in accordance with yet another exemplary embodiment with external probe;

[0038] FIG. 14 depicts an illustration of a splitter-modulator module in accordance with an exemplary embodiment;

[0039] FIG. 15A depicts a process for fabricating the splitter-modulator module in accordance with an exemplary embodiment;

[0040] FIG. 15B depicts a process of fabricating the splitter-modulator module in accordance with an exemplary embodiment;

[0041] FIG. 15C depicts a process of fabricating the splitter-modulator module in accordance with an exemplary embodiment;

[0042] FIG. 16 depicts a miniaturized, handheld LCI system in accordance with an exemplary embodiment; and

[0043] FIG. 17A depicts operation of a miniaturized, handheld LCI system in accordance with an exemplary embodiment; and

[0044] FIG. 17B depicts operation of a miniaturized, handheld LCI system in accordance with another exemplary embodiment.

Detailed Description of an Exemplary Embodiment

[0045] Disclosed herein, in several exemplary embodiments are high-sensitivity low coherence interferometric (LCI) systems (instruments) for optical metrology, which in an exemplary embodiment are miniaturized for use in a variety of sensing and monitoring applications, including, but not limited to, trace chemical sensing, optical properties and changes thereof, medical sensing such as non-invasive glucose monitoring and others. In an exemplary embodiment, the instrument is miniaturized, using integrated optics components such as waveguides, splitters and modulators on a single substrate such as, but not limited to, a LiNbO3 (Lithium Niobate) chip. The exemplary embodiments may also involve the use of a “circulator” type of optical component, including of a polarizing beam splitter and quarterwave plate, which can be combined with the light source and detector into a miniature module that prevents optical feedback into the light source while doubling the detected light. Alternatively, instead of the polarizing beam splitter and quarter wave plate one or more isolators and a waveguide coupler may be employed in a similar module to accomplish the same purpose. Disclosed herein in the exemplary embodiments are multiple methodologies and associated systems employed to derive information from the magnitude and/or phase of an interferometric signal.

[0046] It will be appreciated that while the exemplary embodiments described herein are suitable for the analysis in comparatively highly scattering, i.e. optically heterogeneous biological samples, optically homogeneous (that is, low-scattering or entirely non-scattering) samples also may be analyzed provided suitable implementations of the embodiments of the invention are employed. It may be further appreciated that the methods discussed herein generally do not allow an absolute measurement of the glucose concentration, but rather a relative measurement from a given baseline. Therefore, calibration to establish a baseline is required. For instance, for one exemplary embodiment, a calibration strip is employed to facilitate calibration. Other methodologies, such as using a sample of known index of refraction, or known glucose concentration may also be employed. The particular glucose concentration in the sample may be determined by any previously known procedure, which allows the determination of the absolute glucose concentration.

[0047] It should noted that the light wavelengths discussed below for such methods are in the range of about 300 to about several thousand nanometers (nm), that is in the spectral range from near ultraviolet to near infrared light. In an exemplary embodiment, for the sake of illustration, a wavelength of about 1300 nm is employed. The term “light” as used herein is not to be construed as being limited or restricted to the visible spectral range.

[0048] It is well known that the presence of glucose affects the light scattering properties of tissue, and the refractive index denoted as n of the Interstitial Fluid (ISF) and the refractive index of the scattering centers in the tissue—cell membranes, cellular components and protein aggregates. It is also known that a near infrared (NIR) light of a few milliwatts optical power can penetrate the skin harmlessly, whether being delivered directly from the air, a fiber, some appropriate waveguide, or some combination thereof. In the NIR range the refractive index of ISF is about 1.348–1.352, whereas the refractive index of cellular membranes and protein aggregates ranges from about 1.350 to 1.460. Since the tissue under the skin is highly scattering, the light is scattered in all directions, and only a small amount, the so-called ballistic photons, is captured. It is the light from the ballistic photons captured that are employed to produce an interferometric signal. Raising the glucose concentration in the ISF raises the ISF refractive index by approximately 1.52×10⁻⁶per each mg/dl. Furthermore, the physiological delay of blood glucose transfer from blood to the ISF is of the order of only 2–5 minutes, which makes monitoring of blood glucose in ISF highly practical especially compared to existing methodologies that employ the aqueous humor of the eye.

[0049] It will also be noted that for a homogeneously scattering medium for which a specific property such as the refractive index is to be measured, it is sufficient to probe at a single depth, as the desired information can be obtained from the phase of the interferometric signal, presumed to be independent of the amplitude. In this case, an instrument as described herein can be configured for measurement at a single depth. However, if desired, to probe for inhomogeneities (local changes of absorption, reflection, or refractive index), the instrument may be configured to measure both the amplitude and the phase of the interferometric signal as functions of depth. Described herein in a first exemplary embodiment is a system configured to probe at a fixed depth, while later embodiments may be employed for measurement at variable depths and for general imaging purposes. In any case, emphasis is placed on miniaturization, portability, low power and low cost.

[0050] Finally, it will also be appreciated that while the exemplary embodiments disclosed herein are described with reference and illustration to glucose measurements, applications and implementations for determination of other characteristics of analytes may be understood as being within the scope and breadth of the claims. Furthermore, the methodology and apparatus of several exemplary embodiments are also non-invasive, and thereby eliminate the difficulties associated with existing invasive techniques. In particular, with respect to detection of glucose concentration, the LCI systems of several exemplary embodiments are configured to be non-invasive, avoiding painful lancets and the like.

[0051] Another important consideration is that, as a tool, particularly for medical diagnostic applications, the LCI system of the exemplary embodiments is preferably configured to be easily portable, and for use by outpatients it must be small. Once again, for the purpose of illustration, the LCI system 10 will be described in the context of non-invasive glucose monitoring, capable of detecting the glucose concentration in the dermis by just touching the patient's skin with the instrument. Moreover, the LCI system 10 is configured to be readily hand-held to facilitate convenient measurements by a patient without additional assistance in any location.

[0052] To facilitate appreciation of the various embodiments of the invention reference may be made to FIG. 1, depicting an all-fiber low-coherence interferometer (LCI) system and the mathematical equations developed herein. Referring also to FIG. 4A, in an exemplary embodiment, an LCI system 10 includes, but is not limited to two optical modules: a source-detector module 20a and a splitter-modulator module 40a, and associated processing systems 60. The source-detector module 20a including, but not limited to, a broad-band light source 22, such as a super luminescent diode (SLD) denoted hereinafter as source or SLD, attached to a single-mode fiber 23 or waveguide, an isolator 24 configured to ensure that feedback to the broad band light source 22 is maintained at less than a selected threshold. The source-detector module 20a also includes an optical detector 28.

[0053] The splitter-modulator module 40a includes, but is not limited to, a waveguide input 41, a waveguide output 43, a splitter/coupler 50, and two waveguide light paths: one light path, which is denoted as the reference arm 42, has adjustable length 1r with a reflecting device, hereinafter a mirror 46 at its end; the other light path, which is denoted as the sensing arm 44, allows light to penetrate to a distance z in a medium/object and captures the reflected or scattered light from the medium. It will be appreciated that the captured reflected or scattered light is likely to be only the so-called “ballistic photons”, i.e., those that are along the axis of the waveguide. Provision is also made for one or more modulators 52, 54 in each of the reference arm 42 and sensing arm 44 respectively.

[0054] Continuing with FIG. 4B as well, in another exemplary embodiment, the source-detector module 20b includes, but is not limited to, a polarized broad-band light source 22, attached to a single-mode fiber 23. The source-detector module 20b also includes a polarizing beam splitter 25 with an quarter wave plate 26 employed to ensure a selected polarization configured to facilitate ensuring that feedback to the broad band light source 22 is maintained at less than a selected threshold. The source-detector module 20b also includes an optical detector 28.

[0055] The splitter-modulator module 40b of this embodiment includes, but is not limited to, a waveguide inputs/output 45, a Y-splitter-combiner 51, and the two waveguide arms: reference arm 42, and sensing arm 44. Once again, provision is also made for one or more modulators 52, 54 in each of the reference arm 42 and sensing arm 44 respectively.

[0056] It will be appreciated that while certain components have been described as being in selected modules, e.g., 20, 40, such a configuration is merely illustrative. The various components of the LCI system 10 may readily be distributed in one or more various modules e.g., 20, 40 as suits a given implementation or embodiment. Furthermore, in an exemplary embodiment the waveguide arms 42, 44 and/or fibers 23 are configured for single-transverse-mode transmission, and preferably, but not necessarily, polarization-maintaining waveguides or fibers. Furthermore it will be appreciated that in any of the exemplary embodiments disclosed herein the waveguide and/or fiber tips of each component joined are configured e.g., angled-cleaved in a manner to minimize reflection at the junctions.

[0057] In order to perform the prescribed functions and desired processing, as well as the computations therefore (e.g., the computations associated with detecting and utilizing the interference signal, and the like), the LCI system 10, and more particularly, the processing system 60, may include, but is not limited to a computer system including central processing unit (CPU) 62, display 64, storage 66 and the like. The computer system may include, but not be limited to, a processor(s), computer(s), controller(s), memory, storage, register(s), timing, interrupt(s), communication interface(s), and input/output signal interfaces, and the like, as well as combinations comprising at least one of the foregoing. For example, computer system may include signal input/output for controlling and receiving signals from the source-detector module 20 as described herein. Additional features of a computer system and certain processes executed therein may be disclosed at various points herein.

[0058] The processing performed throughout the LCI system 10 may be distributed in a variety of manners as will also be described at a later point herein. For example, distributing the processing performed in one ore more modules and among other processors employed. In addition, processes and data may be transmitted via a communications interface, media and the like to other processors for remote processing, additional processing, storage, and database generation. Such distribution may eliminate the need for any such component or process as described or vice versa, combining distributed processes in a various computer systems. Each of the elements described herein may have additional functionality that will be described in more detail herein as well as include functionality and processing ancillary to the disclosed embodiments. As used herein, signal connections may physically take any form capable of transferring a signal, including, but not limited to, electrical, optical, or radio.

[0059] The light reflected from the reference mirror 46 (Electric field E_r) in the reference arm 42 and the light reflected or scattered from depth z within the biological sample (Electric field E_s) in the sensing arm 44 are combined at the optical detector 28, whose output current is proportional the combined electric fields. For example, in one instance, the output of the detector is proportional to the squared magnitude of the total electric field E_t=E_r+E_s.

[0060] The detector current I_dis given by:
I_d=η|E_r+E_s|²=I_r+I_s+2√{square root over (I_rI_s)}|G(τ)|cos 2πν_oτ (1)
where η is the detector quantum efficiency (typically <1), I_r=ηE_r·E_r* is the detector current due to E_ralone, I_s=ηE_s·E_s* is the detector current due to E_salone, and the * represents the complex conjugate. E_r·E_r* and E_s·E_s* represent the optical power in the reflected reference field and reflected sensing field, respectively. The quantity τ is the time delay between the reference field E_rand sensing field E_s, and is given by:

[0061]

\begin{matrix} τ = \frac{l_{r}}{c} - \frac{z}{c / n} = \frac{l_{r} - l_{s}}{c} = \frac{Δ l}{c} & (2) \end{matrix}

where l_s=nz and Δl=l_r−l_sand where Δl is the optical path difference between the reference l_rand sensing l_sarms 42, 44, z is the selected or desired target depth in the biological sample, n is the index of refraction in the sample, and c is the speed of light. Also in Equation (1), ν_ois the center frequency of the light source 22, and G(τ) it the cross-correlation function between the reference and sensing fields. Its magnitude is given by:

[0062]

\begin{matrix} \langle G (τ) \rangle = \exp [- {(\frac{πΔ v τ}{2 \sqrt{\ln 2}})}^{2}] & (3) \end{matrix}

where Δν is the FWHM (full width half maximum) frequency bandwidth of the light source 22.

[0063] The last term in Equation (1), the interference term, is the quantity of interest denoted as i_o:
_o(τ)=2√{square root over (I_rI_s)}|G(τ)|cos 2πν_oτ (4)

[0064] It is convenient to express the interference term i_o, in terms of the center wavelength λ_oand the path difference Δl associated with the interferometer, instead of the frequency and time delay. Therefore, using ν_oλ_o=c, where c is the speed of light in vacuum, Δν may be written in terms of the wavelength FWHM bandwidth Δλ, to obtain:

[0065]

\begin{matrix} i_{o} (Δ l) = 2 \sqrt{I_{r} I_{s}} \langle G (Δ l) \rangle \cos ϕ_{s} where ϕ_{s} = \frac{2 π}{λ_{o}} Δ l and & (5) \\ \langle G (Δ l) \rangle = \exp [- {(\frac{Δ l}{L_{c}})}^{2}] & (6) \end{matrix}

where L_cis the coherence length of the light source and is given by

[0066]

\begin{matrix} L_{c} = \frac{2 \sqrt{\ln 2}}{π} \frac{λ_{o}^{2}}{Δλ} = 0.44 \frac{λ_{o}^{2}}{Δλ} . & (7) \end{matrix}

[0067] A plot of the envelope function G(Δl) and if the interference signal G(Δl)cos φ_sis shown in FIGS. 2A and 2B respectively, for an interferometer with a light source 22 having center wavelength λ_o=1.3 μm and FWHM bandwidth Δλ=60 nm (coherence length L_c=12.4 μm). The detected interference signal exhibits a maximum when the interferometer is balanced, i.e., when the path difference Δl=0. As the system 10 becomes increasingly unbalanced, e.g., Δl≠0, the interference signal exhibits maxima and minima of decreasing amplitude over a range determined by Δl.

[0068] It will be appreciated that the interference signal i_oexhibits significant amplitude only over a spatial window of approximately twice the coherence length L_c. As the optical bandwidth increases, the coherence length L_cdecreases and the spatial measurement window narrows. Thus, LCI provides a means for probing samples at precisely defined locations within the samples.

[0069] It is noteworthy to appreciate that the phase, φ_s, of the interference signal i_ochanges by 2π (from a maximum to a minimum then to another maximum) as Δl varies from 0 to λ_o. Therefore, a small change in Δl results in a substantial phase change. It will be further appreciated that the phase of the interference signal i_ois highly sensitive to small changes of optical properties of the mediums, such as refractive indices, or depth z. Thus, while moderate to large changes may readily be observed by measuring the magnitude of the envelope G(Δl), small changes are best detected by measuring the phase φ_sof the interference signal i_o. It will be further appreciated that all the desired information is contained in the range from 0 to 2π. For values of Δl>λ_o, the interference signal i_ois repetitive. Thus, the range from 0 to 2π as indicated in FIG. 3 is a range for which the desired information can be measured without ambiguity. It may also be noted however, that if the coherence length L_cis short enough that the amplitude difference between the main peak and secondary peaks is measurable, then phase measurement beyond 2π may also be realized.

[0070] Therefore, it will be readily be appreciated that there are two types of information, which can be derived from the interference signal i_o: the envelope G(Δl), or its peak G(Δl=0), which may represent scattering, reflection, and absorption; and the more sensitive changes in cos φ_sdue to small optical property changes in the sample. In order to make any such measurements, it is first preferable to separate the DC components I_rand I_sfrom G(Δl) and cos φ_sin the interferometric signal i_odescribed in Equation (5). In one or more exemplary embodiments of the invention, several methodologies are disclosed for extracting the pertinent information from the interference signal i_odescribed in Equation (5). A first methodology addresses a ratio of an effective light path length and a variation thereof for generating the interference signal, while others addresses detecting the phase, or more specifically a particular phase shift φ_sof the interference signal i_ointroduced in the sample.

[0071] Continuing now with FIGS. 4A and 4B, the first approach uses a periodic ramp applied to one of the modulators 52, 54. Another approach, also called a homodyne methodology employs a sine wave applied to one of the modulators 52, 54. It will be appreciated that while for the purposes of description of one or more exemplary embodiments, a particular modulator in a particular arm of the LCI system 10 is described as including a modulator, other configurations are conceivable. For example, while the description of an exemplary embodiment calls for modulation of the length of the reference arm of the LCI system 10, manipulation of other optical lengths in the LCI system 10 may be employed for establishing the interference signal and the level of modulation required to achieve the particular desired result.

[0072] Using one of the modulators, (m₁52, for example) or an equivalent means, a ramp modulation is applied to one of the interferometer arms, the reference arm 42, for example, changing l_rover a distance from −b to a over a time period T, such that:

[0073]

\begin{matrix} Δ l = \frac{a}{T} t - b for t < 0 < T and if periodic, Δ l (t + T) = Δ l (t) . & (8) \end{matrix}

[0074] This yields:

[0075]

\begin{matrix} G (Δ l) \cos (\frac{2 π}{λ_{o}} Δ l) = \exp [- {(\frac{\frac{a}{T} t - b}{L_{c}})}^{2}] \cos [\frac{2 π}{λ_{o}} (\frac{a}{T} t - b)] & (9) \end{matrix}

and,
i_o(t)=2√{square root over (I_rI_s)} cos(2πf_ct−φ_c (10)
where

[0076]

f_{c} = \frac{a}{λ_{o} T} and ϕ_{c} = \frac{2 π b}{λ_{o}} .

[0077] The resultant of the modulation represents a sine wave of frequency f_cwith an arbitrary phase φ_cdetermined by b, which is amplitude-modulated (AM) by the G(Δl) envelope function, now also a function of time. FIG. 5 depicts a plot of the function in equation (9) for a single ramp sweeping over ±2L_cfor the light source example used earlier, with a=4L_cand b=2Lc), and for T=1.9 ms, we get f_c=20 KHz. When the ramp function used for modulation is periodic, this signal repeats with the periodicity of the ramp function. Advantageously, the interference signal may be readily envelope-detected in similar fashion to an AM receiver signal, to obtain G(Δl), and peak-detected, preferably in an iterative fashion, to yield G(Δl=0), which can then be digitized for further processing.

[0078] As discussed earlier, the most sensitive interferometric information can be derived from the phase variations of the interference signal i_oresulting from small changes in material properties. For example, referring to Equation (2), a small change Δn in the refractive index n, at a depth z in the sample, will result in a corresponding change in Δl by the amount zΔn. To analyze such a change by a ramp modulation method as described above, the ramp described above in Equation (8), is applied to the modulator m₁52, but, in addition, an extra optical path length zΔn is also included to facilitate addressing the phase shift in the interference signal i_o, yielding:
Δl=a/Tt−b+zΔn, for 0≦t≦T and Δl(t+T)=Δl(t) (11)
and writing G(Δl)˜1 for a˜λ_o, yields:

[0079]

\begin{matrix} i_{o} (t) = 2 \sqrt{I_{r} I_{s}} \cos (2 π f_{c} t - ϕ_{c} + ϕ_{n}), where ϕ_{n} = \frac{2 π}{λ_{o}} z Δ n . & (12) \end{matrix}

[0080] This interference signal i_ois illustrated in FIGS. 9A and B for a 1 mm depth, where FIG. 7A corresponds to Δn=0, and FIG. 7B corresponds to Δn=0.0003, clearly indicating the phase shift due to Δn. The calculated phase shift φ_nis 1.5 rad (about 80 deg). For glucose concentration monitoring applications such a phase shift would correspond to a glucose concentration of about 197 mg/dl (near the upper limit in humans). Advantageously, this example provides a clear indication of the sensitivity of the phase measurement.

[0081] The phase shift of interference signal i_omay readily be detected using a homodyne process 101 such as depicted in the block diagram of FIG. 6, which produces output signals proportional to the sine and cosine of the phase shift to be measured. In an exemplary embodiment, the change in Δl produced by the ramp modulation generates the interferometric signal i_oin the form of a sine wave, cos(ω_ct+φ_c+φ_n) as in Equation (12), with an additional phase shifting component, denoted φ_nresulting from the change in the refractive index, Δn. At the same time, the parameter a of the ramp modulator 114 is applied to a voltage control oscillator (VCO) 120, configured to generate a sinusoidal function at a frequency ω_cproportional to the magnitude a from the ramp modulation. In an embodiment, a closed loop controller is employed to facilitate determination of φ_nas will be described herein. The sinusoidal function, cos(ω_ct) 122 from the VCO 120 is phase shifted with phase angle φ, as depicted at phase shifting process 124. The phase angle φ is set to π/2, to facilitate formation of quadrature sinusoidal signals. The phase-shifted sinusoid 126 is applied to a first multiplying demodulator 130 configured for the loop closure and denoted as a lock-in amplifier. The sinusoidal function cos(ω_ct) 122 from the VCO 120 (without the phase shift φ) is also applied to a second multiplying demodulator 132. Similar to the earlier embodiments, the interferometric signal i_ofrom the optical detector 28 is amplified and filtered through a narrow band filter 110 around ω_c, to extract the signal cos(ω_ct+φ_c+φ_n), which is then applied to both multipliers 130 and 132 for multiplication by the sinusoidal signal 122 and the phase shifted sinusoidal signal 126. Using the relation cos(x) cos (y)=0.5 cos(x+y)+0.5 cos(x−y), it will be appreciated that the output of the first multiplier 130 contains a term at frequency 2ω_cand a DC term cos(φ_c+φ_n−φ). Advantageously, the DC term cos(φ_c+φ_n−φ) may be readily extracted by low-pass filtering as depicted at filtering process 140, and becomes sin(φ_c+φ_n) if φ=π/2. Similarly, for the second multiplier 132, following similar low-pass filtering at process 142, the dc term is cos(φ_cφ_n). Thus, after filtering the following two dc output signals are obtained:
from the first modulator: I_LP=√{square root over (I_rI_s)} sin(φ_c+φ_n) and (13)
from the lower modulator: I_LP=√{square root over (I_rI_s)} cos(φ_c+φ_n). (14)

[0082] It will be readily appreciated that both signals i.e., sine and cosine, contain the desired information, but the first provides better sensitivity to small Δn changes. Preferably, both signals e.g., Equations 13 and 14, are used with their signs being compared to determine the quadrant in which the phase is located to properly obtain the index change. For example, for phase shifts determined to be between 0 and π/2, both signals are positive and increase or decrease in opposite manner. For phase shifts between π/2 and π, they have opposite signs. For phase shifts between π and 3π/2, both are negative, and for phase shifts between 3π/2 and 2π, they have opposite sign.

[0083] Moreover, while the sine signal is preferably employed because it exhibits greater sensitivity for determining the phase (for example for small phase shifts), it should also be appreciated that advantageously, when the phase is precisely determined, the cosine term provides a precise measure of the magnitude. In other words, with the phase known, the precise peak of the particular period of the cosine term may be ascertained. This approach of determining phase and then the magnitude facilitates determination of the magnitude with greater precision and accuracy than can be achieved with envelope detection methodologies often employed for peak detection.

[0084] Therefore, by measuring the two outputs and observing their signs, the phase and magnitude can be unambiguously ascertained for the range from 0 to 2 π. It is noteworthy to appreciate that this 2π limitation applies only if the variations of G(Δl) [i.e., G(t)] are so small that the difference between the main peak (the absolute peak) and the adjacent ones is not readily measurable. If that difference is measurable, then it is possible to use the above algorithm, together with envelope detection, to ascertain the absolute peak, as well as the secondary peaks, and thus phase shifts in excess of 2π without ambiguity. Moreover, this would be possible if the coherence length L_cis short enough to enable the relative measurements of the peaks and the identification of the absolute peak (see FIG. 3 as an example).

[0085] The index of refraction change Δn can be retrieved in several ways. One methodology involves setting the value of b=0 (and thus φ_c=0) and measuring the values for Equations (13) or (14) directly. Another methodology is to adjust/manipulate the value of b, using a feedback implementation, to result in a change to φ_cuntil the sine function of Equation (13) is nulled. When a null is achieved φ_c=−φ_n, which directly yields:

[0086]

\begin{matrix} Δ n = - \frac{b}{z} & (15) \end{matrix}

where z is the selected target probing depth and b is the magnitude of the change in optical path length introduced to balance the interferometer of the LCI system 10. For example, if the value of φ_crequired for nulling is ascertained to be 5° (0.0873 radians) for λ_o=1.3 μm, it may be ascertained that (from φ=2πb/λ_o) the value of b for nulling to be 0.018 μm, resulting in Δn=1.8×10⁻⁵at z=1 mm. In one exemplary embodiment, a voltage proportional to b is transmitted to an optional VCO 150 to generate a sine wave of frequency proportional to b. That frequency may then be converted to a digital signal by means of an optional zero-crossing detector 152 for further utilization or display.

[0087] The expression of Equation (15) also indicates that greater sensitivity is obtained by probing deeper into the sample. Conversely, due to scattering and absorption in the probed sample e.g., tissue, the magnitude of the returned signal (Electric field E_s) and thus the interferometric signal i_o, is reduced at greater depth in the sample. However, it is also well known that if a signal is periodic, that signal may be retrieved by means of a high-gain narrow-bandwidth amplifier (at the expense of measurement time). Similarly, it can be retrieved by means of autocorrelation, again, at the price of measurement time. Thus, it may be inferred that, by using the phase detection methodology disclosed herein, the measurement sensitivity is, substantially independent of the magnitude of the interferometric signal i_o.

[0088] Turning now to another exemplary embodiment of the invention, in this case, the phase shift of the interference signal i_ois determined utilizing a sinusoidal modulation technique. A small sine wave signal of amplitude δ and frequency f is applied to another modulator 54 (m₂for example in FIGS. 4A and 4B in this instance). It causes Δl to change as δ cos ωt, where ω=2πf, and forms a modulation pattern on the balanced interferometric signal such that it exhibits a “wiggle” or rides back and forth at the frequency f over the peak as shown by the heavy portion of the waveform in FIG. 8A. When the interferometer becomes unbalanced by virtue of a change in a material property, the modulation pattern (the heavy portion of the waveform) shifts resulting in the asymmetry as shown in FIG. 8B by virtue of the optical path change zΔn. Therefore, in an exemplary embodiment, for maintaining a controllable reference phase, provision is made to apply a fixed voltage to modulator m₁52 with a fixed displacement b and resulting in a phase shift φ_cas in the previous case.

[0089] Therefore, the general form of □l at a given depth z with the sinusoidal modulation can be written as:
Δl=zΔn+b+δ cos ωt (16)

[0090] and the general form of φ_sis:

[0091]

\begin{matrix} ϕ_{s} = \frac{2 π}{λ_{o}} z Δ n + \frac{2 π}{λ_{o}} b + \frac{2 π}{λ_{o}} δ \cos ω t = ϕ_{n} + ϕ_{c} + ϕ_{δ} \cos ω t & (17) \end{matrix}

where ω=2 πf. and φ_δ=2πδ/λ_o.

[0092] Taking the cosine of the equation and converting yields:
cos φ_s=cos(φ_n+φ_c)cos(φ_δ cos ωt)−sin(φ_n+φ_c)sin(φ_δ cos ωt). (18)

[0093] The resultant is a phase-modulated signal, which can be expanded into the Fourier-Bessel series, using:

[0094]

\begin{matrix} \cos (\frac{2 π}{λ_{o}} δ \cos ω t) = J_{o} (\frac{2 π}{λ_{o}} δ) + \sum_{n} {(- 1)}^{n} 2 J_{2 n} (\frac{2 π}{λ_{o}} δ) \cos 2 n ω t & (19) \\ \sin (\frac{2 π}{λ_{o}} δ \cos ω t) = \sum_{n} {(- 1)}^{n + 1} 2 J_{2 n - 1} (\frac{2 π}{λ_{o}} δ) \sin (2 n - 1) ω t & (20) \end{matrix}

where J_n(x) is the Bessel function of the first kind of argument x and integer order n.

[0095] Substituting Equations (18), (19) and (20) into Eq (5) for cos φ_s, the interferometric signal now includes a DC term J_o, and an infinite number of harmonics J₁, J₂, J₃, etc. of decreasing magnitude. We keep the lowest terms by filtering through a low-pass filter, a filter centered at f and one centered at 2f, respectively. We obtain the following outputs:
I_oo=2√{square root over (I_rI_s)}J_o(φ_δ)cos(φ_c+φ_n) DC term (21)
i_o1(t)=4√{square root over (I_rI_s)}J₁(φ_δ)sin(φ_c+φ_n)cos ωt Fundamental term (22)
i_o2(t)=4√{square root over (I_rI_s)}J₂(φ_δ)cos(φ_c+φ_n)cos 2ωt Second harmonic term (23)
where G(Δl)≈1, since Δl<<L_c.

[0096] It is noteworthy to appreciate that the DC term I_ooincludes the same information with respect to the argument of interest φ_nas the second harmonic term i_o2(t) and may be employed for ascertaining φ_n. However, it is not as stable as i_o2(t) and is subject to drift and other environmental factors. Furthermore, the DC term is not readily separable from the other detector DC components I_rand I_s. On the other hand, i_o1(t) and i_o2(t) are locked to the oscillator frequency, which can be controlled accurately. Therefore, it is preferred that i_o1(t) and i_o2(t) are selected as the interferometer outputs. The amplitude of the fundamental and second harmonic terms may be ascertained, by various methodologies including, but not limited to, lock-in detection, multiplication and filtering as follows:
I_o1=2√{square root over (I_rI_s)}J₁(φ_δ)sin(φ_c+φ_n) (24)
I₀₂=2√{square root over (I_rI_s)}J₂(φ_δ)cos(φ_c+φ_n). (25)

[0097] FIG. 9 depicts an implementation of a sinusoidal modulation and detection or homodyne methodology shown generally as 200. Elements similar to earlier embodiments are similarly referenced but incremented by one hundred. A source 210 generates a sine wave denoted cos ωt of frequency f and second one of frequency 2f locked to the first one (for example, by means of a frequency doubler). In an exemplary embodiment a frequency f of about 100 KHz is employed, however, it will be evident that other frequencies may readily be utilized depending upon a particular implementation and application of the LCI system 10. In fact, the frequency selected for modulation is not critical, and may be selected in accordance with a given implementation. The frequency selected is only limited by practical design constraints of a selected implementation e.g., size of components, processing requirements to implement functions, and the like. The sine wave at frequency f is fed to modulator m₂54 and to the upper multiplier 230 (e.g., a lock-in amplifier), and the sine wave at frequency 2f is applied only to the lower multiplier 232. A ramp generator 214 similar to that employed in the above embodiments, applies a dc voltage to the m₁modulator 52 to produce the path length bias b, which is adjustable. In an exemplary embodiment, the path length bias b is adjusted to between 0 and λ_oto produce a phase shift φ_cbetween 0 and 2π. The detected interference signal i_ois amplified and passed through narrow-band filters 210 to produce signals at ω and at 2ω, and these two signals are fed to the multipliers 230 and 232 respectively as shown. The output of the multipliers, after low-pass filtering at 240 and 242, are dc signals which correspond to the parameters in Equations (24) and (25).

[0098] It will be readily appreciated that except for the Bessel function coefficients, which are constants whose values are determined by their arguments 2πδ/λ_o, the functions in Equations (24) and (25) are identical to those in Equations (13) and (14) which were obtained from the ramp modulation scheme described earlier. However, advantageously, sinusoidal modulation is more robust and does not need special precautions to generate and control the reference frequency f and the π/2 electronic phase shift as in the previous case. In an exemplary embodiment the selection of the operating frequency f is arbitrary. However, it is preferable to keep it above the shot noise frequency range, i.e., above 120 KHz, where the shot noise is not significant.

[0099] It should also be noted that δ needs not be large in order to maximize J₁. It will be appreciated that the maximum value of J₁(2πδ/λ_o) occurs when (2πδ/λ_o)=1.85, i.e., when δ=0.29λ_o, giving J₁(1.85)=0.582 and J₂(1.85)=0.311. Equations (24) and (25) may then be utilized in a fashion similar to that employed with Equations (13) and (14) in order to measure the phase shift due to the change in index of refraction φ_nwithout ambiguity between 0 to 2π.

[0100] Similarly, the value of Δn may also be obtained from the values in Equations (24) and (25), the known values of J₁, J₂, and the magnitude of the I_sand I_r(which may be obtained via a calibration), and setting/driving φ_cto zero. Also as before, in an exemplary embodiment, an implementation employing feedback may be employed to adjust the bias b until the net phase shift is zero. Thereby, the resultant Δn=−b/z, may readily be obtained. Furthermore, using the value of b to set the frequency of a VCO 150 and zero crossing detector 152, a digital output indicative of the value of b may be obtained to facilitate the computation of Δn.

[0101] It should be further appreciated that the implementation depicted in FIG. 9 includes optional provisions to measure the peak of the interferometric signal i_ounder the application of a ramp to provide information on reflection and scattering coefficients. It also contains provision to apply a signal to an optional additional modulator (such as a fiber coil around a PZT drum) to facilitate ranging measurements by means of extra delays applied to the reference arm 42.

[0102] In another exemplary embodiment, once the interferometric signal i_oand phase thereof is ascertained, for a selected target depth z, additional LCI interferometric signals i_ocorresponding to other target depths may be acquired. Furthermore, if desired, in order to obtain averaged distributions of the LCI signal vs. depth, multiple scans corresponding to multiple target depths may be employed. As disclosed herein, there are several methodologies and exemplary LCI systems that may be employed to acquire an interferometric signal i_ocorresponding to selected depths. In one exemplary embodiment, the modulator m₁52 may be employed to add an additional offsets denoted as Δ to the reference arm 42 corresponding to a group of target depth variations Δz in the vicinity of target depth z. This approach is readily implemented employing the LCI system 10 of FIGS. 4A and 4B. It will be appreciated that the extent of such variations in the target depth Δz are a function of the geometry of the waveguide arm and modulators employed. Additional details regarding the waveguide arms 42, and 44 and the configuration of the modulators m₁52 an M₂54 are addressed at a later point herein. In addition, an extension module (described at a later point herein) may also be implemented to facilitate variations in target depth and depth scans.

[0103] Sensitivity is limited by noise, once the measurement interferometric signal i_ois at or about the noise level it becomes impossible to discern the measurement from noise. Therefore, appreciation of the impact and contributions of noise components becomes instructive with respect to practical considerations of measurement limitations. The detection of the information in Equation (24) is limited by fluctuation noise. The minimum detectable signal is reached when it is equal to the noise, i.e., when the signal-to-noise ratio (SNR) is equal to unity. The noise power is expressed in terms of the photocurrent variance σ_i², which consists of the detector noise σ_r², the photon shot noise σ_s², and for the case of a broadband source 22, the excess photon noise σ_e². Hence, the total noise power is:
σ_i²=σ_r²+σ_s²+σ_e² (26)

[0104] The SNR (for signal I_o1for example), is given by:

[0105]

\begin{matrix} SNR = \frac{I_{o1 (rms)}^{2}}{σ_{i}^{2}}, where I_{01 (rms)} = \frac{I_{o1}}{\sqrt{2}} & (27) \end{matrix}

[0106] The detector noise power is simply the thermal noise due to the input resistance of the receiver. It is given by σ_r²=4kTB, where k is Boltzmann's constant (k=1.38×10⁻²³J/°K), T is the absolute temperature, and B is the bandwidth of the measurement. For a system having 1 KHz system at room temperature (T=300°K), its value is 1.66×10⁻¹⁷W.

[0107] The shot noise, or the noise due to the random arrivals of the photons on the detector 28 from a monochromatic source, obeys Poisson statistics. For a one-ohm resistor, it is given by σ_s²=2eI_dcB, where e is the electronic charge (1.6×10⁻¹⁹coulombs) and I_dcis the average detector DC current. If the total power incident on the detector 28 is of the order of about 1 mW and the detector quantum efficiency is of the order of unity, then I_dcis of the order of 1 mA, and for the same detection bandwidth, the shot noise contribution in a one-ohm input resistor is of the order of about 3.2×10⁻¹⁹W.

[0108] The excess intensity noise from a broadband source 22 is a Bose-Einstein process. It is given by (Rollins) σ_e²=(1+V²)I_dc²B/Δν, where V is the degree of polarization of the light source 22 and Δν its frequency bandwidth. From ν_oλ_o=c, where c is the speed of light in vacuum, the frequency bandwidth is given by cΔλ/λ_o². For a source 22 with single polarization (V=0), center wavelength of 1.3 μm, FWHM wavelength bandwidth of 60 nm (Δν=1.07×10¹³Hz) and the same detector current, input resistance and bandwidth as used previously, the excess intensity noise is about σ_e²=10⁻¹⁶W.

[0109] All the noise components can be reduced by reducing the electrical bandwidth B. However, it is noteworthy to appreciate that, except at very low total optical power, the excess intensity noise is by far the dominant noise for a broadband source, followed by thermal noise, then the shot noise. Thus, for moderate optical power:

[0110]

\begin{matrix} SNR \approx \frac{I_{o1 (rms)}^{2}}{σ_{e}^{2}}, where σ_{e}^{2} = I_{dc}^{2} \frac{B}{Δ v} & (28) \end{matrix}

indicating for a broad band source only the excess intensity noise need be considered to evaluate sensitivity. In addition, analysis of the excess intensity noise permits determination of the advantages of employing a phase based techniques of detection/measurement of the interferometric signal i_o.

[0111] As an example, suppose it is desirable to determine the minimum detectable index change at a depth z in a given specimen with the 1.3-μm light source, 1 KHz detection bandwidth, etc. Assume I_r˜1 mA, I_s˜0.1 mA, J₁=0.582 and let Δn be small enough that the sine term in Equation (24) can be replaced by its argument. Then, for SNR=1, the resultant from Equation (28) yields:
(zΔn)_min≈4×10⁻⁶μm. (29)

[0112] As mentioned earlier, the minimum detectable index change depends inversely on the probing depth. At a depth of 1 mm, i.e., 1000 mm, the minimum index change is Δn_(min)˜4×10⁻⁹. At 2 mm probing depth then Δn_(min)˜2×10⁻⁹. It also depends on the square root of the detection bandwidth, so a factor of three reduction is obtained by reducing the bandwidth from 1 KHz to 100 Hz. The improvement in sensitivity of the phase measurement approach over methodologies employing amplitude stems from the phase factor sin(φ_n) in the expression for the phase signal, which does not exist in the amplitude case. At low signal levels, it multiplies the amplitude signal by the factor 2πzΔn/λ_oin which the quantity 2πz/λ_ois of the order of 6,000 for the examples used in the calculations.

[0113] The interference signal i_orepeats itself after the argument of the sine function in Equation (20) changes by 2π, thus when the coherence length L_cis (larger than about 10 λ_o) such that the maximum detectable phase change without ambiguity is 2π, equivalent to Δl=λ_o. This yields:
(zΔn)_max=λ_o (30)

[0114] Thus, for example, at a 1.3-μm wavelength, Δn_(max)˜1.3×10⁻³at 1-mm probing depth, and Δn_(max)˜6.5×10⁻⁴at 2-mm probing depth. Just as the probing depth can be used to control the sensitivity, it can also be used to control the maximum detectable change. The range of the instrument is the ratio of the maximum to minimum detectable signals. For the examples discussed herein, this range is 3×10⁵or 55 dB. For a larger index change, methodology that may be employed is to use the envelope detection method or a combination of phase and envelope detections as described earlier. Note that, since these results are independent of the overall magnitude of the interferometric signal, it is sufficient to make these measurements at only one depth. The depth needs to be changed only to change the scale of the LCI system 10.

[0115] It is well known that the refractive index change Δn in the dermis due to the presence of glucose is ˜1.52×10⁻⁶per milligram per deciliter (mg/dl). Therefore, by establishing variation in the index of refraction, a glucose determination may be determined. In an exemplary embodiment, a comparison of current measurements for index of refraction and a baseline measurement based on calibration and/or normalization for a particular patient yields an accurate means of determining the current glucose concentration.

[0116] Assuming a linear dependence, this gives a glucose concentration C of
C=6.58×10⁵Δn mg/dl (31)

[0117] With the numbers given above, it will be appreciated that glucose levels from about 0.0026 mg/dl to about 855 mg/dl can be measured at 1-mm probing depth. Similarly, glucose levels from about 0.0013 mg/dl to about 428 mg/dl at 2-mm probing depth. The acceptable glucose concentration in humans ranges from about 70 mg/dl to about 170 mg/dl. Therefore, it is evident that the methodologies disclosed provide a significant benefit when applied particularly to glucose monitoring.

[0118] Correlations of glucose concentrations between LCI-based non-invasive measurements employing such methodologies and invasive measurements in the range of about 80–95% can be achieved. It will further be appreciated that when compared to other analytes, glucose, sodium chloride (NaCl), potassium chloride (KCl), and urea produced the highest changes in ISF refractive index with glucose exhibiting the most pronounced effect.

[0119] Referring once again to FIGS. 4A and 4B, broadband light sources including, but not limited to, SLD's are laser type structures configured and designed to operate substantially without feedback, e.g., of the order of less than 10⁻³, preferably less than 10⁻⁴, more preferably less than 10⁻⁵. In the presence of feedback, the spectrum of the SLD light source 22 may be distorted, the coherence is significantly increased and the spectrum can exhibit very large ripples and even lasing spikes, and thereby may become lasers. Therefore, to prevent distortion and maintain spectral integrity, low coherence, and broadband characteristics, reflections back into the light source 22 are avoided to maintain a broadband light source 22. Thus, in an exemplary embodiment of the LCI system, isolation is provided to alleviate feedback to the light source 22.

[0120] Continuing with FIGS. 4A and 4B, in an exemplary embodiment, the source-detector module 20a, 20b, is configured to prevent the reflected interferometer light from reaching the SLD light source 22 and upsetting its operation. The SLD source 22 is designed and configured such that it is linearly-polarized. SLDs and lasers are “heterostructures” semiconductor devices consisting of a thin “active” layer sandwiched between two “cladding” layers of lower refractive index, all epitaxially grown on a single crystal substrate 23. One such process for fabrication is known as MOCVD (metalorganic chemical vapor deposition). One of the cladding layers is p-doped, and the other is n-doped. The substrate 23 is typically n-doped, and the n-cladding layer is the first to be deposited on it. The structure forms a p-n semiconductor junction diode, in which the active layer is caused to emit light of energy equal to its bandgap upon the application of an electric current.

[0121] The structure is called heterostructure because the active and clad layers are made of different material. This is in contrast with ordinary diodes in which the p-n junction is formulated between similar materials of opposite doping. The use of heterostructure has made it possible to confine the electrical carriers to within the active region, thus providing high efficiency and enabling operation at room temperature. In many heterostructures, light is emitted in both TE polarization (the electric field in the plane of the layer) and TM polarization (electric field perpendicular to the layer).

[0122] However, useful effects are obtained when the active layer is sufficiently thin such that quantum mechanical effects become manifest. Such thin layers are called “quantum well” (QW) layers. Furthermore, the active layer can be “strained”, i.e., a slight mismatch (of about 1%) with respect to the substrate crystal lattice can be introduced during the deposition of the QW layer. The strain can modify the transition characteristics responsible for light emission in beneficial ways. In particular, the light is completely polarized in the TE mode if the strain is compressive. Thus, it is now possible to make a linear polarized laser or broadband SLD by compressive strain of the active layer. In an exemplary embodiment, such a linearly-polarized light source 22 is employed.

[0123] In one exemplary embodiment, as depicted in FIG. 4A, the light from the light source 22 is directed through an isolator 24 configured to transmit light in one direction, while blocking light in the opposite direction. The light is directed to a splitter/coupler 50 of the splitter-modulator module 40a. The source-detector module 20a also contains a detector 28 to receive from the splitter/coupler 50.

[0124] In another exemplary embodiment as depicted in FIG. 4B, the linearly-polarized light from the SLD light source 22 is collimated with lenses 27 and applied to a splitter 25. If a basic 50/50 splitter 24 is employed, half of the returned light goes to the detector 28 and the other half is directed to the SLD light source 22. Once again, in this configuration an isolator 24 may be employed to prevent feedback to the light source 22. Similarly, as stated earlier, in another exemplary embodiment, the splitter 25 is a polarizing beam splitter 25 operating in cooperation with a quarter wave plate 26, employed to prevent feedback light from reaching the light source 22. The polarizing beam splitter 25 facilitates the elimination of feedback to the SLD light source 22 by redirecting substantially all the reflected light from the splitter-modulator module 40b to the detector 28.

[0125] The splitter 25 transmits the horizontally polarized light to the quarter wave plate 26, which coverts the light to another polarization, (for example, circular polarization). Likewise, the returning, circularly polarized light is received by the quarter wave plate 26 and is reconverted to a linear polarization. However, the linear polarization opposite, for example, vertical. The vertically polarized light is transmitted to the polarizing beam splitter 25, which directs all of the light to the detector 28. Advantageously, this approach transmits substantially all of the light i.e., the interference signal, to the detector 28. Whereas embodiments employing the isolator 24 transmits approximately half of the light to the detector 28.

[0126] The polarizing beam splitter 25 is a device that transmits light of one polarization (say the horizontal, or TE-polarized SLD light) and reflects at 90° any light of the other polarization (e.g., vertical or TM-polarized). The quarter-wave plate 26 is a device that converts a linearly polarized incident light to circular polarization and converts the reflected circularly-polarized light to a linearly-polarized of the other polarization which is then reflected at a 90° angle by the polarizing beam splitter 25 to the detector 28. Therefore, essentially all the light transmitted by the light source 22 is re-polarized and transmitted to the splitter-modulator module 40b and all the reflected light from the sample and reflecting device 48 is deflected by the polarizing beam splitter 25 to the detector 28. Advantageously, this doubles the light received at the detector 28 relative to the other embodiments, and at the same time minimizes feedback to the SLD light source 22.

[0127] In an exemplary embodiment an SLD chip for the light source 22 has dimensions of approximately 1 mm×0.5 mm×0.1 mm (length×width×thickness), and emits a broadband light typically of up to 50 mW upon the application of an electric current of the order of 200–300 mA. The light is TE-polarized if the active layer is a compressively strained QW. The FWHM spectrum is of the order of 2% to 3% of the central wavelength emission. A SLD light source 22 with 1.3 μm center wavelength emission and operating at 10 mW output power at room temperature would have a bandwidth of about 40 nm and would require about 200 mA of current. In an exemplary embodiment, for continuous wave (cw) operation at room temperature, the SLD light source 22 may be mounted on an optional thermoelectric cooler (TEC) 32 a few millimeters larger than the SLD light source 22 chip to maintain the temperature of the light source 22 within its specified limits. It will be appreciated that the SLD light source 22 and associated TEC 32 peripherals in continuous operation would have the largest power consumption in the LCI system 10. However, without the TEC 32, the SLD junction temperature would rise by several degrees under the applied current and would operate at reduced efficiency.

[0128] Advantageously, in yet another exemplary embodiment, the utilization of a TEC 70 may readily be avoided without incurring the effects of significant temperature rise by pulsed operation of the SLD light source 22. Pulsed operation has the further advantage of reducing the SLD electrical power requirement by a factor equal to the pulsing duty cycle. Moreover, for selected applications of digital technology and storage, only a single pulse is sufficient to generate an interference signal and retrieve the desired information. Therefore, for example, with pulses of duration 10 μs and 1% duty factor, the LCI system 10 of an exemplary embodiment can average 1000 measurements per second without causing the SLD light source 22 temperature to rise significantly. Thus, for low power consumption, the LCI system 10 should preferably be designed for the SLD light source 22 to operate in a pulsed mode with a low duty cycle and without a TEC 32. In such a configuration the source-detector module 20 would be on the order of about 2 centimeters (cm)×2 cm×1 cm.

[0129] The splitter-modulator module 40a, and 40b of an exemplary embodiment includes a splitter/coupler 50 and Y-splitter/combiner 51 respectively, with a “reference” arm 42 and a “sensing” arm 44, the reference arm 42 having a slightly longer optical path (for example, 1 to 3 mm for measurements in biological tissues) than the sensing arm 44. The optical path difference between the two arms 42, 44 is configured such that the LCI system 10 balanced for the chosen probing depth z. Provision is also made to include a modulator m₁52 and m₂54 in the reference arm 42 and sensing arm 44 respectively.

[0130] In an exemplary embodiment, the splitter/coupler 50, Y-splitter/combiner 51 reference arm 42 and a sensing arm 44 are formed as waveguides in a substrate. However, other configurations are possible, including but not limited to separate components, waveguides, optical fiber, and the like. The substrate 23 for this module should preferably, but not necessarily, be selected such that the waveguides of the arms 42, 44 and modulators 52, 54 can be fabricated on/in it by standard lithographic and evaporation techniques. In one exemplary embodiment, the waveguides of the arms 42, 44 are fabricated by thermal diffusion of titanium or other suitable metal that increases the index of refraction of the substrate, evaporated through masks of appropriate width for single transverse-mode operation. In another exemplary embodiment, the waveguides are formed by annealed proton exchange in an acid bath. This process raises the refractive index in the diffusion region, thus creating a waveguide by virtue of the refractive index contrast between the diffusion region and the surrounding regions. In an exemplary embodiment, is lithium niobate (LiNbO3) is employed as a substrate 23. It will be appreciated that other possible materials, namely ferroelectric crystals, may be utilized such as lithium tantalite (LiTaO3) and possibly indium phosphide depending on configuration and implementation of the LCI system 10.

[0131] Lithium niobate is a ferroelectric crystal material with excellent optical transmission characteristics over a broad wavelength range from the visible to the infrared. It also has a high electro-optic coefficient, i.e., it exhibits a change of refractive index under the application of an external electric field. The refractive index change is proportional to the electric field. The speed of light in a transparent solid is slower than in vacuum because of its refractive index. When light propagates in a waveguide built into the electro-optic material, an applied electric field can alter the delay in the material, and if the electric field is time-varying, this will result in a phase modulation of the light. The LiNbO3 material is very stable, the technology for making it is mature, and LiNbO3 modulators, which can be compact and are commercially available.

[0132] In an exemplary embodiment, the high electro-optic coefficient (refractive index change with applied electric field) of lithium niobate is exploited to facilitate implementation of a modulator, such as modulators m₁52 and m₂54. In this embodiment, a modulator is implemented on or about the waveguide arms 42, 44, by depositing metal electrodes 56, 58 in close proximity to the waveguide arms. In one embodiment, the metal electrodes 56, 58 are deposited on the sides of the waveguide arms 42, 44. In another, the metal electrodes 56, 58 may be deposited on the waveguide arms 42, 44 with an appropriate insulation layer, in a selected region. FIGS. 4A and 4B also shows a diagrammatic depiction of a modulators m₁52, m₂54 in each arm 42, 44 fabricated by depositing metal films (electrodes) 56 on the outside the waveguides and a larger “common” electrode 58 between them. Modulation with modulator m₁52 is obtained by applying a voltage between the upper electrode 56 and the common electrode 58, and modulation with modulator m₂54 is obtained by applying a voltage between the lower 56 and the common electrodes 58. The change of refractive index with applied voltage results in a delay or a change of optical path between for the modulated arm 52, 54. For a given applied voltage, the optical path change depends on the length of the electrodes 56, 58.

[0133] FIG. 14 depicts an illustration of a splitter-modulator module 40b with a Y-splitter 51 and two modulators 52, 54 integrated on a LiNbO3 substrate 23. One method of making the Y-splitter 51 (or splitter/combiner 50 of splitter-modulator module 40a) and waveguide arms 42, 44 is by diffusing titanium or another suitable metal into a substrate 23 at high temperature. Another method of fabrication is by proton exchange in an acid bath. In an exemplary embodiment, titanium and a lithium niobate substrate 23 are employed. The process of fabricating the module 40b (or 40a) is illustrated in FIGS. 15A–C. In the diffusion process, the waveguide pattern is etched in a mask and a thin layer of titanium is vacuum-deposited onto the substrate 23 through the mask. The substrate 23 is then heated in an oven at about 900–1000 degrees C. to diffuse the titanium into the lithium niobate substrate 23. The index of refraction of the diffusion region is slightly higher than that of the surrounding material, and this constitutes waveguides in which light is guided in the diffusion region by virtue of its higher refractive index (just as in an optical fiber where the light propagates in the higher index core). Following diffusion, the metal electrodes 56 and 58 for the modulator(s) 52, 54 are deposited on the sides as shown, with a small spacing d between them. Application of a voltage V between one of the outer electrodes 56 and the negative center electrode 58 establishes an electric field of value V/d across the waveguide e.g. reference arm 42 and/or sensing arm 44. In an exemplary embodiment, the width of the waveguide is approximately 3–5 microns, and the spacing d is only a few more microns wider.

[0134] The refractive index change due to the electro-optic effect is given by

[0135]

\begin{matrix} Δ n = - \frac{1}{2} n_{o}^{3} r \frac{V}{d} & (32) \end{matrix}

where n_ois the refractive index, and r is the electro-optic coefficient. The phase shift of a light of wavelength λ propagating in a LiNbO3 modulator is given by

[0136]

\begin{matrix} Δϕ = π \frac{L}{λ} n_{o}^{3} r \frac{V}{d} & (33) \end{matrix}

where L is the length of the modulator electrodes 56, 58. In the context of the LCI systems 10 disclosed herein, this corresponds to an optical path length change of

[0137]

\begin{matrix} Δ l = \frac{1}{2} n_{o}^{3} rL \frac{V}{d} & (34) \end{matrix}

[0138] Typical material properties are:
r=11.3×10−12 m/V
n_o=2.35

[0139] To obtain larger scale modulations, it will be appreciated that an increase in the voltage on/or the length of the modulator will result in larger changes in the index of refraction by the modulator, resulting in an increased variation of the corresponding phase delay. For example, with a configuration of d=10 microns, an applied voltage of only 3.6 volts is sufficient to yield a value of Δl or b (as discussed above) of 1.3 microns (the wavelength of the light discussed in the examples above). This illustrates that a modulator with a range equivalent to the wavelength λ (for example) 1.3 microns may readily be achieved employing the configuration described.

[0140] In an exemplary embodiment, the reference arm 42 is terminated in an evaporated mirror (metal or quarter-wave stack) 46, and the sensing arm 44 is terminated in an anti-reflection (AR) coating, or is covered with an index-matching agent 48 that prevents or minimizes reflection from the end of the sensing arm 44 when placed in contact with the object to be measured. In such a configuration splitter-modulator module 40 would be on the order of about 2 cm×2 cm×0.5 cm.

[0141] Referring now to FIG. 16, a miniaturized, optionally handheld, LCI system 10 is depicted in accordance with an exemplary embodiment. In an exemplary embodiment, the LCI system 10 is packaged in a small enclosure 12 and includes, but is not limited to, various modules including, but not limited to source-detector module 20a, 20b, splitter-modulator module 40a, 40b and may include one or more additional extension, adapter or interface modules such as 80, 90, and 92 (See FIGS. 4A, 4B and 10–13) or even calibration strip 70. In addition, also optionally packaged within the enclosure may be processing system 60, including processor 62 (not shown in this view) associated controls 63 e.g., keys, selectors, pointers, and the like, display 64, data media 66, as well as communication interfaces 65, and the like as well as rechargeable batteries. Therefore, in one exemplary embodiment the LCI system 10 as packaged in enclosure 12 should be comparable in size to that of a typical cell phone or a Personal Digital Assistant (PDA), i.e., about 4 cm×6 cm×1 cm. to readily facilitate handheld operation.

[0142] Continuing with FIG. 16, it should also be appreciated as mentioned earlier, that various portions of the LCI system 10, and particularly, processing system 60 may be enclosed within the enclosure 12, or associated with an external processing unit 14, or remotely located, such as with a computer processing system 60 in another facility 16. In yet another exemplary embodiment, the LCI system 10 may also include communication interfaces 65, including wireless interfaces (e.g., infrared, radio frequency, or microwave transmitter/receiver) similar to modern computers, cell phones, PDAs, and the like to enable communication, including, but not limited to Internet communication, with external systems 14 and remote facilities 16. For example, as a non-invasive glucose monitor and controller, a sensing portion including the source-detector module 20a, 20b and splitter-modulator module 40a, 40b can be detachable, in the form of a wrist band or wrist watch for continuous monitoring, while the rest of the remainder of the LCI system 10 may be in a patient's pocket, separate computer, at a doctor's office, and the like.

[0143] Referring now to FIGS. 17A and B, to illustrate operation of the LCI system 10, as a non-invasive glucose monitor, the instrument is placed against the skin at some location such as the finger, the back of the hand, or the arm. It should be noted that the finger presents a desirable location for sensing, since it has no hair and few other features that may interfere with the light path. However, other locations are possible. The LCI system 10 would rapidly measure and determine the glucose concentration, (or a multitude of measurements can be made and averaged over a few seconds). A display 66 may also be utilized to provide visual information with respect to measurement to a patient. Furthermore, the LCI system 10 could be coupled to a dispenser embedded in the patient for real-time control and administration of medications such as, insulin.

[0144] The phase associated with a selected length of the reference arm is pre-calibrated to correspond to a set distance (about 1 to 3 mm) under the skin. The spot size for the light at the tip of the sensing fiber or waveguide of the sensing arm 44 is on the order of a few microns. The LCI system 10 may readily be calibrated by placing a strip of known refractive index (or, in the case of a glucose monitor, known glucose content,) and appropriate thickness at the sensing end of the splitter-modulator module 40 prior to performing a measurement. FIG. 10 depicts the LCI system of FIGS. 4A and 4B with a calibration strip in place. The calibration strip 70 can serve the dual purpose of calibration and refractive index matching. Its placement in contact with the splitter-modulator module 40a, 40b does not affect the reference arm 42, since the reference arm light does not penetrate it due to the presence of the end mirror 46. The calibration strip 70 and associated processing may be configured such that the LCI system 10 provides a first reading when the calibration strip 70 is not in contact with the LCI system 10 and a corrected reading when in contact with the calibration strip. Furthermore, the calibration strip may be configured as a disposable item.

[0145] The configuration described above with reference to FIGS. 4A and 4B is convenient to use when the instrument can be placed directly in contact with the sample to provide a reading for a selected depth. Some applications may require the probing depth to be dynamic to enable locating a feature. For example, in medical diagnostics or imaging, the operator may need to probe for features such as tumors, characterized by large changes of optical properties (absorption, reflection, or refractive index change due to a different density). Some other (medical) applications may require a probe to be inserted into the body or object under study. For example, employing an expansion to the embodiments disclosed herein with a fiber probe with a catheter and guide wire to facilitate internal diagnostics and imaging. FIGS. 11A–11C depict an adapter and several expansion or extension modules 90, 92, which can be attached to the LCI system 10 of FIGS. 4A and 4B to provide additional versatility and functionality. FIG. 11A, depicts an adapter 80, configured, in one exemplary embodiment as a short section of waveguides 82, preferably, but not necessarily, made of the same material as the splitter-modulator 40a, 40b, with mirror 46 and AR coating 48, which can be attached to the splitter-modulator 40a, 40b (with matching fluid) to operate as an interface for various extension modules 90, 92. The purpose of the extension module 90 is to provide for adequate lengths of the reference and sensing arms 42, 44 while using a minimum of space, and for adjusting the length of the reference arm 42 and/or sensing arm 44 to enable probing at various depths. The length of the arms 42, 44 can be adjusted in any number of ways, including mechanically changing an air gap between two sections of the reference arm, moving the mirror 46, actually modifying the length of the arm, and the like, as well as combinations including at least one of the foregoing. A preferred way to manipulate the length of an arm 42, 44, in this instance the reference arm 42, in order to maintain small size, accuracy, and stability, is to perform this operation electromechanically.

[0146] Referring now to FIGS. 11B and 11C, in yet another exemplary embodiment, an extension modules 90 and 92 including windings of two lengths of single-mode fibers 94, 96, preferably a polarization maintaining fiber (PMF), (reference and sensing arms respectively) on two drums 98a and 98b. In one embodiment, the drum for the reference arm 42 is made out of a piezoelectric material such as, but not limited to PZT (lead zirconate titanate). The diameter of the drums is selected to be large enough to prevent radiation from the fibers 94, 96 due to the bending for example, about 3–4 centimeters (cm). The diameter of the fibers 94, 96 with claddings is of the order of 0.12 mm. The application of a voltage to the PZT drum 98a causes it to expand or contract, thus straining the reference fiber 94 (for example) and changing its effective length and thereby the optical path length for the reference arm 42. Therefore, as the total length of the unstrained fiber is increased, the total expansion increases as well. For example, if the strain limit for the fiber 94 is about Δl/l is 10⁻⁴, then it requires a 10-meter length of fiber 94 to provide for about a 1 mm extension. Advantageously, a length tens of meters is relatively easy to achieve if the fiber 94 is not too lossy. In the 1.3 μm to 1.55 μm wavelength range, the absorption in optical fibers 94, 96 is of the order of 0.2 dB/Km. There for the losses associated with a 10 meter length would be quite small. Thus, the approach of using a voltage applied a piezoelectric drum e.g., 98a wound with a fiber 94 coil is an effective means to provide changes of several millimeters in the optical path length of the reference arm 42.

[0147] Continuing with FIGS. 11B and 11C, the extension module 90 is configured to provide the extension of the reference and sensing arms 42 and 44 as described above and interfaces with an adapter 80 to facilitate depth profiling. Extension module 92 also includes an evaporated metal mirror 46 to terminate the reference arm 42, while the sensing arm 44 is terminated with a fiber probe 97 configured to facilitate probing such as may include a guidewire and catheter.

[0148] FIGS. 12 and 13 depict various implementations of the extended instrument starting from the base configuration depicted in FIGS. 4A and 4B and using the adapter and the extension modules 80, 90, and 92. FIG. 12 depicts a configuration of an exemplary embodiment where in addition to the source-detector module 20a, 20b and splitter modulator module 40a, 40b and extension module 90 and adapter 80 are employed. This configuration facilitates probing at various depths as well as facilitating depth profile scanning. FIG. 13 depicts a configuration of another exemplary embodiment where in addition to the source-detector module 20 and splitter modulator module 40 and extension module 92 including an external probe 97 are employed. This configuration facilitates probing either at a distance from the device or remote probing such as with a catheter and guidewire. FIG. 12 depicts a configuration of an exemplary embodiment where in addition to the source-detector module and splitter modulator module 40a, 40b and extension module 90 and adapter 80 are employed. This configuration facilitates probing at various depths as well as facilitating depth profile scanning.

[0149] The disclosed invention can be embodied in the form of computer, controller, or processor implemented processes and apparatuses for practicing those processes. The present invention can also be embodied in the form of computer program code containing instructions embodied in tangible media 66 such as floppy diskettes, CD-ROMs, hard drives, memory chips, or any other computer-readable storage medium, wherein, when the computer program code is loaded into and executed by a computer, controller, or processor 62, the computer, controller, or processor 62 becomes an apparatus for practicing the invention. The present invention may also be embodied in the form of computer program code as a data signal 68 for example, whether stored in a storage medium, loaded into and/or executed by a computer, controller, or processor 62 or transmitted over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the computer program code is loaded into and executed by a computer 62, the computer 62 becomes an apparatus for practicing the invention. When implemented on a general-purpose processor the computer program code segments configure the processor to create specific logic circuits.

[0150] It will be appreciated that the use of first and second or other similar nomenclature for denoting similar items is not intended to specify or imply any particular order unless otherwise stated.

[0151] While the invention has been described with reference to an exemplary embodiment, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential scope thereof. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed as the best mode contemplated for carrying out this invention, but that the invention will include all embodiments falling within the scope of the appended claims.

Claims

What is claimed is:

1. A method for diagnosing a human condition from a characteristic of an analyte in a biological sample, the method comprising:

directing broadband light by means of a sensing light path at the biological sample, at a target depth defined by said sensing light path and a reference light path;

receiving said broadband light reflected from the biological sample by means of said sensing light path;

directing said broadband light by means of said reference light path at a reflecting device;

receiving said broadband light reflected from said reflecting device by means of said reference light path;

modulating at about a distance of a center wavelength of said broadband light an effective light path length of at least one of said reference light path and said sensing light path with respect to said target depth;

detecting said broadband light resulting from said broadband light reflected from said target depth in the biological sample and said broadband light reflected from said reflecting device, to generate a signal indicative of an interference of said broadband light reflected from said target depth in the biological sample and said broadband light reflected from said reflecting device, wherein said signal is a product of a cosine function;

determining a phase component that is an argument of the cosine function;

determining a change in the phase component angle of said signal resulting from a change of said effective light path length that is itself induced by a change of a property of the biological sample;

determining a variation in an index of refraction from said change in the phase component of said signal;

determining the characteristic of the analyte in the biological sample from said variation in said index of refraction; and

diagnosing a human condition from the characteristic.

2. The method of claim 1 further comprising varying an effective light path length of at least one of said reference light path and said sensing light path to define an other target depth, said other target depth is at a distance of about at least a coherence length of broadband light from said target depth.

3. The method of claim 2 wherein said varying an effective light path length of at least one of said reference light path and said sensing light path employs at least one of a modulator comprising metallic electrodes disposed at an optical waveguide and a modulator formed with an optical fiber wound on a piezoelectric drum forming at least a portion of at least one of said reference light path and said sensing light pat.

4. The method of claim 1 wherein said modulating comprises sinusoidal modulating.

5. The method of claim 1 wherein said modulating comprises ramp modulating.

6. The method of claim 1 wherein said target depth is defined as a difference between an effective light path length of said reference light path and an effective light path length of said sensing light path.

7. The method of claim 1 wherein said reflecting device comprises a fixed reflecting device.

8. The method of claim 1 wherein at least one of said reference light path and said sensing light path comprises at least one of an optical fiber and a waveguide.

9. The method of claim 1 wherein said modulating said effective light path length is achieved with at least one of a modulator comprising metallic electrodes disposed at an optical waveguide and a modulator formed with an optical fiber wound on a piezoelectric drum forming at least a portion of at least one of said reference light path and said sensing light path.

10. The method of claim 1:

wherein said modulating includes a phase parameter; and

further comprising relating said change in said phase component to said phase parameter, wherein said phase parameter is indicative of said variation in said index of refraction.

11. The method of claim 10 wherein said phase parameter is indicative of said variation in said index of refraction as a ratio of said phase parameter to said target depth.

12. The method of claim 1 further including calibrating at least one of said reference light path and said sensing light path by adjusting effective light path length of at least one of said reference light path and said sensing light path based on a sample exhibiting properties including at least one of known refractive index and scattering coefficient.

13. The method of claim 1 wherein the characteristic of the analyte in the biological sample includes glucose concentration.

14. A system for determining a characteristic of an analyte in a biological sample, the system comprising:

a broadband light source for providing a broadband light;

a sensing light path receptive to said broadband light from said broadband light source, said sensing light path configured to direct said broadband light at the biological sample and to receive said broadband light reflected from the biological sample;

a reflecting device;

a reference light path receptive to said broadband light from said broadband light source, said reference light path configured to direct said broadband light at said reflecting device and to receive said broadband light reflected from said reflecting device, said reference light path and said sensing light path cooperating to define a target depth;

a modulator associated with at least one of said reference light path and said sensing light path, said modulator for modulating at about a distance of a center wavelength of said broadband light an effective light path length of said at least one of said reference light path and said sensing light path;

a detector receptive to said broadband light resulting from said broadband light reflected from the biological sample and said broadband light reflected from said reflecting device, said detector generating a signal indicative of an interference of said broadband light reflected from said target depth in the biological sample and said broadband light reflected from said reflecting device, wherein said signal is a product of a cosine function; and

processing means configured to determine a phase component that is an argument of the cosine function, said processing means further configured to determine a change in the phase component of said signal resulting from a change of said effective light path length that is itself induced by a change of a property of the biological sample, a variation in an index of refraction from said change in the phase component of said signal, and the characteristic of the analyte in the biological sample from said variation in said index of refraction.

15. The system of claim 14 further comprising

means for varying an effective light path length of at least one of said reference light path and said sensing light path to define an other target depth, said other target depth is at a distance of about at least a coherence length of said broadband light from, said target depth.

16. The system of claim 15 wherein said means for varying comprises a modulator comprising metallic electrodes disposed at an optical waveguide and a modulator formed with an optical fiber wound on a piezoelectric drum forming at least a portion of at least one of said reference light path and said sensing light path.

17. The system of claim 14 wherein said modulator comprises a sinusoidal modulator.

18. The system of claim 14 wherein said modulator comprises a ramp modulator.

19. The system of claim 14 wherein said effective light path length of said sensing light path is different from said effective light path length of said sensing light path by a difference, said target depth defined as said difference.

20. The system of claim 14 wherein said reflecting device comprises a fixed reflecting device.

21. The system of claim 14 wherein at least one of said reference light path and said sensing light path comprises at least one of an optical fiber and a waveguide.

22. The system of claim 14 wherein said modulator comprises at least one of metallic electrodes disposed at an optical waveguide and a plezoelectric drum with an optical fiber wound thereon forming at least a portion of at least one of said reference light path and said sensing light path.

23. The system of claim 14 wherein:

said modulator comprises modulating including a phase parameter; and

said processing means further configured to relate said change in said phase component to said phase parameter, wherein said phase parameter is indicative of said variation in said index of refraction.

24. The system of claim 14 further including a calibrating strip for calibrating at least one of said reference light path and said sensing light path by adjusting effective light path length of at least one of said reference light path and said sensing light path based on a sample exhibiting properties including at least one of known refractive index.

25. The system of claim 14 where in the characteristic of the analyte in the biological sample includes glucose concentration.

Patent Citations (124)

Patent	Date	Inventor	Cited By
US4883953(A)	1989-11-01	Koashi	Applicant
US5202745(A)	1993-04-01	Sorin	Applicant
US5321501(A)	1994-06-01	Swanson	Applicant
US5341205(A)	1994-08-01	McLandrich	Applicant
US5383467(A)	1995-01-01	Auer	Applicant
US5398681(A)	1995-03-01	Kupershmidt	Applicant
US5434791(A)	1995-07-01	Koko	Applicant
US5459570(A)	1995-10-01	Swanson	Applicant
US5465147(A)	1995-11-01	Swanson	Applicant
US5491552(A)	1996-02-01	Knuttel	Applicant
US5501226(A)	1996-03-01	Petersen	Applicant
US5507288(A)	1996-04-01	Bocker	Applicant
US5549114(A)	1996-08-01	Petersen	Applicant
US5565986(A)	1996-10-01	Knuttel	Applicant
US5582171(A)	1996-12-01	Chornenky	Applicant
US5710630(A)	1998-01-01	Essenpreis	Applicant
US5830145(A)	1998-11-01	Tenhoff	Applicant
US5835215(A)	1998-11-01	Toida et al.	Applicant
US5835642(A)	1998-11-01	Gelikonov	Applicant
US5836877(A)	1998-11-01	Zavislan	Applicant
US5847827(A)	1998-12-01	Fercher	Applicant
US5867268(A)	1999-02-01	Gelikonov	Applicant
US5877856(A)	1999-03-01	Fercher	Applicant
US5883717(A)	1999-03-01	DiMarzio	Applicant
US5892583(A)	1999-04-01	Li	Applicant
US5905572(A)	1999-05-01	Li	Applicant
US5920390(A)	1999-07-01	Farahi	Applicant
US5921926(A)	1999-07-01	Rolland	Applicant
US5943133(A)	1999-08-01	Zeylikovich	Applicant
US5956355(A)	1999-09-01	Swanson	Applicant
US5962852(A)	1999-10-01	Knuettel	Applicant
US5991697(A)	1999-11-01	Nelson	Applicant
US5994690(A)	1999-11-01	Kulkarni	Applicant
US6002480(A)	1999-12-01	Izatt	Applicant
US6006128(A)	1999-12-01	Izatt	Applicant
US6014214(A)	2000-01-01	Li	Applicant
US6020963(A)	2000-02-01	DiMarzio	Applicant
US6037579(A)	2000-03-01	Chan	Applicant
US6053613(A)	2000-04-01	Wei	Applicant
US6057920(A)	2000-05-01	Fercher	Applicant
US6069698(A)	2000-05-01	Ozawa	Applicant
US6072765(A)	2000-06-01	Rolland	Applicant
US6111645(A)	2000-08-01	Tearney	Applicant
US6124930(A)	2000-09-01	Fercher	Applicant
US6134003(A)	2000-10-01	Tearney	Applicant
US6141577(A)	2000-10-01	Rolland	Applicant
US6144449(A)	2000-11-01	Knuettel	Applicant
US6152875(A)	2000-11-01	Hakamata	Applicant
US6160826(A)	2000-12-01	Swanson et al.	Applicant
US6175669(B1)	2001-01-01	Colston et al.	Applicant
US6191862(B1)	2001-02-01	Swanson et al.	Applicant
US6198540(B1)	2001-03-01	Ueda et al.	Applicant
US6201608(B1)	2001-03-01	Mandella et al.	Applicant
US6208415(B1)	2001-03-01	De Boer et al.	Applicant
US6219055(B1)	2001-04-01	Bhargava	Applicant
US6226089(B1)	2001-05-01	Hakamata	Applicant
US6233055(B1)	2001-05-01	Mandella et al.	Applicant
US6252666(B1)	2001-06-01	Mandella et al.	Applicant
US6268921(B1)	2001-07-01	Seitz et al.	Applicant
US6282011(B1)	2001-08-01	Tearney et al.	Applicant
US6288784(B1)	2001-09-01	Hitzenberger et al.	Applicant
US6304373(B1)	2001-10-01	Zavislan	Applicant
US6307633(B1)	2001-10-01	Mandella et al.	Applicant
US6307634(B2)	2001-10-01	Hitzenberger et al.	Applicant
US6330063(B1)	2001-12-01	Knuettel et al.	Applicant
US6351325(B1)	2002-02-01	Mandella et al.	Applicant
US6370422(B1)	2002-04-01	Richards-Kortum	Applicant
US6377349(B1)	2002-04-01	Fercher	Applicant
US6381015(B1)	2002-04-01	Sonehara et al.	Applicant
US6381025(B1)	2002-04-01	Bornhop	Applicant
US6381490(B1)	2002-04-01	Ostrovsky	Applicant
US6384915(B1)	2002-05-01	Everett et al.	Applicant
US6385358(B1)	2002-05-01	Everett et al.	Applicant
US6390978(B1)	2002-05-01	Irion et al.	Applicant
US6407872(B1)	2002-06-01	Lai et al.	Applicant
US6419360(B1)	2002-07-01	Hauger et al.	Applicant
US6421164(B2)	2002-07-01	Tearney et al.	Applicant
US6423956(B1)	2002-07-01	Mandella et al.	Applicant
US6430455(B1)	2002-08-01	Rebello	Applicant
US6437867(B2)	2002-08-01	Zeylikovich et al.	Applicant
US6441356(B1)	2002-08-01	Mandella et al.	Applicant
US6445939(B1)	2002-09-01	Swanson et al.	Applicant
US6445944(B1)	2002-09-01	Ostrovsky	Applicant
US6456769(B1)	2002-09-01	Furusawa et al.	Applicant
US6466713(B2)	2002-10-01	Everett et al.	Applicant
US6469489(B1)	2002-10-01	Bourquin et al.	Applicant
US6477403(B1)	2002-11-01	Eguchi et al.	Applicant
US6485413(B1)	2002-11-01	Boppart et al.	Applicant
US6496267(B1)	2002-12-01	Takaoka	Applicant
US6498948(B1)	2002-12-01	Ozawa et al.	Applicant
US6501551(B1)	2002-12-01	Tearney et al.	Applicant
US6507747(B1)	2003-01-01	Gowda et al.	Applicant
US6519076(B2)	2003-02-01	Fisher et al.	Applicant
US6522913(B2)	2003-02-01	Panescu et al.	Applicant
US6525862(B2)	2003-02-01	Fisher et al.	Applicant
US6527708(B1)	2003-03-01	Nakamura et al.	Applicant
US6538817(B1)	2003-03-01	Farmer et al.	Applicant
US6546272(B1)	2003-04-01	MacKinnon et al.	Applicant
US6549801(B1)	2003-04-01	Chen et al.	Applicant
US6552796(B2)	2003-04-01	Magnin et al.	Applicant
US6564087(B1)	2003-05-01	Pitris et al.	Applicant
US6564089(B2)	2003-05-01	Izatt et al.	Applicant
US6570659(B2)	2003-05-01	Schmitt	Applicant
US6571117(B1)	2003-05-01	Marbach	Applicant
US6577394(B1)	2003-06-01	Zavislan	Applicant
US6615071(B1)	2003-09-01	Casscells, III	Applicant
US2001/0047137(A1)	2001-11-01	Moreno	Applicant
US2003/0023152(A1)	2003-01-01	Abbink	Applicant
US2003/0023170(A1)	2003-01-01	Gardner et al.	Applicant
US2003/0028100(A1)	2003-02-01	Tearney	Applicant
US2003/0055307(A1)	2003-03-01	Elmaleh	Applicant
US2003/0076508(A1)	2003-04-01	Cornsweet	Applicant
US2003/0137669(A1)	2003-07-01	Rollins	Applicant
US2003/0171691(A1)	2003-09-01	Casscella	Applicant
EP831312(A1)	1998-03-01		Applicant
EP831312(B1)	2001-10-01		Applicant
WO92/19930	1992-04-01		Applicant
WO97/32182	1997-09-01		Applicant
WO9932897(A2)	1997-12-01		Applicant
WO9957507(A1)	1998-04-01		Applicant
WO1120602(A2)	2000-08-01		Applicant
WO2002080767(A1)	2002-03-01		Applicant
WO2003010510(A2)	2002-07-01		Applicant
WO2003088817(A2)	2003-10-01		Applicant

Non-Patent Literature (34)

PCT International Search Report, PCT/US2005/015373, Apr. 5, 2005.Applicant
PCT International Search Report, PCT/US2005/015558, Apr. 5, 2005.Applicant
PCT International Search Report, PCT/US2005/015372, Apr. 5, 2005.Applicant
Optical Coherence tomography—principles and applications, by A.F. Fercher, W. Drexler, C.K. Hitzenberger and T. Lasser, Institute of Physics Publishing, Reports on Progress in Physics, Jan. 20, 2003, pp. 239-303.Applicant
Low Coherence Interferometric Fibre Multiplexed Sensor Systems Using an Integrated-Optical Configuration, by A.J. Rogers, M.J. Plissi, Department of Electronic and Electrical Engineering, King's College London, UK Spie vol. 2510 XP-002112588.Applicant
Integrated Optic Error Detecting Circuit Using Ti:Linbo3 Interferometric Light Modulaters by Hiroshi Haga, Masato Ohta, Masayuki Izutsu and Tadasi Sueta. Faculty of Engineering Science, Osaka University Toyonaka, Osaka 560, Japan IOOOC-ECOC '85.Applicant
Recent progress in fibre optic low-coherence interferometry by Yun-Jiang Rao and David A. Jackson Applied Optics Group, Physics Dept., University of Kent at Canterbury, Kent CT2 7NR, UK Meas. Sci. Technol. 7 (1996) pp. 981-999.Applicant
Joseph M. Schmitt, “Optical Coherence Tomography (OCT): A Review,” IEEE J. Select Topics in Quant. Elect. vol. 5., No. 4, Jul/Aug 1999, pp. 1205-1215.Applicant
Gerard A. Alphonse, “Design of High-Power Superluminescent Diodes with Low Spectral Modulation,” Proceedings of SPIE, vol. 4648, pp. 125-138 (2002).Applicant
Andrew M. Rollins and Joseph A. Izatt, “Optical Interferometer Designs for Optical Coherence Tomography,” Optics Lett. vol. 24, No. 21, Nov. 1, 1999, pp. 1484-1486.Applicant
Rinat O. Esenaliev, Kirill V. Larin and Irina V. Larina, “Noninvasive Monitoring of Glucose Concentration with Optical Coherence Tomography,” Optics Lett. vol. 28, No. 13, Jul. 1, 2001, pp. 992-994.Applicant
Kirill V. Larin, Moshen S. Eledrisi, Massoud Motemedi, and Rinat O. Esenaliev, “Noninvasive Blood Glucose Monitoring with Optical Coherence Tomography—A Pilot Study in the Human Subjects,” Diabetes Care, vol. 25, No. 12, Dec. 2002, pp. 2263-2267.Applicant
H.C. Casey, Jr. and M.B. Panish, Heterostructure Lasers Parts A and B, Academic Press, New York, 1978.Applicant
R.W. Waynant, V.M. Chenault, Overview of Non-Invasive Fluid Glucose Measurement Using Optical Techniques to Maintain Glucose Control in Diabetes Mellitus, LEOS Newsletter, vol. 12, No. 2, Apr. 1998.Applicant
R.J. McNichols, G.L. Cote, Optical Glucose Sensing in Biological Fluids: An Overview, Journal of Biomedical Optics, vol. 5, 2000, pp. 5.Applicant
A. Dogariu and G. Popescu, “High-Resolutionspatial and Spectral Characterization of Optical Fields,” Optics & Photonics News, Dec. 2002, p. 21.Applicant
M. Kohl, M. Cope, M. Essenpries, D. Bocker, “Influence of Glucose Concentration on Light Scattering in Tissue-Simulating Phantoms,” Optics Letters, vol. 19, No. 2, Dec. 15, 1994, pp. 2170.Applicant
J. T. Bruulsema, et al, “Correlation Between Blood Glucose Concentration in Diabetics and Noninvasively Measured Tissue Optical Scattering Coefficient,” Optics Letters, vol. 22, No. 3, Feb. 1, 1997, pp. 190.Applicant
C. Lentner, ed. Geigy Scientific Tables, Ciba-Geigy, Basle, Switzerland, 1984, pp. 69.Applicant
J.M. Schmidt, S.H. Xiang, K.M.Yung, “Speckle in Optical Coherence Tomography” J. of Biomedical Optics, vol. 4, No. 1, Jan. 1999, pp. 95-105.Applicant
A. M. Rollins, J.A. Izatt, Optimal Interferometer Designs for Optical Coherence Tomography, Optics Letters, vol. 24 (21), Nov. 1, 1999, pp. 1484.Applicant
S.L. Jacques, Skin Optics, Oregon Medical Laser Center News, Jan. 1998.Applicant
J.M. Schmitt, G. Kumar, Optical Scattering Properties of Soft Tissue: A Discrete Particle Model, Applied Optics, vol. 37, No. 13, May 1st, 1998, pp. 2788.Applicant
Aristide Dogariu and Gabriel Popsecu, “Measuring the Phase of Spatially Coherent Polychromatic Fields,” School of Optics, University of Central Florida, vol. 89, No. 24, Dec. 9, 2002, pp. 1-2.Applicant
John S. Maier et al., “In Vivo Study of Human Tissues with a Portable Near-Infrared Tissue Spectrometer,” http://lfd.uiuc.edu/spie95/jm/jmspie95.html, pp. 1-8.Applicant
Valery V. Tuchin, “Light Propagation in Tissues with Controlled Optical Properties,” J. Biomedical Optics 2(04), 1997, pp. 401-417.Applicant
Kirill V. Larin, Massoud Motamedi, Taras V. Ashitkov, and Rinat O. Esenaliev, “Specificity of Noninvasive Blood Glucose Sensing Using Optical Coherence Tomography Technique: A Pilot Study,” Physics in Medicine and Biology, Phys. Med. Biol. 48 (2003) pp. 1371-1390.Applicant
Akira Ishimaru, “Diffusion of Light in Turbid Material,” Applied Optics, vol. 28, No. 12, Jun. 15, 1989, pp. 2210-2215.Applicant
Sid Bennett and Steven R. Emge, “Fiber Optic Rate Gyro for Land Navigation and Platform Stabilization,” Sensors Exp 1994, Cleveland, Ohio, Sep. 20, 1994.Applicant
A F. Fercher, K. Mengedoht, and W. Werner, “Eye-Length Measurement By Interferometry with Partially Coherent Light,” Optic Letters, Optical Society of America, vol. 13 No. 3, Mar. 1988, pp. 186-188.Applicant
B.L. Danielson and C.D. Whittenberg, “Guided-wave Reflectometry with Micrometer Resolution,” Applied Optics, vol. 26, No. 14, Jul. 15, 1987.Applicant
Kazumasa Takada, Itaru Yokohama, Kazumori Chida, and Juichi Noda, “New Measurement System for Fault Location in Optical Waveguide Devices Based on an Interferometric Technique,” Applied Optics, vol. 26, No. 9., May 1, 1987, pp. 1603-1606.Applicant
Robert C. Youngquist, Sally Carr and D.E.N. Davies, Optical Coherence-Domain Reflectometry: A New Optical Evaluation Technique, Optics Letters, vol. 12, No. 3, Mar. 1987.Applicant
Harry Delcher, “Continuous Measurement of Glucose in Interstitial Fluid for Extended Time Periods,” SpectRx Inc., pp. 1-8.Applicant

Abstract