Method for non-invasive measurement of an analyte

The invention relates to a method and apparatus for measuring an analyse in blood and other tissue. The method employs Raman spectroscopy wherein a temperature probe within the tissue is excited with electromagnetic radiation. The method results in enhanced spectral information related to concentration of the analyte.

There has long been considerable interest in the non-invasive monitoring of body chemistry. There are 16 million American diabetics, all of whom would benefit from a method for non-invasive measurement of blood glucose levels. Using currently accepted methods for measuring blood glucose levels, many diabetics must give blood five to seven times per day to adequately monitor their insulin requirements. With a non-invasive blood glucose measurement, closer control could be imposed and the continuing damage, impairment and costs caused by diabetes could be minimized.

Blood oximetry is an example of an application of electronic absorption spectroscopy to non-invasive monitoring of the equilibrium between oxygenated and deoxygenated blood (U.S. Patent No. 5,615,673, issued April 1, 1997). Similarly, vibrational spectroscopy is a reliable mode of quantitative and qualitative ex vivo analysis for complex mixtures, and there are reports of in vitro applications of this method to metabolically interesting analyses (S. Y. Wang et al., 1993, Analysis of metabolites in aqueous solution by using laser Raman spectroscopy, Applied Optics 32(6):925-929; A. J. Berger et al., 1996, Rapid, noninvasive concentration measurements of aqueous biological analytes by near-infrared Raman spectroscopy; Applied Optics 35(1):209-212). Infrared measures, such as vibrational absorption spectroscopy, have been applied to skin tissue, but with success limited by unavailability of suitable light sources and detectors at crucial wavelengths, and by heating of the tissue due to the absorption of incident radiation (U.S. Patent No. 5,551,422, see also R. R. Anderson and J. A. Parrish, 1981, The Optics of Human Skin, J. Investigative Dermatology 77(1):13-19). Previous attempts to provide methods for non-invasive blood glucose monitoring are summarized in U.S. Patent No. 5,553,616, issued on September 10, 1996. German Patent No. 19538372 discloses a device for determining blood glucose using Raman spectroscopy of the eye, and a focusing system that directs the spectroscopy at the aqueous humour. European Patent Application No. 0587008 discloses a surface-enhanced Raman spectroscopy immunoassay in which analytes are combined with Raman-active labels in in vitro samples.

One aspect of the invention provides a method of measuring an analyte in a tissue of a subject comprising. <ol id="ol0001" compact="compact" ol-style=""><li>(a) contacting the tissue with electromagnetic radiation having a first excitation wavelength both in the presence of and in the absence of electromagnetic radiation having a second excitation wavelength, wherein the first excitation wavelength is substantially equal to an absorption wavelength of a temperature probe within the tissue able of making resonance Raman scattering occur and the second excitation wavelength is substantially equal to an absorption wavelength of the analyte, able of making resonance Raman scattering occur and wherein the temperature probe and the analyte are sufficiently proximate to one another that energy deposited into one by absorption of radiation is transferred to the other;</li><li>(b) collecting the Raman spectra emitted by the tissue; and</li><li>(c) analyzing the collected spectra by comparing the Raman spectra emitted by the temperature probe in response to the first excitation wavelength in the presence and in the absence of the second excitation wavelength to determine a concentration of analyte present in the tissue.</li></ol>

In one embodiment, the concentration of analyte is determined in accordance with the following relationship: <maths id="math0001" num=""><math display="block"><msub><mi>C</mi><mi>A</mi></msub><mo>=</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="italic">T</mi><mo>/</mo><mfenced separators=""><msub><mi>C</mi><mi>P</mi></msub><mo>⁢</mo><msub><mi>I</mi><mi>A</mi></msub><mo>⁢</mo><msub><mi mathvariant="italic">ε</mi><mi>A</mi></msub><mo>⁢</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="italic">t</mi></mfenced><mo>;</mo></math><img file="EP1037553B1_D0001.tif" /></maths> wherein C_A is the concentration of analyte; wherein ΔT is the temperature shift between the temperature, in Kelvin, associated with Raman scattering emitted by the temperature probe in response to the first and second excitation wavelengths, and the temperature associated with Raman scattering emitted by the temperature probe in response to the first excitation wavelength in the absence of the second excitation wavelength; wherein C_p is the heat capacity of the tissue; wherein I_A is the intensity of the electromagnetic radiation of the second excitation wavelength; wherein ε_A is the absorption coefficient for a molecule of analyte; and wherein Δt is the duration of contact with the second excitation wavelength.

Another aspect of the invention provides a method of measuring an analyte in a blood-containing tissue of a subject comprising: <ol id="ol0002" compact="compact" ol-style=""><li>(a) contacting the tissue with electromagnetic radiation having a first excitation wavelength while the tissue is replete with blood and while the tissue is depleted of blood, wherein the first excitation wavelength is substantially equal to an absorption wavelength of a temperature probe within the tissue able of making resonance Raman scattering occur and wherein the temperature probe and the analyte are sufficiently proximate to one another that energy deposited into one by absorption of radiation is transferred to the other;</li><li>(b) collecting the Raman spectra emitted by the tissue; and</li><li>(c) analyzing the collected spectra by measuring the Raman spectra associated with the analyte and determining the difference between the Raman spectra collected in the blood replete and blood depleted states to determine a concentration of analyte present in the tissue.</li></ol>

Examples of a temperature probe include, but are not limited to, hemoglobin, carboxy-hemoglobin, myoglobin, melanin, and bilirubin. Examples of an analyte include, but are not limited to, glucose, creatinine, pyruvate, drugs, blood gases and urea.

The invention also provides an apparatus for measuring an analyte in a tissue in a subject comprising: <ol id="ol0003" compact="compact" ol-style=""><li>(a) a first light source that emits light having a first excitation wavelength that is substantially equal to an absorption wavelength of a temperature probe within the tissue ;</li><li>(b) a second light source that emits light having a second excitation wavelength that is substantially equal to an absorption wavelength of the analyte able of making resonance Raman scattering occur;</li><li>(c) a detector that is disposed tn receive Raman spectra scattered by the tissue in response to irradiation with the first and second excitation wavelengths and that produces output signals representative of the scatter Raman spectra; and</li><li>(d) a computer coupled to receive the output signals produced by the detector that processes the output signals to derive a value indicative of the concentration of analyte in the tissue from a comparison of output signals representative of the Raman spectra scattered in response to the first and second excitation wavelengths.</li></ol>

BRIEF DESCRIPTION OF THE DRAWINGS

<ul id="ul0001" list-style="none"><li>Figure 1 is a block diagram representing one embodiment of an apparatus of the invention.</li><li>Figure 2 is a block diagram representing a second embodiment of an apparatus of the invention.</li></ul>

DETAILED DESCRIPTION

The invention disclosed herein provides a method and apparatus for measuring an analyte in blood or other tissue of a subject. The method can be performed non-invasively, and provides an improved signal-to-noise ratio over currently available non-invasive techniques. The signal-to-noise ratio is enhanced by techniques, which enhance the signal, reduce the noise, reduce background, and/or unravel spectral congestion. The method achieves the improved signal-to-noise ratio by using Raman spectroscopy and selecting incident wavelengths that excite temperature probes present in tissues of the body. A temperature probe is a component that is sufficiently proximate to an analyte of interest that energy deposited into either the temperature probe or the analyte by absorption of radiation is transferred to the other.

In one embodiment, excitation of the temperature probe is used to achieve resonance enhancement of the Raman spectrum of the temperature probe, which makes the Raman signal more detectable over background fluorescence. In this embodiment, the Raman scattering produced by the temperature probe in response to an excitation wavelength that is absorbed by the temperature probe is measured with and without coincident excitation with a wavelength that is absorbed by the analyte. The shift in Raman scattering created by exciting the analyte is indicative of the amount of analyte present in the tissue.

In another embodiment, interference from background fluorescence is reduced by enhancing the anti-Stokes component of the Raman spectrum of the analyte. In this embodiment, the Raman scattering produced by the analyte is analyzed following excitation with a wavelength that is absorbed by the temperature probe. The excitation of the temperature probe induces a temperature increase that enhances the anti-Stokes component of the Raman spectrum of the analyte.

Raman scattering is the term applied to the phenomenon whereby light that is scattered off a sample of matter is shifted in wavelength from the incident wavelength. The amount of wavelength shift depends on the vibrational motions the matter is capable of undergoing and this provides a sensitive measure of molecular structure. When the incident light is chosen to coincide with a wavelength the molecules of interest absorb, resonance Raman scattering occurs. The interaction between the incident radiation and the molecule is intensified to the extent that the degree of Raman scattering is a thousand to a million times more intense than if the incident light is chosen not to coincide with a molecular absorption. This means that when blue-green light is directed onto a sample of material containing blood (hemoglobin absorbs blue-green light), it scatters 99.9% of the incident photons not absorbed by the sample, without a shift in wavelength of the photons. The other 0.1% of the photons are shifted in wavelength in ways which depend on the structure, local chemical and physical environment of the hemoglobin, vibrational motions of the hemoglobin molecules and on the temperature of the sample.

If another wavelength is used, such as a longer wavelength (e.g., 700 nm or longer) that does not correspond to a wavelength strongly absorbed by hemoglobin or another temperature probe, then the percentage of photons scattered with shifted wavelengths is about a thousand times smaller. That is, typically 99.9999% are scattered without any shift in wavelength when there is no resonance enhancement. Such an enhancement in the strength of Raman scattering is caused by the absorption of the wavelength of the incident light.

This means if a mixture of substances is exposed to an incident wavelength, then those components of the mixture having a resonance enhancement will dominate the production of scattered photons with a shifted wavelength, i.e., the Raman scattering. This is what occurs in blood when blue-green light is used as a probe due to the resonance enhancement of Raman scattering from hemoglobin. The blue to blue-green spectral range is only absorbed preferentially by hemoglobin among blood constituents, and so presents a special opportunity for blood based temperature methods. Hemoglobin is a general temperature probe which, in one example, can be used to detect the temperature rise due to infrared absorption by glucose, such as anhydrous D-glucose (DAG).

The co-localization of glucose and hemoglobin is advantageous. Temperature gradients are inevitable, so the closer the temperature probe to the site of absorption, the stronger the effect on the temperature probe. The temperature is obtained by measuring and comparing the amount of scattered radiation having shorter wavelengths than the incident light, referred to as anti-Stokes scattering, to the amount of light scattered having wavelengths longer than the incident beam, i.e., Stokes scattering. The ratio of the two intensities for a given wavelength shift is directly related to how one calculates an actual temperature. This effect allows one to measure the temperature of the hemoglobin, and the immediate surrounding region, in a completely non-invasive manner. Moreover, it is not necessary to determine the actual temperature. By detecting the change occurring in one or the other of the Stokes or anti-Stokes scattering, one can calibrate the relative change in Raman scattering to a known glucose concentration.

Definitions

All scientific and technical terms used in this application have meanings commonly used in the art unless otherwise specified. As used in this application, the following words or phrases have the meanings specified.

As used herein, "tissue" means any portion of an organ or system of the body, including, but not limited to, skin, capillary beds, blood, muscle, breast and brain.

As used herein, "temperature probe" refers to a component present in the same tissue as the analyte. When energy is deposited into either the temperature probe or the analyte by absorption of radiation, energy is transferred from one to the other.

As used herein, "Raman spectra associated with" a given component refers to those emitted Raman spectra which one skilled in the art would attribute to that component. One can determine which Raman spectra are attributable to a given component by irradiating that component in a relatively pure form, and collecting and analyzing the Raman spectra emitted by the component in the relative absence of other components.

As used herein, "blood replete" refers to a state in which blood flow through a tissue is unobstructed by, for example, vasoconstriction induced by cooling or the application of pressure. The blood replete state can be enhanced by conditions which increase vasodilation, such as warming.

As used herein, "blood depleted" refers to a state in which blood flow through a tissue is substantially restricted and blood volume is minimized. A blood depleted state can be achieved by, for example, cooling and/or applying pressure to the tissue. Pressure can be applied directly using static force, such as by applying a macroscopic object against the tissue (e.g., flat optical window, directed stream of gas), or in a dynamic manner, such as by applying ultrasonic or other frequency acoustic waves to the tissue.

Methods of the Invention

The invention provides a method of measuring an analyte in a tissue of a subject. In one embodiment, the method comprises contacting the tissue with electromagnetic radiation having a first excitation wavelength. The first excitation wavelength is substantially equal to an absorption wavelength of a temperature probe within the tissue. The temperature probe and the analyte are sufficiently proximate to one another that energy deposited into one by absorption of radiation is transferred to the other. Examples of a temperature probe include, but are not limited to, hemoglobin, carboxy-hemoglobin, myoglobin, melanin and bilirubin. Examples of an analyte include, but are not limited to, glucose, urea, creatinine, pyruvate, tyrosine, tryptophan, bicarbonate, electrolytes, lactic acid, drugs, and blood gases such as O₂, CO₂ and NO. In a preferred embodiment, the temperature probe is hemoglobin and the analyte is glucose. The method further comprises collecting the Raman spectra emitted by the tissue, and analyzing the collected spectra to determine the amount of analyte present in the tissue. In one embodiment, the analyzing comprises measuring the Raman spectra associated with the temperature probe.

Optionally, the method can further comprise contacting the tissue with electromagnetic radiation having a second excitation wavelength. The tissue is contacted with the first and second excitation wavelengths simultaneously. The second excitation wavelength is substantially equal to an absorption wavelength of the analyte, and the analyzing comprises comparing the spectra emitted in response to the first and second excitation wavelengths. Information about the analyte can then be derived from the effect of the second excitation wavelength on the emitted spectra as compared to the spectra emitted in response to the first excitation wavelength in the absence of the second excitation wavelength. The Raman spectra analyzed can include both Stokes and anti-Stokes spectra. In one embodiment, the analyzing comprises determining the concentration of analyte from a temperature shift effected in the temperature probe by the second excitation wavelength in accordance with the following relationship: <maths id="math0002" num=""><math display="block"><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">A</mi></msub><mo mathvariant="normal">=</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">T</mi><mo mathvariant="normal">/</mo><mfenced separators=""><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">p</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">I</mi><mi mathvariant="normal">A</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">ε</mi><mi mathvariant="normal">A</mi></msub><mo>⁢</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">t</mi></mfenced><mo mathvariant="normal">;</mo></math><img file="EP1037553B1_D0002.tif" /></maths> wherein C_A is the concentration of analyte; wherein ΔT is the temperature shift between the temperature, in Kelvin, associated with Raman scattering emitted by the temperature probe in response to the first and second excitation wavelengths, and the temperature associated with Raman scattering emitted by the temperature probe in response to the first excitation wavelength in the absence of the second excitation wavelength; wherein C_P is the heat capacity of the tissue; wherein I_A is the intensity of the electromagnetic radiation of the second excitation wavelength; wherein ε_A is the absorption coefficient for a molecule of analyte; and wherein ΔT is the duration of contact with the second excitation wavalength.

In another embodiment, the analyzing comprises measuring the Raman spectra associated with the analyte. The Raman spectra analyzed can comprise anti-Stokes spectra, which are the Raman scattered spectra having wavelengths shorter than the excitation wavelength. The intensity of anti-Stokes scattering is very sensitive to the temperature of the sample. The higher the temperature, the greater the intensity of the anti-Stokes scattered light, at the expense of the Stokes scattered light. Measurements of either or both of the Stokes and anti-Stokes components can be used to obtain a signal that is proportional to the amount of analyte. Anti-Stokes components of the spectra, although typically weaker than Stokes components, are useful because they are less likely than the Stokes components to be obstructed by interfering background. Stokes features can be used, but are preferably avoided when fluorescence or other background radiation is problematic.

Tissue Modulation

In preferred embodiments, the method of the invention is performed in conjunction with tissue modulation. Tissue modulation, as used herein, refers to manipulating the tissue to which the method is applied so that measurements can be made in both blood replete and blood depleted states. The difference between measurements taken in the blood replete and blood depleted states provides a measure indicative of components in the blood while minimizing the effects of extraneous spectroscopic signals due to calluses, dirt, soap residue and other sources associated with the surrounding tissue. When tissue modulation is employed, the analyzing step includes determining the difference between the Raman spectra collected in the blood replete and blood depleted states.

In preferred embodiments, the tissue is blood, such as blood circulating in the capillary bed of the fingertip. Other tissues can be used, such as ear lobe, muscle, skin, breast or brain. The subject is preferably a vertebrate, such as a mammal, bird, reptile or fish. Examples of mammals include, but are not limited to, human, bovine, porcine, ovine, murine, equine, canine, and feline. In a most preferred embodiment, the subject is human.

The method can further include the step of determining blood volume in the tissue of interest. In this embodiment, the analyzing comprises determining the amount of analyte per unit of blood volume. Blood volume can be measured by contacting the tissue modulated region with light chosen to have a wavelength at an isosbestic point for the oxy-deoxyhemoglobin binary equilibrium, i.e. 805 nm or 580 nm. The amount of light that returns is directly related to the volume of blood in the contacted region. I_I is the intensity of isosbestic points wavelength radiation injected into the tissue in question. I_R is the intensity of isosbestic points wavelength radiation returned from the tissue in question. The ratio, I_R/I_I, which normalizes for the variation in blood volume over successive applications of tissue modulation, will often be constant for a series of measurements that occur on the same region of tissue with the same lasers and detectors. Consequently, this factor can be left out in the equations that follow, and measurement of concentration can be in units of blood volume defined by I_R/I_I.

These units of blood volume (or other tissue volume) will also be evident in the units of e, so that the interrogated volume need not appear explicitly in the equations that relate a temperature rise to the intensity of the incident excitation, the concentration of the temperature probe, the heat capacity of the medium and the duration of exposure. The volume will be implicitly included in the numerical value of ε, and can be either measured empirically or estimated from the known extinction coefficients of the relevant substances in <i>in vitro</i> settings. This is possible due to the objective of determining relative changes in analyte concentration. The spatial overlap between the blood volume interrogated region and the analyte interrogated region need only maintain relative constancy over the course of the measurements to provide sufficient precision in the estimate of temperature change and thereby analyte concentration. The regions would be expected to overlap considerably, although the different wavelengths used for different tasks would have different penetration depths. Given other measurement-related constraints, it is not necessary to achieve complete overlap.

In one embodiment, the first excitation wavelength is about 647 to about 400 nm. Preferably, the wavelength is about 550 to about 400 nm. In more preferred embodiments, the wavelength is about 588, 514, 473, 532 or 457 nm. The electromagnetic radiation to be used as the excitation wavelength in the method preferably has a bandwidth of about 0.03 nm to about 0.000003 nm. In one embodiment, the second excitation wavelength is about 0.98, about 1.42, about 1.89, about 2.15, or about 9 to about 11 µm.

In a preferred embodiment, the invention provides a method of measuring blood glucose concentration in a subject. The method comprises contacting tissue of the subject with electromagnetic radiation having a first excitation wavelength in the blue to blue-green range, and while the tissue of the subject is in a blood replete state. In one embodiment, the first excitation wavelength is about 550 to about 400 nm. The method further comprises collecting Raman spectra emitted by the tissue in response to the first excitation wavelength, wherein the Raman spectra comprise Stokes and anti-Stokes spectra. The collected spectra are then analyzed to determine the spectra associated with hemoglobin. The method further comprises contacting tissue of the subject with electromagnetic radiation having a second excitation wavelength of about 0.98, about 1.42, about 1.89, about 2.15, or about 9 to about 11 µm, and while the tissue of the subject is in a blood replete state. The tissue is contacted with the first and second excitation wavelengths simultaneously. The Raman spectra, including Stokes and anti-Stokes spectra, emitted by the tissue in response to the second excitation wavelength are collected and analyzed to determine the spectra associated with glucose. The above steps are also performed while the tissue is in a blood depleted state, either before, after or alternately with, performance of the steps while the tissue is in a blood replete state. The method further comprises determining the net spectra emitted in response to the first and second excitation wavelengths, wherein the net spectra comprise the difference between the spectra obtained in the blood replete state and the blood depleted state, and determining a first ratio of anti-Stokes to Stokes spectra emitted in response to the first excitation wavelength and a second ratio of anti-Stokes to Stokes spectra emitted in response to the second excitation wavelength. The second ratio is divided by the first ratio to obtain a value representative of the concentration of glucose.

In one embodiment, the analyzing comprises determining the concentration of glucose from a temperature shift effected in hemoglobin by the second excitation wavelength in accordance with the following relationship: <maths id="math0003" num=""><math display="block"><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">G</mi></msub><mo mathvariant="normal">=</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">T</mi><mo mathvariant="normal">/</mo><mfenced separators=""><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">p</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">I</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">ε</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">t</mi></mfenced><mo mathvariant="normal">;</mo></math><img file="EP1037553B1_D0003.tif" /></maths> wherein C_G is the concentration of glucose; wherein ΔT is the temperature shift between the temperature, in Kelvin, associated with Raman scattering emitted by hemoglobin in response to the first and second excitation wavelengths, and the temperature associated with Raman scattering emitted by hemoglobin in response to the first excitation wavelength in the absence of the second excitation wavelength; wherein C_P is the heat capacity of the tissue; wherein I_G is the intensity of the electromagnetic radiation of the second excitation wavelength; wherein ε_G is the absorption coefficient for a molecule of glucose; and wherein Δt is the duration of contact with the second excitation wavelength.

In another preferred embodiment, the invention provides a method of measuring blood glucose in a subject. The method comprises contacting tissue of the subject with electromagnetic radiation having a first excitation wavelength in the blue to blue-green range, and while the tissue of the subject is in a blood replete state. In one embodiment, the wavelength is about 550 nm to about 400 nm. The method further comprises collecting Raman spectra, including anti-Stokes spectra emitted by the tissue, and analyzing the collected spectra to determine the spectra associated with glucose. The above steps are also performed while the tissue is in a blood depleted state, either before, after or alternately with, performance of the steps while the tissue is in a blood replete state. The method further comprises determining the net spectra emitted in response to the first excitation wavelength, wherein the net spectra comprise the difference between the spectra obtained in the blood replete state and the blood depleted state, to obtain a value representative of the concentration of glucose. Limiting the analysis to the anti-Stokes portion of the emitted spectra eliminates much of the background found in the Stokes portion of the spectra.

Apparatus

The invention provides an apparatus for measuring an analyte in a sample. The sample can be a tissue in a subject, including, for example, skin, muscle, capillary beds, blood, breast or brain. A diagram representing an apparatus of the invention is shown in Figure 1. The apparatus comprises a first light source 100, a detector 120, and a signal processor 130. The apparatus is suitable for use with the methods of the invention disclosed herein. Suitable components and specific embodiments of the apparatus can be adapted from other Raman spectroscopy systems known in the art (see, e.g., U.S. Patent Nos. 5,553,616; 5,510,894; 5,615,673; and 5,551,422).

The first light source 100 emits electromagnetic radiation having a first excitation wavelength, wherein the first excitation wavelength is substantially equal to an absorption wavelength of a temperature probe within the sample. Preferably, the light source is a laser. Examples of lasers suitable for use in producing the first excitation wavelength include, but are not limited to, external cavity diode lasers, CO₂ lasers and semiconductor lasers. In a preferred embodiment, the apparatus further comprises a second light source 110 that emits electromagnetic radiation having a second excitation wavelength that is substantially equal to an absorption wavelength of the analyte. Examples of lasers suitable for use in producing the second excitation wavelength include, but are not limited to, argon ion, krypton ion and frequency-doubled YAG lasers. In one embodiment, the first excitation wavelength is about 550 to about 400 nm, and the second excitation wavelength is about 0.98, 1.41, 1.89, 2.15, or about 9 to about 11 µm.

The detector 120 is a photosensitive device that is disposed to receive Raman spectra emitted by the sample. The detector 120 can include a wavelength selection device 210. In one embodiment, the wavelength selection device 210 is a wavelength selection spectrograph, preferably utilizing holographic transmission grating and charge coupled device (CCD) array detection. Scattered light can be pre-filtered prior to entering the spectrograph slit using a holographic notch filter. In one embodiment, the filter has a band blocking width of about 250 cm^-1 (or about 1 nm), centered at the corresponding excitation wavelength. In another embodiment, the wavelength selection device 210 is a single channel detector. Examples of a single channel detector include, but are not limited to, a photodiode such as an avalanche photodiode, and a photomultiplier tube. Light entering the single channel detector can be filtered, for example, using a single notch filter or dielectric stack. In preferred embodiments, the detector 120 is cooled using, for example, liquid nitrogen or another suitable cooling method known in the art. The detector 120 produces output signals representative of the Raman spectra scattered by the sample in response to irradiation with the first and/or second light sources 100, 110.

The processor 130 is coupled to receive the output signals produced by the detector 120. In one embodiment, the processor 130 is a computer coupled to receive the output signals produced by the detector. The computer processes the output signals to derive a value indicative of the concentration of analyte in the tissue from a comparison of output signals representative of the Raman spectra scattered in response to the first and second excitation wavelengths. In one embodiment, the computer processes the output signals in accordance with the following relationship: <maths id="math0004" num=""><math display="block"><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">A</mi></msub><mo mathvariant="normal">=</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">T</mi><mo mathvariant="normal">/</mo><mfenced separators=""><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">p</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">I</mi><mi mathvariant="normal">A</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">ε</mi><mi mathvariant="normal">A</mi></msub><mo>⁢</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">t</mi></mfenced><mo mathvariant="normal">;</mo></math><img file="EP1037553B1_D0004.tif" /></maths> wherein C_A, ΔT, C_P, I_A, ε_A and Δt are as defined above.

Another embodiment of the apparatus is shown in Figure 2. Light emitted by the first light source 100 passes through an amplitude modulator 230 to first and second reflectors 240, 250 such that the light joins the path of light emitted by the second light source 110. The light then irradiates the sample 200. Light scattered by the sample reaches a wavelength selection device 210 that includes a filter tuned to transmit a selected Raman feature. One such feature includes wavelengths of 1008 nm to 1029 nm for skin keratin excited with 785 nm. In one embodiment, the wavelength selection device 210 has a 250 µm slit and uses a holographic transmission grating.

The selected signal is then received by a detector 120. In one embodiment, the detector 120 is a charge coupled device (CCD) array, which can be cooled with liquid nitrogen. Alternatively, the detector 120 can comprise a single channel detector. The detector 120 optionally includes spatial characteristics, such as a quadrant spatial layout, allowing simultaneous use of spatial (Hadmaard/other/tissue modulation) encoding/modulation.

The signal then passes to a phase sensitive amplifier 220 that comprises a lock-in amplifier/gated integrator set to demodulate the signal stimulated by modulator 230. The processor 130 then takes analog demodulated signal from the phase sensitive amplifier 220 and performs digitization, storage and data processing. In addition, the processor 130 can provide synchronization with concurrent tissue modulation or spatial encoding.

Examples of temperature probes include, but are not limited to, hemoglobin, carboxy-hemoglobin, myoglobin, melanin and bilirubin. Examples of analytes include, but are not limited to, glucose, urea, creatinine, pyruvate, tyrosine, tryptophan, bicarbonate, electrolytes, lactic acid, drugs and blood gases such as O₂, CO₂ and NO. In one embodiment, the temperature probe is hemoglobin and the analyte is glucose.

EXAMPLES

The following examples are presented to illustrate the present invention and to assist one of ordinary skill in making and using the same. The examples are not intended in any way to otherwise limit the scope of the invention.

Example 1: Use of Resonance Enhancement of Hemoglobin to Obtain an <i>in Vivo</i> Glucose Measurement

For a given molecule, it is possible to identify various wavelengths of light, from the infrared to the ultraviolet parts of the electromagnetic spectrum, which would be preferentially absorbed, compared to the total absorption by all other components of the blood and tissues. Anhydrous d-glucose (DAG) has a molecular extinction coefficient at or near the 9-11 µm wavelength region which is significantly larger than that of most proteins, salts, fats and even water, comprising human tissues and blood. A number of sugars, partially phosphorylated sugars and other intermediates in sugar metabolism absorb at various wavelengths within the same wavelength range.

The concentrations of a majority of these sugars and near sugars are closely correlated with each other by the equilibria associated with sugar metabolism. Individual sugars and related components absorb infrared light at distinctly different wavelengths within and beyond the wavelength range mentioned above. It has not been possible, however, to detect this excitation using conventional means which rely on measurements directed to the same wavelength as that which is incident on the sample. For example, detection of backscattered light at the same wavelength, known as Rayleigh scattered radiation, is not selective since it contains contributions from the blood volume as well as the tissue. Either of these scattering centers can be sufficiently variable to render a nonselective, total power measurement meaningless.

Because absorption of light deposits energy into the illuminated region, it effects a rise in temperature. The amount of temperature rise is proportional to the number of photons absorbed, the heat capacity of the absorber and immediately surrounding region and the relative rates of energy absorption and energy conduction into the surrounding tissues and fluids. By measuring the temperature rise, on exposing a certain blood or tissue sample <i>in vivo</i> to radiation of an appropriate wavelength, one can detect the absorption of the absorbing species. Thus, the measurement of absorption of light, for example through absorption by DAG upon illumination by the correct specific wavelength(s), reduces to the problem of measuring the associated temperature rise. This example describes a simple, versatile method for obtaining non-invasive temperature measurements.

For the specific case of DAG, the method depends on obtaining measurements of Raman scattered radiation from the erythrocytes in which anhydrous d-glucose is localized. The technique can be applied to a variety of analytes and temperature probes. Free DAG can be found in blood in the cytoplasmic volume of the erythrocytes. Exposing the blood volume through the skin via the capillary beds to infrared radiation at the correct wavelength, (a specific wavelength in the 9-10 µm wavelength range for example), so that it will be preferentially absorbed by DAG, will cause the intraerythrocyte materials to increase in temperature. One substance in close proximity to the point of absorption, and thereby heated by the energy deposited into the DAG, is hemoglobin.

Blood volume can be measured by contacting the tissue modulated region with light chosen to have a wavelength at an isosbestic point for the oxy-deoxyhemoglobin binary equilibrium, i.e. 805 nm or 580 nm. The amount of light that returns is directly related to the volume of blood in the contacted region. I_I is the intensity of isosbestic points wavelength radiation injected into the tissue in question. I_R is the intensity of isosbestic points wavelength radiation returned from the tissue in question. The ratio, I_R/I_I, which normalizes for the variation in blood volume over successive applications of tissue modulation, will often be constant for a series of measurements that occur on the same region of tissue with the same lasers and detectors. Consequently, this factor is left out in the equations that follow and measurement of concentration will be in units of blood volume defined by I_R/I_I.

A normalizing factor relating the quality of the optical contact between the instrumental optical system and the tissue region of interest is obtained by contacting the region with a wavelength that is not absorbed by any component of the region, and detecting the amount of the light which returns. I_N is the intensity of non-absorbed wavelength radiation injected into the tissue in question. Is is the intensity of non-absorbed wavelength radiation returned from the tissue in question. It should be pointed out that the ratio, I_S/I_N, which normalizes for the variation in blood volume over successive applications of tissue modulation, will often be constant for a series of measurements that occur on the same region of tissue with the same lasers and detectors. Consequently, this factor is left out in the equations that follow and measurement of concentration is in units of blood volume defined by I_S/I_N.

<u style="single">Step 1:</u> During this step, there is no external pressure applied to the tissue being examined. A Raman spectrum, Stokes and anti-Stokes, is obtained by optically contacting the skin of a finger tip and the tissue/capillary bed below the skin with narrow spectral bandwidth light. The light that scatters due to the interaction of this "exciting" light with the tissue is collected and analyzed by wavelength. The exciting light is chosen so that the intensity of the Raman spectrum of hemoglobin is resonance enhanced.

The Stokes part of the Raman spectrum obtained during step 1 is referred to as R_o(λ)_s1. The anti-Stokes part of the Raman spectrum obtained during step1 is referred to as R_o(λ)_as1.

The Raman scattered light, i.e. wavelength shifted light, has its origin in the interaction of the incident light, i.e. the exciting light, with the various chemical species in the mixture. In addition to any fluorescence the absorption may have induced, part of the Raman spectrum occurs at longer wavelengths than the incident light. This is called Stokes scattering. In contrast to fluorescence following the absorption of a single photon, another part of the Raman scattered light also can be observed at shorter wavelengths. This is called anti-Stokes scattering. The intensity of this scattering to shorter wavelengths is very sensitive to the temperature of the sample. The higher the temperature, the greater the intensity of the anti-Stokes scattered light at the expense of the Stokes scattered photons. The ordinary Raman scattering process, which involves the interaction of the vibrations of the molecules with the incident light, tends to produce one wavelength shifted photon of light for every 10⁶ incident photons. In contrast, ordinary fluorescence produces one photon for every 10³ incident photons.

Temperature

The following equation relates the intensity of a generic Stokes (I_S) and anti-Stokes (I_AS) components of a vibrational mode having frequency (υ), to the temperature (T) of the sy stem, the frequency (υ₀) of the exciting radiation and q, the ratio of the detector response at the two wavelengths corresponding to υ_AS and υ_S. h is Planck's constant and k is Boltzmann's constant. <maths id="math0005" num=""><math display="block"><msub><mi mathvariant="normal">I</mi><mi>AS</mi></msub><mo mathvariant="normal">/</mo><msub><mi mathvariant="normal">I</mi><mi mathvariant="normal">S</mi></msub><mo mathvariant="normal">=</mo><mrow><mo mathvariant="normal">(</mo><msub><mi mathvariant="normal">υ</mi><mi>AS</mi></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">υ</mi><mn mathvariant="normal">0</mn></msub><mo mathvariant="normal">)</mo><mn mathvariant="normal">4</mn><mo mathvariant="normal">/</mo><msup><mfenced separators=""><msub><mi mathvariant="normal">υ</mi><mi mathvariant="normal">S</mi></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">υ</mi><mn mathvariant="normal">0</mn></msub></mfenced><mn mathvariant="normal">4</mn></msup></mrow><mi mathvariant="normal">q exp</mi><mfenced separators=""><mi mathvariant="normal">hυ</mi><mo mathvariant="normal">/</mo><mi>kT</mi></mfenced></math><img file="EP1037553B1_D0005.tif" /></maths>

Glucose has a strong Raman feature at 200 cm^-1. A temperature rise of ≈ 12 °C, which would occur on a transient basis, and is below the static temperature known to cause histological damage, would cause a 10% change in the ratio I_AS/I_S. The transient temperature rise due to absorption of radiation, ΔT, is proportional to the total amount (energy) of radiation absorbed, E_A. and the net heat capacity of the medium, C_P. The heat capacity of human tissue is on the same order as water. Based on an irradiated volume on the order of microliters, absorbing only on the order of 10⁰-10¹ microjoules is sufficient to cause such a temperature rise. This amount of energy can be produced by CO₂ lasers in a bandpass around 9-10 µm in wavelength, which happens to be absorbed by glucose and a number of its immediate metabolic progeny. There are other wavelengths in the near-infrared spectral region, such as 980, 1410, 1890, 2150 nm, that are also thought to be selectively absorbed by glucose and are more conveniently produced using semiconductor lasers. <maths id="math0006" num=""><math display="block"><mi mathvariant="normal">ΔT</mi><mo mathvariant="normal">=</mo><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">E</mi><mi mathvariant="normal">A</mi></msub></math><img file="EP1037553B1_D0006.tif" /></maths>

The energy absorbed, within the irradiated volume, is proportional to the concentration of absorbing species, i.e. glucose molecules, C_G, the intensity of the radiation at the wavelength in question, and the absorption coefficient per such molecule. E_A is expressed as: <maths id="math0007" num=""><math display="block"><msub><mi mathvariant="normal">E</mi><mi mathvariant="normal">A</mi></msub><mo mathvariant="normal">=</mo><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">I</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">ε</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">t</mi></math><img file="EP1037553B1_D0007.tif" /></maths>

Glucose and hemoglobin tend to be spatially localized in the erythrocytes such that energy deposited into one by absorption of radiation will tend to be transferred to the other (and other components of the localized region), preferentially affecting the spectra of each. The energy absorbed is proportional to the intensity of the incident radiation, the duration (Δt) the absorbers are exposed, the concentration of absorbers (glucose, C_G), and the extinction coefficient of the absorbers at the wavelength of the incident radiation, e.g. ε_G for glucose.

Beer's Law applies to the absorption of radiation by tissue in the following relative sense. The effective pathlength for a particular wavelength can be held constant during the spectral interrogation, and the incident radiation can be maintained at a constant intensity. This is true, independent of the fact that light in tissue propagates in the multiple scattering limit (R. R. Anderson and J. A. Parrish, 1981, The Optics of Human Skin, J. Investigative Dermatology 77(1):13-19). In addition, the temporal synchronization of the various exposures and measurements can be held constant so the relative absorption is dependent only on the relative concentration of the absorbers. Thus the relative temperature change of the hemoglobin is dependent on the concentration of glucose.

Resonance Enhancement

There is a range of wavelengths spanning the blue-green part of the visible spectrum, and to a weaker extent even beyond the 600 nm range, which correspond to electronic absorption of the hemoglobin in erythrocytes. Excitation of Raman spectra of hemoglobin, either oxygenated or deoxyhemoglobin, using incident light in these ranges (particularly in the blue-green) results in as much as a 10³-10⁶ -fold enhancement of the Raman intensity. This enhancement is due to the electronic interaction intensifying the scattering process and it makes Raman spectroscopy a particularly sensitive probe of such molecules, particularly in the presence of a fluorescence background.

<u style="single">Step 2:</u> During this step, there is no external pressure applied to the tissue being examined. Another resonance enhanced Raman spectrum, Stokes and anti-Stokes, is obtained by optically contacting the skin of a finger tip with narrow spectral bandwidth light. Simultaneously, the tissue is contacted with a light source having intensity I_G, capable of preferentially being absorbed by glucose over other species in the tissue/blood.

The Stokes part of the Raman spectrum obtained during step 2 is referred to as R_o(λ)_s2, and the anti-Stokes part of the Raman spectrum obtained during step 2 as R_o(λ)_as2.

<u style="single">Step 3:</u> During this is step, external pressure is applied to the tissue being examined. A Raman spectrum, Stokes and anti-Stokes, is obtained by optically contacting the skin of a finger tip and the tissue/capillary bed below the skin with narrow spectral bandwidth light. The light that scatters due to the interaction of this "exciting" light with the tissue is collected and analyzed by wavelength. The exciting light is chosen so that the intensity of the Raman spectrum is resonance enhanced for the temperature probe.

The Stokes part of the Raman spectrum obtained during step 3 is referred to as R_p(λ)_s3, and the anti-Stokes part of the Raman spectrum obtained during step 3 as R_p(λ)_as3.

<u style="single">Step 4</u>: During this step, external pressure is applied to the tissue being examined. Another resonance enhanced Raman spectrum, Stokes and anti-Stokes, is obtained by optically contacting the skin of a finger tip with narrow spectral bandwidth light. Simultaneously, the tissue is contacted with a light source of intensity, I_G, capable of preferentially being absorbed by glucose over other species in the tissue/blood.

The Stokes part of the Raman spectrum obtained during step 4 is referred to as R_p(λ)_s4, and the anti-Stokes part of the Raman spectrum obtained during step 4 as R_p(λ)_as4. Let: <maths id="math0008" num=""><math display="block"><mfenced separators=""><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">1</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">3</mn></mrow></msub><mo mathvariant="normal">)</mo><mo mathvariant="normal">/</mo><mo mathvariant="normal">(</mo><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">1</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">3</mn></mrow></msub></mfenced><mo mathvariant="normal">=</mo><mi mathvariant="normal">A</mi></math><img file="EP1037553B1_D0008.tif" /></maths><maths id="math0009" num=""><math display="block"><mfenced separators=""><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">1</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">3</mn></mrow></msub><mo mathvariant="normal">)</mo><mo mathvariant="normal">/</mo><mo mathvariant="normal">(</mo><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">1</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">3</mn></mrow></msub></mfenced><mo mathvariant="normal">=</mo><msup><mfenced separators=""><msub><mi mathvariant="normal">υ</mi><mi>AS</mi></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">υ</mi><mn mathvariant="normal">0</mn></msub></mfenced><mn>4</mn></msup><mo mathvariant="normal">/</mo><msup><mfenced separators=""><msub><mi mathvariant="normal">υ</mi><mi mathvariant="normal">S</mi></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">υ</mi><mn mathvariant="normal">0</mn></msub></mfenced><mn mathvariant="normal">4</mn></msup><mo>⁢</mo><mi mathvariant="normal">q exp</mi><mfenced separators=""><mi mathvariant="normal">hυ</mi><mo mathvariant="normal">/</mo><msub><mi>kT</mi><mn>0</mn></msub></mfenced><mn>.</mn></math><img file="EP1037553B1_D0009.tif" /></maths>

Where T₀ corresponds to the temperature, in Kelvin, associated with the Raman scatterers without the preferential excitation of glucose.

Similarly, let: <maths id="math0010" num=""><math display="block"><mfenced separators=""><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">2</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">4</mn></mrow></msub><mo mathvariant="normal">)</mo><mo mathvariant="normal">/</mo><mo mathvariant="normal">(</mo><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">2</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">4</mn></mrow></msub></mfenced><mo mathvariant="normal">=</mo><mi mathvariant="normal">B</mi></math><img file="EP1037553B1_D0010.tif" /></maths><maths id="math0011" num=""><math display="block"><mfenced separators=""><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">2</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">4</mn></mrow></msub><mo mathvariant="normal">)</mo><mo mathvariant="normal">/</mo><mo mathvariant="normal">(</mo><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">2</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi mathvariant="normal">s</mi><mo>⁢</mo><mn mathvariant="normal">4</mn></mrow></msub></mfenced><mo mathvariant="normal">=</mo><msup><mfenced separators=""><msub><mi mathvariant="normal">υ</mi><mi>AS</mi></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">υ</mi><mn mathvariant="normal">0</mn></msub></mfenced><mn>4</mn></msup><mo mathvariant="normal">/</mo><msup><mfenced separators=""><msub><mi mathvariant="normal">υ</mi><mi mathvariant="normal">S</mi></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">υ</mi><mn mathvariant="normal">0</mn></msub></mfenced><mn mathvariant="normal">4</mn></msup><mo>⁢</mo><mi mathvariant="normal">q exp</mi><mfenced separators=""><mi mathvariant="normal">hυ</mi><mo mathvariant="normal">/</mo><msub><mi>kT</mi><mn>0</mn></msub></mfenced><mn>.</mn></math><img file="EP1037553B1_D0011.tif" /></maths>

Where T corresponds to the temperature, in Kelvin, associated with the Raman scatterers with the preferential excitation of glucose.

So that: <maths id="math0012" num=""><math display="block"><mi mathvariant="normal">B</mi><mo mathvariant="normal">/</mo><mi mathvariant="normal">A</mi><mo mathvariant="normal">=</mo><mi mathvariant="normal">C exp</mi><mfenced separators=""><mi mathvariant="normal">ΔT</mi><mo mathvariant="normal">/</mo><msub><mi>TT</mi><mn mathvariant="normal">0</mn></msub></mfenced></math><img file="EP1037553B1_D0012.tif" /></maths><maths id="math0013" num=""><math display="block"><mrow><mi>ln</mi><mo>(</mo><mi mathvariant="normal">B</mi><mo mathvariant="normal">/</mo><mi mathvariant="normal">A</mi><mo>)</mo></mrow><mo mathvariant="normal">=</mo><mi>ln</mi><mfenced><mi mathvariant="normal">C</mi></mfenced><mo mathvariant="normal">+</mo><mfenced separators=""><mi mathvariant="normal">ΔT</mi><mo mathvariant="normal">/</mo><msub><mi>TT</mi><mn mathvariant="normal">0</mn></msub></mfenced></math><img file="EP1037553B1_D0013.tif" /></maths>

Thus the measured ratio B/A can be used in the following: <maths id="math0014" num=""><math display="block"><mi mathvariant="normal">ΔT</mi><mo>=</mo><msub><mi>TT</mi><mn mathvariant="normal">0</mn></msub><mrow><mi>ln</mi><mo>(</mo><mi mathvariant="normal">B</mi><mo mathvariant="normal">/</mo><mi mathvariant="normal">A</mi><mo>)</mo></mrow></math><img file="EP1037553B1_D0014.tif" /></maths> with <maths id="math0015" num=""><math display="block"><mi>ΔT</mi><mo>=</mo><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">E</mi><mi mathvariant="normal">A</mi></msub><mo mathvariant="normal">=</mo><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">I</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">ε</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">t</mi></math><img file="EP1037553B1_D0015.tif" /></maths> or <maths id="math0016" num=""><math display="block"><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">T</mi><mo mathvariant="normal">/</mo><mfenced separators=""><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">I</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><msub><mi mathvariant="normal">ε</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><mi mathvariant="normal">Δ</mi>&lt;mtable /&gt;<mi mathvariant="normal">t</mi></mfenced><mo>=</mo><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">G</mi></msub><mn>.</mn></math><img file="EP1037553B1_D0016.tif" /></maths>

It is possible to calculate the concentration of glucose, C_G, directly from the set of measurements, with many interfering influences dividing themselves out along the way.

Measurements of relative temperature rise can be empirically correlated with independently measured glucose concentrations. The above discussion demonstrates an algorithm and shows that there is a rigorous basis for empirically applying two variations of the method for noninvasive analyte quantitation.

Example 2: Use of Anti-Stokes Raman Spectra to Obtain an <i>in Vivo</i> Glucose Measurement

The strength of the anti-Stokes component of the glucose spectrum is dependent on the intensity of exciting radiation, which excites hemoglobin and heats glucose. The intensity of anti-Stokes scattering is very sensitive to the temperature of the sample. The higher the temperature, the greater the intensity of the anti-Stokes scattered light, at the expense of the Stokes scattered light. Anti-Stokes components of the spectra are useful because they are less likely than the Stokes components to be obstructed by interfering background. In this Example, an excitation wavelength that excites hemoglobin is used to heat glucose, thereby enhancing the anti-Stokes spectra associated with glucose.

<u style="single">Step 1:</u> During this step. there is no external pressure applied to the tissue being examined. A Raman spectrum, Stokes and anti-Stokes, is obtained by optically contacting the skin of a finger tip and the tissue/capillary bed below the skin with narrow spectral bandwidth light. The light that scatters due to the interaction of this "exciting" light with the tissue is collected and analyzed by wavelength. The exciting light is chosen to have a wavelength that is absorbed by hemoglobin. The objective is to induce as little fluorescence as possible while injecting sufficient energy into the hemoglobin. This will heat the glucose in the immediate vicinity of the excited hemoglobin and cause the intensity of anti-Stokes features to be observed. Low frequency modes, such as about 200 to about 600 cm^-1, are suitable for this purpose.

The Stokes part of the Raman spectrum obtained during step 1 is referred to as R₀(λ)_s1, and the anti-Stokes part of the Raman spectrum obtained during step 1 as R₀(λ)_as1.

<u style="single">Step 2:</u> During this step, external pressure is applied to the tissue being examined. A Raman spectrum, Stokes and anti-Stokes, is obtained by optically contacting the skin of a finger tip and the tissue/capillary bed below the skin, or another analogous section of a human body, with narrow spectral bandwidth light. The light that scatters due to the interaction of this "exciting" light with the tissue is collected and analyzed by wavelength. The exciting light has a wavelength that is absorbed by hemoglobin.

The Stokes part of the Raman spectrum obtained during step 2 is referred to as R_p(λ)_s2, and the anti-Stokes part of the Raman spectrum obtained during step 1 as R_p(λ)_as2.

An objective is to effect a temperature change discernable in the Raman spectra, thereby improving detectability of the anti-Stokes components and obtaining a different form of signal modulation of Stokes components all leading to better qualitative and quantitative analysis.

One strategy is to use the Stokes features for analysis. If fluorescence or other background radiation precludes use of Stokes features, however, the following calculation can be used with the anti-Stokes components: <maths id="math0017" num=""><math display="block"><mrow><mo mathvariant="normal">(</mo><msub><mi mathvariant="normal">R</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">1</mn></mrow></msub><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi></msub><mo>⁢</mo><msub><mfenced><mi mathvariant="normal">λ</mi></mfenced><mrow><mi>as</mi><mo>⁢</mo><mn mathvariant="normal">2</mn></mrow></msub><mo mathvariant="normal">)</mo><mo mathvariant="normal">=</mo><mi mathvariant="normal">P</mi><mfenced><mi mathvariant="normal">λ</mi></mfenced><mo>⁢</mo><msub><mi mathvariant="normal">C</mi><mi mathvariant="normal">G</mi></msub><mo>⁢</mo><mi mathvariant="normal">σ</mi><mo mathvariant="normal">;</mo></mrow></math><img file="EP1037553B1_D0017.tif" /></maths> where P(λ) is the incident power of the exciting light, C_G is glucose concentration, and σ is the cross-section function. If the values on the left side of the above equation are measured once to obtain a value which can be compared to an independent measurement of C_G, then P(λ) σ can be determined, and the instrument calibrated for use in monitoring subsequent changes in C_G.

Subtracting the results of step 1 from step 2, Stokes or anti-Stokes, gives a number that is proportional to the concentration of glucose.

An apparatus suitable for use in this embodiment of the method is shown in Figure 2. Light emitted by the first light source 100 passes through an amplitude modulator 230 to first and second reflectors 240, 250 such that the light joins the path of light emitted by the second light source 110. The light then irradiates the sample 200. In this particular embodiment, the first excitation wavelength is 785 nm and has a spectral bandwidth of 0.3 cm^-1. The background amplified spontaneous emission is preferably less than 1 part in 10⁸ of the 785 nm portion. If the 785 nm power is at least 100 mW at the skin surface, the Stokes glucose Raman spectral feature can be obtained in about 10 minutes.

Light scattered by the sample reaches a wavelength selection device 210 that is a wavelength selection spectrograph having a filter tuned to transmit a selected Raman feature such as 1008 nm to 1029 nm for skin keratin excited with 785 nm. The wavelength selection device 210 has a 250 µm slit and uses a holographic transmission grating. The scattered light is pre-filtered before entering the spectrograph slit using a holographic notch filter having a band blocking width of about 250 cm^-1 (about 1 nm) centered at 785 nm. The optical density of the filter at 785 nm is preferably at least 6 to 8. The light collection system contains no fiber optics and has an effective f number of 1.4. This can be accomplished using known optical components such as conventional 35 mm camera lenses. The focussing lenses for Raman excitation and hemoglobin absorption can be from about 2 cm to about 18 cm, depending on the size of the region to be probed.

The selected signal is then received by a detector 120. The detector 120 is a charge coupled device (CCD) array, which can be cooled with liquid nitrogen. The detector 120 can comprise a single channel detector and optionally includes spatial characteristics, such as a quadrant spatial layout, allowing simultaneous use of spatial (Hadmaard/other/tissue modulation) encoding/modulation.

The signal then passes to a phase sensitive amplifier 220 that comprises a lock-in amplifier/gated integrator set to demodulate the signal stimulated by modulator 230. The processor 130 then takes analog demodulated signal from the phase sensitive amplifier 220 and performs digitization, storage and data processing. In addition, the processor 130 can provide synchronization with concurrent tissue modulation or spatial encoding. The processor is preferably a computer.

Those skilled in the art will appreciate other variations and modifications that can be used for the method and apparatus disclosed herein.