+

WO1998034097A1 - Determination of the ratio of absorption coefficients at different wavelengths in a scattering medium - Google Patents

Determination of the ratio of absorption coefficients at different wavelengths in a scattering medium Download PDF

Info

Publication number
WO1998034097A1
WO1998034097A1 PCT/GB1998/000297 GB9800297W WO9834097A1 WO 1998034097 A1 WO1998034097 A1 WO 1998034097A1 GB 9800297 W GB9800297 W GB 9800297W WO 9834097 A1 WO9834097 A1 WO 9834097A1
Authority
WO
WIPO (PCT)
Prior art keywords
medium
light
ratio
wavelengths
characteristic
Prior art date
Application number
PCT/GB1998/000297
Other languages
French (fr)
Inventor
Matthias Kohl
Russell Watson
Mark Cope
Original Assignee
University College London
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University College London filed Critical University College London
Priority to EP98902110A priority Critical patent/EP0956495A1/en
Priority to JP53264098A priority patent/JP2001509893A/en
Priority to AU58738/98A priority patent/AU5873898A/en
Publication of WO1998034097A1 publication Critical patent/WO1998034097A1/en

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/314—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths

Definitions

  • This invention relates to the determination of the ratio of absorption coefficients at different wavelengths in a scattering medium, and to the determination of the relative concentrations of chromophores in the medium.
  • the invention provides a method of determming the ratio of the absorption coefficients ⁇ a ( ⁇ j) and ⁇ ( ⁇ _) of a light scattering medium at two different wavelengths comprising a) passing light (as herein defined) through the medium b) measuring a characteristic X of the light at two wavelengths ⁇ j and ⁇ which is affected by the absorption coefficient during its passage through the medium c) effecting changes ⁇ J ⁇ and ⁇ j ) in the absorption coefficient at the two wavelengths; d) repeating step (b) e) determining the ratio of the absorption coefficients ⁇ ,) and ⁇ j at said two wavelengths from said changes therein or from quantities representative thereof and from the characteristic X and from the transport scattering coefficients of the medium ⁇ ,) and ⁇ . ⁇ _) at the two wavelengths.
  • light as used herein includes the near infrared spectrum.
  • the said changes in the absorption coefficient may be effected by a common cause, for example by varying the concentration of a chromophore in the scattering medium, the quantities representative of the changes in the absorption coefficient being the absorption coefficients ⁇ .( ⁇ j),
  • these changes can be deliberately generated by injecting an exogenous dye (e.g. Indocyanine green) intravenously or modifying haemoglobin oxygenation or volume.
  • an exogenous dye e.g. Indocyanine green
  • the method may comprise first adding to the medium the chromophore where concentration is varied.
  • the ratio of the absorption coefficients ⁇ ,( ⁇ ) and ⁇ j ) may be determined from the ratio of the chromophore absorption coefficients ⁇ j ), ⁇ j ) and from the ratio of the transport scattering coefficients of the scattering medium ⁇ s '( _) and ⁇ s 'Q _).
  • the characteristic measured in step (b) may be the attenuation A of the light during its passage through the medium.
  • the light may be intensity modulated or pulsed light, the measured characteristic X being the phase or depth of the modulation or the time of flight (transmission time) of the light through the medium respectively.
  • the invention also provides a method of determining the variation with wavelength of the absorption coefficient of a scattering medium (the spectral characteristic) comprising performing the method as set forth above at a number of wavelengths across a waveband of interest, and in step (e) thereof determining the ratio of the absorption coefficients relative to that at a reference wavelength ⁇ ., the variation of said ratios with wavelength being the said spectral characteristic.
  • the invention provides a method of determining the ratio of chromophore concentrations of known spectral characteristics (as herein defined) in a light scattering medium comprising performing the method of steps (a) to (d) as set forth above at not less than n different wavelengths where n is the number of chromophores, one of said wavelengths being the reference wavelength ⁇ , performing step (e) to determine the ratio - ⁇ ( ⁇ )/ ⁇ .( ⁇ -) of the medium for each of the other wavelengths, and determining the relative concentrations of the chromophores in the medium from said ratios.
  • spectral characteristic we mean the variation of ⁇ t with ⁇ .
  • the method will yield the absolute spectral characteristics and hence absolute concentrations of the other chromophores.
  • the invention also includes apparatus for use in the method as set forth above.
  • Figure 2a shows the optical properties of a liquid tissue-simulating phantom used in an experiment described hereafter.
  • the absorption coefficient is that of water ( ⁇ w ) and the added dye ( ⁇ ).
  • the dotted line shows the true values derived from the water and dye absorption of the phantom ( ⁇ w + ⁇ d in figure 2a).
  • the measured attenuation changes that were induced by a small increase in ⁇ were used to estimate ⁇ ⁇ y ⁇ OOnm) by assuming a constant ⁇ ,' (dashed line, calculated from Eq. 6 below) and by including the wavelength dependence of ⁇ , 1 (solid line, Eq. 7).
  • Figure 3 shows values of the transport scattering coefficient ⁇ ,' measured on the heads of seven volunteers. The values are an estimate from measured TPSFs (data from Matcher et al 1996).
  • Figure 4 shows values of ⁇ s ' measured on the calf muscle of eleven volunteers; ⁇ ,' measured on the head of 7 volunteers. The values are an estimate from measured TPSFs (data from Matcher et al 1996).
  • Figure 5 shows apparatus according to the invention.
  • Changes in the chromophore concentration ie. changes in the absorption coefficient ( ⁇ J in a light scattering medium alter the intensity of the diffusely transmitted light collected at the surface of the medium.
  • ⁇ J in a light scattering medium alters the intensity of the diffusely transmitted light collected at the surface of the medium.
  • the spectral characteristic of the tissue absorption coefficient can be calculated from the known spectral characteristic of an injected dye absorption coefficient and the change in the measured spectrum of the intensity (attenuation), phase shift or mean time.
  • the relative concentrations of its chromophores notably the concentrations of oxy- and deoxy - haemoglobin (HbO 2 and Hb) can be calculated.
  • An application is the monitoring of the oxygen saturation (HbO 2 /(HbO 2 + Hb)) of blood in tissue.
  • the derivative of the attenuation A (defined as the logarithm of the ratio of incident and detected intensity) with respect to change in ⁇ is:
  • measurements of dA/dd ⁇ can be used to estimate the product of absorption and scattering coefficient.
  • the spectral shape of ⁇ can then be calculated from measurements of dAf ⁇ t at different wavelengths.
  • Eq. 4 states that the ratio of the absorption coefficients at different wavelengths ( ⁇ h_) can be estimated from measurements of attenuation changes ( ⁇ A) at these wavelengths, the ratio of the changes in absorption coefficient and the wavelength dependence of -'.
  • the spectra of both ⁇ .' and ⁇ ⁇ are shown in Fig. 2a.
  • a CCD spectrometer in connection with a halogen light source was employed to measure the reflected light intensity of the phantom (volume 300 ml).
  • Attenuation spectra can be corrected for the wavelength dependence of ⁇ s ':
  • Figure 3 and 4 show ⁇ . of the head and the calf muscle from volunteers obtained from measured temporal point spread functions (TPSFs).
  • TPSFs temporal point spread functions
  • the absolute ⁇ varies by up to 30% for different volunteers, however it has a similar wavelength_ characteristic. Therefore, a correction of the attenuation spectra for ⁇ s ' variations is feasible.
  • a pure change in absorption coefficient occurs and the wavelength dependence of the absorption coefficient change is known.
  • Absorption changes can be induced by an exogenous dye or by changes of haemoglobin concentration. Change in Hb or HbO 2 can be generated by small changes in the inspired gas. Alternatively, the absorption coefficient can be reduced by dilution (for example by injecting saline).
  • the dominant chromophores are deoxy- and oxyhemoglobin and water, the absorption spectra of which are known.
  • the absorption coefficient at a wavelength ⁇ can be written as
  • ⁇ n (unit: mmXmM "1 ) is the extinction coefficient of the n ⁇ chromophore and c n is its concentration (unit: mM).
  • ⁇ _ is a reference wavelength
  • an iterative technique can be used to calculate the relative chromophore concentrations c, c R (where C R signifies the concentration of a reference chromophore).
  • the number of wavelengths has to be equal to or higher than the number of chromophores in the medium.
  • the absorption spectrum of the dye used to induce attenuation changes has to be known.
  • the dye could for example be Indocyanine green.
  • the mean transit time is given by
  • the modulation depth can be written as
  • Figure 5 shows diagrammatically apparatus according to one embodiment of the invention in which changes in attenuation of light passing through the medium are measured.
  • a medium 10 under examination receives light from a white light source 12 via an optical fibre 14.
  • Diffused light is received by an optical fibre, and is taken to a wavelength dispersive device 18 and thence to a detector 20.
  • the output of the detector 20 is fed to a PC or other data processor 22.
  • the wavelength dispersive device for example a spectrometer comprising an optical slit and a diffraction grating divides up the white light issuing from the medium into its constituent wavelengths, those of interest then being selected by output slits or by being detected by respective elements if the detector 20 is of the staring-array type.
  • the detector 20 may be a charge-coupled device or an array of photodiodes.
  • the intensity of the light source 12 is kept constant, it is unnecessary to provide the detector 20 with a measure of the input light intensity. Ail necessary data is contained in changes in intensity arising from imposed changes in ⁇ and ⁇ .
  • a wavelength selective device such as a rotatably mounted filter wheel with interference filters with two or more coloured filter elements for passing the wavelengths of interest.
  • the wavelengths of interest are typically in the range 650nm to lOOOnrn for tests on tissue.
  • wavelengths of up to about 2000nm (2 ⁇ m) are appropriate.
  • the wavelength dispersive or selective device should have a bandwidth of less than 10 nm.
  • the white light source 12 is replaced by switchable or otherwise selectable monochromatic light sources (eg laser diodes or LEDs) of different wavelengths, the wavelength dispersive or selective device 18 then not being required.
  • Laser diodes are preferred to LEDS because the wavelength of the emitted light is more monochromatic.
  • a pulsed laser diode is used as the source 12 to provide pulsed light of a specific wavelength to the medium 10.
  • the wavelength dispersive device 18 is omitted and the detector 20 is replaced by means to determine the time of arrival of the pulse.
  • a signal is provided via line 24 from the drive circuitry of the laser diode to indicate the time of initiation of the pulse to the detector so that ⁇ t> may be established.
  • light source 12 is an intensity modulated monochromatic source with a choice of frequencies, either switchable laser diodes or LEDs, or a white light source with several alternative filters.
  • the detector 20 is phase-sensitive, and receives a reference signal from the input to the sample via a line 26. Phase detection may conveniently be by a lock-in technique in which the output from the medium is compared with the input from line 26 by accurately multiplying the detector signal by both the reference signal and the reference signal shifted in phase by 90°.

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

A method of determining the ratio of the absolute absorption coefficients νa(μ2) and νa(μ1) of a light scattering medium at two different wavelengths comprising: a) passing light (as herein defined) through the medium; b) measuring a characteristic X of the light at two wavelengths μ1 and μ2 which is affected by the absorption coefficient during its passage through the medium; c) effecting changes Δνa(μ1) and Δνa(μ2) in the absorption coefficients at the two wavelengths; d) repeating step (b); e) determining the ratio of the absolute absorption coefficients νa(μ1) and νa(μ2) at said two wavelengths from the said changes therein Δνa(μ1) and Δνa(μ2) or from quantities representative thereof, from the characteristic X, and from the ratio of transport scattering coefficients of the medium νs'(μ1) and νs'(μ2) at the two wavelengths.

Description

DETERMINATION OF THE RAΗO OF ABSORPΗON COEFFICIENTS AT DIFFERENT WAVELENGTHS IN A SCATTERING MEDIUM
This invention relates to the determination of the ratio of absorption coefficients at different wavelengths in a scattering medium, and to the determination of the relative concentrations of chromophores in the medium.
There are a number of industrial processes in which measurement of the ratio of chromophore concentrations in a scattering medium would (if available) be an important indicator for process quality control. Examples are the measurement of water, protein and fat content in liquids and solids in bioreactors and food processing, and the measurement of relative dye concentrations in dyeing processes, which relative concentrations are critical to the consistency of the final colour achieved in the dyed product.
Other industrial applications would be the measurement of the ratio of chromophore concentrations in tissue samples, for example relative oxy-and deoxy haemoglobin concentrations. Measurement of such quantities in vivo would also be useful in support of clinical practice, as would measurement of the oxygen saturation of blood in tissue.
It is an object of at least the preferred embodiment of the invention to provide a method of determining relative concentrations of chromophores in a scattering medium, and also to provide a method of determining relative absorption coefficients of the medium itself, since this may be sufficient for some applications.
In a first aspect the invention provides a method of determming the ratio of the absorption coefficients μa(λj) and μ (λ_) of a light scattering medium at two different wavelengths comprising a) passing light (as herein defined) through the medium b) measuring a characteristic X of the light at two wavelengths λj and λ^ which is affected by the absorption coefficient during its passage through the medium c) effecting changes ΔμJλ and Δμ^λj) in the absorption coefficient at the two wavelengths; d) repeating step (b) e) determining the ratio of the absorption coefficients μ^λ,) and μ^λj at said two wavelengths from said changes therein or from quantities representative thereof and from the characteristic X and from the transport scattering coefficients of the medium μ ,) and μ.\λ_) at the two wavelengths.
The term "light" as used herein includes the near infrared spectrum.
The said changes in the absorption coefficient may be effected by a common cause, for example by varying the concentration of a chromophore in the scattering medium, the quantities representative of the changes in the absorption coefficient being the absorption coefficients Δμ.(λj),
Figure imgf000004_0001
In tissue, these changes can be deliberately generated by injecting an exogenous dye (e.g. Indocyanine green) intravenously or modifying haemoglobin oxygenation or volume.
The method may comprise first adding to the medium the chromophore where concentration is varied.
The ratio of the absorption coefficients μ,(λι) and μ^λj) may be determined from the ratio of the chromophore absorption coefficients Δμ^λj), Δμ^λj) and from the ratio of the transport scattering coefficients of the scattering medium μs'( _) and μs'Q _).
In one embodiment the characteristic measured in step (b) may be the attenuation A of the light during its passage through the medium.
In another embodiment the light may be intensity modulated or pulsed light, the measured characteristic X being the phase or depth of the modulation or the time of flight (transmission time) of the light through the medium respectively.
The invention also provides a method of determining the variation with wavelength of the absorption coefficient of a scattering medium (the spectral characteristic) comprising performing the method as set forth above at a number of wavelengths across a waveband of interest, and in step (e) thereof determining the ratio of the absorption coefficients relative to that at a reference wavelength λ., the variation of said ratios with wavelength being the said spectral characteristic.
In a further aspect the invention provides a method of determining the ratio of chromophore concentrations of known spectral characteristics (as herein defined) in a light scattering medium comprising performing the method of steps (a) to (d) as set forth above at not less than n different wavelengths where n is the number of chromophores, one of said wavelengths being the reference wavelength λ,, performing step (e) to determine the ratio -ι(λι)/μ.(λ-) of the medium for each of the other wavelengths, and determining the relative concentrations of the chromophores in the medium from said ratios. By spectral characteristic we mean the variation of μt with λ.
Where the absolute value of the absorption coefficient of one of the chromophores is known (eg water) or can be found by other means, the method will yield the absolute spectral characteristics and hence absolute concentrations of the other chromophores.
The invention also includes apparatus for use in the method as set forth above.
TTie invention will now be described merely by way of example with reference to the accompanying figures, wherein
Figure 1 shows the fractional error in the ratio (^( jJ/μJ j) from the approximation in equations 3 and 4, discussed below. Specifically, it shows fractional error in ratio μ»(λ2) μ*(λι) from approximations of equations 3 and 4 assuming μ£ ) = 0.02 mm-1, μ,1^ = μ,'(λι) = 1 mm"1 and a source detector r = 50 mm.
Figure 2a shows the optical properties of a liquid tissue-simulating phantom used in an experiment described hereafter. The absorption coefficient is that of water (μ^w) and the added dye (μ^). The transport scattering coefficient (μ,1) is based upon calculations from Mie theory for the concentration c, = 2% and the diameter (1 μm) of the spheres.
Figure 2b shows the spectrum of the absorption coefficient of the phantom normalised at λ=800 nm. Specifically, it shows spectrum of the relative absorption coefficient (μ^λ μ^δOOnm)) of the phantom normalised at λ = 800nm. The dotted line shows the true values derived from the water and dye absorption of the phantom (μ^w + μ^d in figure 2a). The measured attenuation changes that were induced by a small increase in μ^ were used to estimate μ λyμ^δOOnm) by assuming a constant μ,' (dashed line, calculated from Eq. 6 below) and by including the wavelength dependence of μ,1 (solid line, Eq. 7).
Figure 3 shows values of the transport scattering coefficient μ,' measured on the heads of seven volunteers. The values are an estimate from measured TPSFs (data from Matcher et al 1996).
Figure 4 shows values of μs' measured on the calf muscle of eleven volunteers; μ,' measured on the head of 7 volunteers. The values are an estimate from measured TPSFs (data from Matcher et al 1996). Figure 5 shows apparatus according to the invention.
Changes in the chromophore concentration ie. changes in the absorption coefficient (ΔμJ in a light scattering medium alter the intensity of the diffusely transmitted light collected at the surface of the medium. Here we describe how the spectral characteristic of the tissue absorption coefficient can be calculated from the known spectral characteristic of an injected dye absorption coefficient and the change in the measured spectrum of the intensity (attenuation), phase shift or mean time.
Once the spectral characteristic of the tissue absorption coefficient has been determined, the relative concentrations of its chromophores, notably the concentrations of oxy- and deoxy - haemoglobin (HbO2 and Hb), can be calculated. An application is the monitoring of the oxygen saturation (HbO2/(HbO2 + Hb)) of blood in tissue.
Diffusion Theory Model
For a semi-infinite halfspace the reflectance at a distance r from the light source can be written as
_ l _.n ^ exp(-μ5ff -p)
R(r) = zQ -(-+μt„) p 2rc-p (1)
In Eq. 1 p = (r2+z0 2)'Λ and z0 = l/μs', where μs' is the transport scattering coefficient. The velocity of light in the medium c = ς/n (where the speed of light in vacuum is c0).
μeII = (3.μ. fa+μ.')) is the effective attenuation coefficient and D = (3.(μ^+μ ))'1 is the diffusion coefficient. The derivative of the attenuation A (defined as the logarithm of the ratio of incident and detected intensity) with respect to change in μ^ is:
dA
■ • (2μ, +μ,' ) (2)
3μ, 2- lnlO 1/p + μ .a
Here it is suggested that measurements of dA/ddμ, can be used to estimate the product of absorption and scattering coefficient. For a known wavelength dependence of μ. the spectral shape of μ can then be calculated from measurements of dAfμt at different wavelengths.
Assuming that scattering dominates absorption (μs'» μ_ Eq. 2 can be approximated to:
Figure imgf000007_0001
Furthermore, for large source-detector distances it is:
(4)
Figure imgf000007_0002
Eq. 4 states that the ratio of the absorption coefficients at different wavelengths (λ h_) can be estimated from measurements of attenuation changes (ΔA) at these wavelengths, the ratio of the changes in absorption coefficient and the wavelength dependence of -'.
μ. (λι ) ΔA(λ ) Δμa(λt ) μ-' (λj ) μa(λι) A(λt) Δμa(λ2)J μs ' (λ2) (5)
To estimate the validity of the approximations leading from Eq. 2 to Eq.s 3 and 4, the errors in the calculated ratio
Figure imgf000007_0003
assuming μ (λ.) = 0.002 mm-1, ^'(λ^) = μ.'(λ_) = 1 mm"1 and a source-detector distance r = 50 mm. In the figure the fractional errors of Eq.s 3 and 4 are shown in solid and dotted line, respectively. The fractional errors increase with increasing differences of
Figure imgf000007_0004
Variations in the magnitude of μ of a factor of 2 result in errors of approximately 10%. Determination of the relative absorption spectrum of a light-scattering phantom
The method suggested above was tested in a liquid, light scattering tissue simulating phantom of known optical properties. Silica spheres (Monospher M1000, Merck, Germany) of 1 μm diameter were suspended in water and served as scattering centres in the phantom. The scattering properties could be derived from Mie theory calculations as the diameter of the spheres and the refractive index of both the spheres and the surrounding medium are known. The calculated scattering cross section of the spheres was compared with measurements from a collimated beam setup. Experimental and theoretical scattering cross sections were in agreement when the wavelength dependence of the refractive index of fused silica (n. = 1.447 + 3890λ"2 - 37 "4 with λ in nm) and water (μ_ = 1.32 + 6878λ"2 - 1.132.10"9 λ4) were employed. The transport scattering coefficient μ.' was calculated for wavelengths between 600 and 1050 nm for these refractive indices and the concentration of the spheres in the phantom (cs = 2% volJvol.). The calculated μ.' decreases approximately linearly with wavelength and a first order regression through the calculated values gave ^,'=0.92 - 0.00047 λfmm"1] (λ in nm). The absorption coefficient of the phantom consisted of two components: the water absorption u^ and the absorption of a dye (S 109564, Id, GB) (μ^_) added at a concentration of cd = 0.0031% volJvol.. The spectra of both μ.' and μλ are shown in Fig. 2a.
A CCD spectrometer in connection with a halogen light source was employed to measure the reflected light intensity of the phantom (volume 300 ml). The ends of the optical fibres that transported the light from the lamp to the phantom and from the phantom to the detector were submerged approximately 2 mm in the phantom at a source-detector distance of r = 38 mm. Two spectra (I0 and I) of the reflected light intensity were recorded before and after the concentration of the dye was changed by Δcd = 2.5 x 10^% . This corresponds to a change in μ of between 4% and 8% depending of wavelength. The spectrum of the attenuation change ΔA(λ)=log[I0(λ)/I(λ)] was calculated. Ignoring the wavelength dependence of scattering gives a first approximation of the spectral characteristic μ λ) according to Eq. 4/5, normalised a) for the relative magnitude of the absorption change (ie. divided by the absorption spectrum of the dye μ^_) and b) for the value at λ = 800 nm:
μa(λ) _ [" ΔA(800πm) μa,d(λ) μa(800nm) ' |_ΔA(λ) μa>d (800nm)j ' ^ This first approximation of u1(λ)(dashed line in Fig. 2b) shows considerable discrepancies with the true /^(λy ^δOOnm) spectrum (dotted line). Correcting the experimental data for the wavelength dependent scattering coefficient of the phantom (μs'(K) = 0.92 mm"1 - 0.00047λ[mm"1].
μ»fl ΔA(800n ) μa.d(λ) μs' (λ) μa (800am) ΔA(λ) μatd(800nm) μs' (800nm)
(7)
gives the spectrum shown plotted with a solid line in Fig. 2b. The experimental data are in good agreement with the expected μ spectrum for λ < 950 mm. For longer wavelengths the experimental data are noisy and underestimate μ . The poor signal to noise ratio for these wavelengths is due to the low sensitivity of the CCD detector and the higher absorption coefficient of the phantom. The likely explanation for the underreading at these wavelengths is, as demonstrated in Fig.l, that the approximations used for the derivation of Eq. 3-5 break down for larger ratios .i(λ)/ .1(800nm).
Calculation of relative chromophore concentrations from the relative absorption spectrum
In the experiment, the determination of ratio of chromophore concentrations and the oxygenisaturation of blood depends on the following points: a) Attenuation spectra can be corrected for the wavelength dependence of μs': Figure 3 and 4 show μ. of the head and the calf muscle from volunteers obtained from measured temporal point spread functions (TPSFs). The absolute μ varies by up to 30% for different volunteers, however it has a similar wavelength_ characteristic. Therefore, a correction of the attenuation spectra for μs' variations is feasible.
b) A pure change in absorption coefficient occurs and the wavelength dependence of the absorption coefficient change is known. Absorption changes can be induced by an exogenous dye or by changes of haemoglobin concentration. Change in Hb or HbO2 can be generated by small changes in the inspired gas. Alternatively, the absorption coefficient can be reduced by dilution (for example by injecting saline).
c) The spectral characteristic of all tissue chromophores is known:
For tissue and wavelengths in the near infrared the dominant chromophores are deoxy- and oxyhemoglobin and water, the absorption spectra of which are known. In general the absorption coefficient at a wavelength λ can be written as
where εn (unit: mmXmM"1) is the extinction coefficient of the nΛ chromophore and cn is its concentration (unit: mM). Once the ratio μj^lμ ^) at least at, n wavelengths has been determined (λ_ is a reference wavelength), an iterative technique can be used to calculate the relative chromophore concentrations c, cR (where CR signifies the concentration of a reference chromophore).
The number of wavelengths has to be equal to or higher than the number of chromophores in the medium. The absorption spectrum of the dye used to induce attenuation changes has to be known. For the determination of the relative oxy- /deoxy haemoglobin concentrations in tissue, the dye could for example be Indocyanine green.
In the discussion above measured changes in the light intensity were used. Equally changes in the time of flight, the phase shift or the modulation depth of an intensity modulated light wave can be exploited.
Thus changes in the mean (transit) time (<t>) of the photons or the phase (Φ) of an intensity light_wave (modulated at the frequency v^) can be used to compute the spectral characteristic of μ_-
The mean transit time is given by
P2 (8)
2c-(D+pχμa - D) and the phase is approximately a linear function of <t>: φ = -2π - vM - < t >
(9)
The derivation of <t> with respect to μ, can be approximated as
d < t > -p f- = - 3 - dμa 4- c (10) μa μa
Therefore,
Figure imgf000011_0001
The equivalent to Eq.5 with respect to measured changes in mean transit time is therefore
Figure imgf000011_0002
The measurements of phase changes ΔΦ, Eq. 12 can be combined with Eq. 9 to give
Figure imgf000011_0003
The modulation depth can be written as
Figure imgf000011_0004
where ψ0 = μe£f.p[l + χ2]I 4 , ψr = -ψ0 .sin (θ/2), ψi = -ψ0 .cos (θ/2), θ = tan_1(χ), ψ„ = μe£f.p and χ = (2πvNi)/(μ1c).
Changes in μa induce changes ΔM. It can be shown that equivalently to Eq.7, 12 and 13 the ratio of absorption coefficients at two different wavelengths can be derived from the measured fractional change δM = ΔM M at these wavelengths according to
Figure imgf000012_0001
Figure 5 shows diagrammatically apparatus according to one embodiment of the invention in which changes in attenuation of light passing through the medium are measured. A medium 10 under examination receives light from a white light source 12 via an optical fibre 14.
Diffused light is received by an optical fibre, and is taken to a wavelength dispersive device 18 and thence to a detector 20. The output of the detector 20 is fed to a PC or other data processor 22.
The wavelength dispersive device, for example a spectrometer comprising an optical slit and a diffraction grating divides up the white light issuing from the medium into its constituent wavelengths, those of interest then being selected by output slits or by being detected by respective elements if the detector 20 is of the staring-array type. The detector 20 may be a charge-coupled device or an array of photodiodes.
Provided the intensity of the light source 12 is kept constant, it is unnecessary to provide the detector 20 with a measure of the input light intensity. Ail necessary data is contained in changes in intensity arising from imposed changes in λ and μ .
As an alternative to the spectrometer, a wavelength selective device may be employed such as a rotatably mounted filter wheel with interference filters with two or more coloured filter elements for passing the wavelengths of interest.
The wavelengths of interest are typically in the range 650nm to lOOOnrn for tests on tissue. When the distance between light source and detector fibers is small (<5mm) wavelengths of up to about 2000nm (2μm) are appropriate. Preferably the wavelength dispersive or selective device should have a bandwidth of less than 10 nm.
In an alternative embodiment the white light source 12 is replaced by switchable or otherwise selectable monochromatic light sources (eg laser diodes or LEDs) of different wavelengths, the wavelength dispersive or selective device 18 then not being required. Laser diodes are preferred to LEDS because the wavelength of the emitted light is more monochromatic.
When the invention is to be implemented by measuring the time of flight <t>, a pulsed laser diode is used as the source 12 to provide pulsed light of a specific wavelength to the medium 10.
The wavelength dispersive device 18 is omitted and the detector 20 is replaced by means to determine the time of arrival of the pulse. A signal is provided via line 24 from the drive circuitry of the laser diode to indicate the time of initiation of the pulse to the detector so that <t> may be established.
For implementing the invention by measuring the phase shift Δφ or depth of modulation of light, light source 12 is an intensity modulated monochromatic source with a choice of frequencies, either switchable laser diodes or LEDs, or a white light source with several alternative filters. The detector 20 is phase-sensitive, and receives a reference signal from the input to the sample via a line 26. Phase detection may conveniently be by a lock-in technique in which the output from the medium is compared with the input from line 26 by accurately multiplying the detector signal by both the reference signal and the reference signal shifted in phase by 90°.
Each feature disclosed in this specification (which term includes the claims) and/or shown in the drawings may be incorporated in the invention independently of other disclosed and/or illustrated features.

Claims

1. A method of determining the ratio of the absolute absorption coefficients ╬╝^╬╗^ and ╬╝^╬╗JJ of a light scattering medium at two different wavelengths comprising
a) passing light (as herein defined) through the medium b) measuring a characteristic X of the light at two wavelengths λj and ^ which is affected by the absorption coefficient during its passage through the medium c) effecting changes Δμl(λ:) and
Figure imgf000014_0001
in the absorption coefficients at the two wavelengths; d) repeating step (b) e) determining the ratio of the absolute absorption coefficients μ.(λx) and μjfa) at said two wavelengths from the said changes therein ΔμJλ and Δμ^λj) or from quantities representative thereof, from the characteristic X, and from the ratio of transport scattering coefficients of the medium μ3'(λ:) and μ3'(7 _) at the two wavelengths.
2. A method as claimed in claim 1 comprising effecting said changes in the absorption coefficients by means of a common cause.
3. A method as claimed in claim 1 or claim 2 comprising effecting said changes in the absorption coefficient
Figure imgf000014_0002
and Δμ [ _) by varying the concentration of an existing chromophore in the scattering medium.
4. A method as claimed in claim 1 or claim 2 comprising effecting said changes in the absorption coefficients Δμ^λJJ and μ (λ_) by adding a chromophore" to the scattering medium.
5. A method as claimed in claim 3 or claim 4 wherein the ratio of the effected absorption coefficient changes
Figure imgf000014_0003
and A╬╝.Q _) is determined from the ratio of the absolute chromophore absorption coefficients ╬╝^╬╗j), ╬╝ldQ _).
6. A method as claimed in any preceding claim wherein the characteristic X measured in step (b) is the attenuation A of the light during its passage through the medium.
7. A method as claimed in any of claims 1 to 5 wherein the light is intensity modulated light, the measured characteristic X being the phase, ac intensity or depth of the modulation of the light through the medium.
8. A method as claimed in claims 1 to 5 wherein the light is a short pulse of light, the measure characteristic X being the time of flight (as mean time of transmission) of the light through the medium.
9. A method of determining the variation with wavelength of the ratio of absorption coefficients of a scattering medium (the relative spectral characteristic) comprising performing the method of any of claims 1 to 8 at a number of wavelengths across a waveband of interest, and in step (e) thereof determining the ratio of the absorption coefficients relative to that at a reference wavelength ╬╗_, the variation of said ratios with wavelength being the said spectral characteristic.
10. A method of determining the ratio of concentrations of chromophores of known spectral characteristic (as herein defined) in a light scattering medium comprising performing the method of steps (a) to (d) of any preceding claim at not less than n different wavelengths where n is the number of chromophores, one of said n wavelengths being the reference wavelength ╬╗-, performing step (e) to determine the ratio ╬╝^╬╗ /╬╝/╬╗-) of the medium for each of the other wavelengths, and determining the ratio of concentrations of the chromophores in the medium from said ratios.
11. A method of determining the absolute concentration of chromophores of known spectral characteristic, as claimed in claim 9, where one of the chromophores is of known concentration.
12. Apparatus for determining the ratio of the absolute absorption coefficient μx(λ_) and μ,(λj) of a light scattering medium at two different wavelengths comprising a) means for passing light (as herein defined) through the medium; b) means for measuring a characteristic X of the light at least two wavelengths λ. and λ^ which is affected by the absorption coefficient during its passage through the medium; c) means for effecting changes Δμx(λj) and Δμ^λ^) in the absorption coefficient at the two wavelengths; d) means for determining the ratio of the absolute absorption coefficients μ_(λj) and μ.Q _) at said two wavelengths from the said changes therein Δμ^λ and Δμ Q _) or from quantities representative thereof, from the characteristic X, and from the ratio of transport scattering coefficients of the medium μ- λj) and μ,'Q _) at the two wavelengths.
13. Apparatus as claimed in claim 12 wherein the means (b) is adapted to measure the attenuation A of the light during its passage through the medium as said characteristic X.
14. Apparatus as claimed in claim 12 wherein the light is intensity modulated light and means (b) is adapted to measure the phase or the intensity or the depth of modulation of the modulated light through the medium as said characteristic X.
15. Apparatus as claimed in claim 12 wherein the light is a short pulse of light, and means (b) is adapted to measure the time of flight (mean time of transmission) of the light through the medium as said characteristic X.
16. Apparatus for determining the ratio of concentrations of chromophores of known spectral characteristic (as herein defined) in a light scattering medium comprising means performing the method of steps (a) to (d) of any of claims 1 to 11 at not less than n different wavelengths where n is the number of chromophores, one of said n wavelengths being the reference wavelength λ., means for performing step (e) to determine the ratio Δμ^λ and Aμ } _) of the medium for each of the other wavelengths, and means for determining the ratio of concentrations of the chromophores in the medium from said ratios.
17. Apparatus as claimed in any of claims 12 to 16 wherein the means (c) is arranged to effect said changes by a common cause.
PCT/GB1998/000297 1997-01-31 1998-01-30 Determination of the ratio of absorption coefficients at different wavelengths in a scattering medium WO1998034097A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP98902110A EP0956495A1 (en) 1997-01-31 1998-01-30 Determination of the ratio of absorption coefficients at different wavelengths in a scattering medium
JP53264098A JP2001509893A (en) 1997-01-31 1998-01-30 Determination of the ratio of absorption coefficients at different wavelengths in a scattering medium
AU58738/98A AU5873898A (en) 1997-01-31 1998-01-30 Determination of the ratio of absorption coefficients at different wavelengt hs in a scattering medium

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GBGB9702018.4A GB9702018D0 (en) 1997-01-31 1997-01-31 Determination of the ratio of optical absorbtion coefficients at different wavelengths in a scattering medium
GB9702018.4 1997-01-31

Publications (1)

Publication Number Publication Date
WO1998034097A1 true WO1998034097A1 (en) 1998-08-06

Family

ID=10806893

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/GB1998/000297 WO1998034097A1 (en) 1997-01-31 1998-01-30 Determination of the ratio of absorption coefficients at different wavelengths in a scattering medium

Country Status (5)

Country Link
EP (1) EP0956495A1 (en)
JP (1) JP2001509893A (en)
AU (1) AU5873898A (en)
GB (1) GB9702018D0 (en)
WO (1) WO1998034097A1 (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001046679A1 (en) * 1999-12-22 2001-06-28 Applied Optics Center Of Delaware, Inc. Method and apparatus for analyzing samples in a clinical analyzer using coherent radiation
EP1136811A1 (en) * 1998-11-05 2001-09-26 Hamamatsu Photonics K.K. Method and device for measuring internal information of scattering absorber
WO2003011127A1 (en) * 2001-08-02 2003-02-13 The Electrode Company Limited Optoelectronic blood analytical apparatus
US6591122B2 (en) * 2001-03-16 2003-07-08 Nellcor Puritan Bennett Incorporated Device and method for monitoring body fluid and electrolyte disorders
CN1297232C (en) * 2004-05-21 2007-01-31 天津大学 Optical-circuit-variable airspace light-dividing differencial wavelength spectometer for detecting tissue content and detection method thereof
CN1299646C (en) * 2004-05-21 2007-02-14 天津大学 Optical-circuit-variable time-domain light-dividing differential wave length spectrometer for detecting tissue content and detection method thereof
US7239902B2 (en) 2001-03-16 2007-07-03 Nellor Puritan Bennett Incorporated Device and method for monitoring body fluid and electrolyte disorders
US7277741B2 (en) 2004-03-09 2007-10-02 Nellcor Puritan Bennett Incorporated Pulse oximetry motion artifact rejection using near infrared absorption by water
WO2007104925A3 (en) * 2006-03-11 2007-12-06 Univ Durham Optical transmissometer and light source and light detector for such optical transmissometer
US7643858B2 (en) 2006-09-28 2010-01-05 Nellcor Puritan Bennett Llc System and method for detection of brain edema using spectrophotometry
US8690864B2 (en) 2007-03-09 2014-04-08 Covidien Lp System and method for controlling tissue treatment
CN103969223A (en) * 2014-04-18 2014-08-06 南京信息工程大学 Device for measuring atmospheric visibility in hazy weather

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103630506B (en) 2012-08-20 2016-10-26 台医光电科技股份有限公司 Detection module and detection device

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2030696A (en) * 1978-08-29 1980-04-10 Kommandiittyhtio Finnpipette O A method for semiquantitative automatic measurement of colour intensity of turbidity of a liquid solution
US4305659A (en) * 1980-03-06 1981-12-15 Baxter Travenol Laboratories, Inc. Photometric apparatus and method
US4329149A (en) * 1980-03-06 1982-05-11 Hach Chemical Company Method for spectrophotometric compensation for colorimetric reagent variation
US4448889A (en) * 1982-06-18 1984-05-15 Instrumentation Laboratory Inc. Fluid analysis
EP0210417A1 (en) * 1985-06-21 1987-02-04 Radiometer A/S A method and an apparatus for determining blood components
US5281396A (en) * 1989-12-01 1994-01-25 Baxter International Inc. Automated test system for evaluating the physical compatibility of intravenous drugs in solutions
EP0619981A1 (en) * 1993-04-12 1994-10-19 Ohmeda Inc. Arterial blood monitoring probe

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2030696A (en) * 1978-08-29 1980-04-10 Kommandiittyhtio Finnpipette O A method for semiquantitative automatic measurement of colour intensity of turbidity of a liquid solution
US4305659A (en) * 1980-03-06 1981-12-15 Baxter Travenol Laboratories, Inc. Photometric apparatus and method
US4329149A (en) * 1980-03-06 1982-05-11 Hach Chemical Company Method for spectrophotometric compensation for colorimetric reagent variation
US4448889A (en) * 1982-06-18 1984-05-15 Instrumentation Laboratory Inc. Fluid analysis
EP0210417A1 (en) * 1985-06-21 1987-02-04 Radiometer A/S A method and an apparatus for determining blood components
US5281396A (en) * 1989-12-01 1994-01-25 Baxter International Inc. Automated test system for evaluating the physical compatibility of intravenous drugs in solutions
EP0619981A1 (en) * 1993-04-12 1994-10-19 Ohmeda Inc. Arterial blood monitoring probe

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6704110B2 (en) 1998-11-05 2004-03-09 Hamamatsu Photonics K.K. Method and apparatus for measuring internal information of scattering medium
EP1136811A1 (en) * 1998-11-05 2001-09-26 Hamamatsu Photonics K.K. Method and device for measuring internal information of scattering absorber
EP1136811A4 (en) * 1998-11-05 2003-01-02 Hamamatsu Photonics Kk Method and device for measuring internal information of scattering absorber
WO2001046679A1 (en) * 1999-12-22 2001-06-28 Applied Optics Center Of Delaware, Inc. Method and apparatus for analyzing samples in a clinical analyzer using coherent radiation
US7236811B2 (en) 2001-03-16 2007-06-26 Nellcor Puritan Bennett Incorporated Device and method for monitoring body fluid and electrolyte disorders
US6591122B2 (en) * 2001-03-16 2003-07-08 Nellcor Puritan Bennett Incorporated Device and method for monitoring body fluid and electrolyte disorders
US7239902B2 (en) 2001-03-16 2007-07-03 Nellor Puritan Bennett Incorporated Device and method for monitoring body fluid and electrolyte disorders
US8457722B2 (en) 2001-03-16 2013-06-04 Covidien Lp Device and method for monitoring body fluid and electrolyte disorders
WO2003011127A1 (en) * 2001-08-02 2003-02-13 The Electrode Company Limited Optoelectronic blood analytical apparatus
US7277741B2 (en) 2004-03-09 2007-10-02 Nellcor Puritan Bennett Incorporated Pulse oximetry motion artifact rejection using near infrared absorption by water
CN1297232C (en) * 2004-05-21 2007-01-31 天津大学 Optical-circuit-variable airspace light-dividing differencial wavelength spectometer for detecting tissue content and detection method thereof
CN1299646C (en) * 2004-05-21 2007-02-14 天津大学 Optical-circuit-variable time-domain light-dividing differential wave length spectrometer for detecting tissue content and detection method thereof
WO2007104925A3 (en) * 2006-03-11 2007-12-06 Univ Durham Optical transmissometer and light source and light detector for such optical transmissometer
US7643858B2 (en) 2006-09-28 2010-01-05 Nellcor Puritan Bennett Llc System and method for detection of brain edema using spectrophotometry
US8690864B2 (en) 2007-03-09 2014-04-08 Covidien Lp System and method for controlling tissue treatment
CN103969223A (en) * 2014-04-18 2014-08-06 南京信息工程大学 Device for measuring atmospheric visibility in hazy weather

Also Published As

Publication number Publication date
JP2001509893A (en) 2001-07-24
AU5873898A (en) 1998-08-25
EP0956495A1 (en) 1999-11-17
GB9702018D0 (en) 1997-03-19

Similar Documents

Publication Publication Date Title
US6078833A (en) Self referencing photosensor
US5402778A (en) Spectrophotometric examination of tissue of small dimension
JP3526652B2 (en) Optical measuring method and optical measuring device
US5324945A (en) Method of nondestructively measuring sugar content of fruit by using near infrared transmittance spectrum
CA1273223A (en) Methods and apparatus for measuring the light absorbance of a fluid medium
US8064051B2 (en) Apparatus and method for spectrophotometric analysis
US5898487A (en) Apparatus and method for determining the concentrations of hemoglobin derivatives
WO1998034097A1 (en) Determination of the ratio of absorption coefficients at different wavelengths in a scattering medium
KR890701053A (en) Spectrophotometric method to quantitatively determine the concentration of lean components in light or radiation scattering environments
US5696580A (en) Method of and apparatus for measuring absorbance, component concentration or specific gravity of liquid sample
CA3025935A1 (en) Photothermal interferometry apparatus and method
WO1999040432A1 (en) Optical measurement of blood hematocrit incorporating a self-calibration algorithm
AU2002349800A1 (en) Spectroscopic fluid analyzer
Jensen et al. Near-infrared transmission spectroscopy of aqueous solutions: Influence of optical pathlength on signal-to-noise ratio
EP0033253B1 (en) Device for the photometrical measurement of complex solutions with varying background noise
Pérez et al. Optical Fiber Sensors for Chemical and Biological
Chen et al. Light transmittance through a region of an intact fruit
Schenkman et al. Near-infrared spectroscopic measurement of myoglobin oxygen saturation in the presence of hemoglobin using partial least-squares analysis
US4652136A (en) Method for the measurement of the difference in color contents of samples
US10241044B2 (en) NDIR glucose detection in liquids
US6995835B2 (en) Method and apparatus for measuring analytes in blood bags
BIRTH RADIOMETRIC MEASUREMENT OF FOOD QUALITY‐A REVIEW
WO2006097437A1 (en) Monitoring of predetermined substances in blood
Räty et al. Resolving optical constants from reflectance of liquids in the UV-visible range
JPH11248622A (en) Urinalysis device

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GE GH GM GW HU ID IL IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG US UZ VN YU ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW SD SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 09364204

Country of ref document: US

ENP Entry into the national phase

Ref country code: JP

Ref document number: 1998 532640

Kind code of ref document: A

Format of ref document f/p: F

WWE Wipo information: entry into national phase

Ref document number: 1998902110

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1998902110

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

WWR Wipo information: refused in national office

Ref document number: 1998902110

Country of ref document: EP

WWW Wipo information: withdrawn in national office

Ref document number: 1998902110

Country of ref document: EP

点击 这是indexloc提供的php浏览器服务，不要输入任何密码和下载