CA2377724A1 - Methods and systems using field-based light scattering spectroscopy - Google Patents

Methods and systems using field-based light scattering spectroscopy Download PDF

Info

Publication number
CA2377724A1
CA2377724A1 CA002377724A CA2377724A CA2377724A1 CA 2377724 A1 CA2377724 A1 CA 2377724A1 CA 002377724 A CA002377724 A CA 002377724A CA 2377724 A CA2377724 A CA 2377724A CA 2377724 A1 CA2377724 A1 CA 2377724A1
Authority
CA
Canada
Prior art keywords
light
region
interest
wavelength
tissue
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
CA002377724A
Other languages
French (fr)
Inventor
Changhuei Yang
Adam P. Wax
Lev T. Perelman
Ramachandra R. Dasari
Michael S. Feld
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Massachusetts Institute of Technology
Original Assignee
Individual
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 Individual filed Critical Individual
Publication of CA2377724A1 publication Critical patent/CA2377724A1/en
Abandoned legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • G01N21/4795Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0062Arrangements for scanning
    • A61B5/0066Optical coherence imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods

Abstract

The present invention relates to systems and methods of field-based light scattering spectroscopy. These systems and methods provide for the diagnosis of tissue by measuring the size and distribution of cellular characteristics . Field based measurements provide phase information resulting from the interaction of scatterers within the material and the incident wavefront. These measurements can be used to provide three dimensional images of tissue .

Description

METHODS AND SYSTEMS USING FIELD-BASED LIGHT SCATTERING
SPECTROSCOPY
RELATED APPLICATIONS
This application claims the benefit of U.S. Application No. 60/200,187 filed April 28, 2000. The entire teachings of the above application is incorporated herein by reference.
GOVERNMENT SUPPORT
The invention was supported, in whole or in part, by a grant P41-RR02594 from the National Institutes for Health. The Government has certain rights in the invention.
BACKGROUND OF THE INVENTION
Optical imaging techniques have been shown to be useful diagnostic tools for a number of biomedical applications. Increasingly, the structure of biological tissue is being studied using the interaction of the tissue with light. For example, using light scatteung spectroscopy (LSS) the size and index of refraction of scatterers comparable in size to the wavelength of light are profiled by measuring the periodicity of the spectra of reflected light. In standard pathological analysis, sensitivity to variations in the size and index of refraction of cell nuclei can give indications of abnormal cell activity. For example, pre-cancerous epithelial cells exhibit nuclear enlargement.

Generally, obtaining spectroscopic information from tissue in vivo is difficult as light scattered from surrounding tissues often obscures the desired optical signal. This can be overcome through a combination of spectroscopy with a biomedical imaging techniques which provide the ability to isolate the light scattered from a specific portion of a biological. sample.. For example, light scattering techniques have been used during endoscopic procedures to measure the size distribution of cell nuclei and refractive index changes in the epithelial linings of the body. In these studies, the intensity of white light back-scattered from the tissue is collected via an optical fiber probe and spectrally analyzed. The cell nuclei behave like Mie scatterers. Such particles exhibit periodic intensity variations with wavelength that are proportional to their sizes (typically 5-15 ~,m) and relative refractive indices. Thus light scattering spectroscopy (LSS) is of interest because changes in the size of cell nuclei and their chromatin content (related to refractive index) are primary indicators of dysplasia, the precursor of cancer, and treatment is most simple and effective when implemented at this early stage. However, intensity based LSS only provides a two-dimensional image.
Optical Coherence Tomography (OCT) is another technique that has been developed as a diagnostic tool for the study of biological tissue especially for noninvasive cross-sectional imaging in biological systems. OCT uses low-coherence interferometry to produce a three-dimensional image of optical scattering from internal tissue microstructures in a way that is analogous to ultrasonic pulse-echo imaging. OCT has longitudinal and lateral spatial resolutions of a few micrometers and can detect reflected signals as small as about 10-1° of the incident optical power.
In OCT, tissue is placed in a Michelson interferometer illuminated by a broadband light source. Due to the limited coherence length of the source (typically 10-15 ~.m), light returning from the interferometer reference arm and light backscattered by internal sample reflections interferes constructively or destructively only when the interferometer arrn optical path lengths are matched to within the source coherence length. Scanning the reference arm length generates a localized interference pattern to appear in the detector current for every internal sample reflection as a function of depth along the sample arm beam. A sample containing many reflection sites distributed along its depth (such as biological tissue) generates a detector current that contains the sum of multiple, overlapping copies of this interference pattern. A map of tissue reflectivity versus depth (called an A-Scan) is obtained by scanning the reference mirror while recording the detector current. The envelope, or outline of the detector current, may be recorded with high dynamic range by scanning the reference mirror at fixed velocity, while demodulating the detector current at the resulting Doppler frequency. Cross-sectional images of tissue bacl~scatter (called B-Scans) are acquired by collecting sequential A-scans while scanning the probe beam across the tissue surface. The resulting two-dimensional datasets are plotted as gray-scale or false-color images. However, OCT does not provide a functional image as the resultant resolution does not provide an image of the cell nuclei, whose size, shape and other characteristics can reveal abnormal cell activity.

The present invention rebates to systems and methods employing field-based Iight scattering spectroscopy to provide three-dimensional images and for determining physical characteristics of biological tissue. The three-dimensional images provide information regarding the characteristics 'of cellular or other structures including, for example, size, shape, texture of the cell nuclei as well as cellular dynamics which can include motility and chaos. Field-based LSS
systems incorporate spectroscopic information into a low coherence interferometer system.
This field-based type of LSS can be differentiated from intensity based LSS as it is able to isolate small regions of a biological sample for study, and yields information regarding the phase, as well as the amplitude, of the scattered field. Because field-based LSS measures phase of the wavefront, and wavefront uniformity, as well as providing a high degree of spatial resolution, field-based LSS can measure the shape and spatial distribution of cell nuclei along with size and refractive index of the nuclei. To ascertain the shape of cell nuclei, measurements of the scattered light's amplitude wavefront, phase wavefront, and wavefront non-uniformity axe combined.

Spatial distribution of cell nuclei is ascertained by three dimensional profiling which allows for imaging of cell nuclei density over the probed region.
In a preferred embodiment, the field-based LSS, system includes an angular light scattering system which measures both the position and angular distribution of a light field. Angular distributions provide information regarding variations in phase front uniformity (coherence of the scattered light) as well as indicate size of the scatterer. This system can include a light source such as a laser and an optical system that provides spatially coherent light having first and second wavelength components that are aligned to illuminate a region of interest of the material or tissue under study. A reference beam having both first and second wavelengths is also directed through the optical system along an optical path having an adjustable path length. An actuator is used to vary the path length at a selected speed such that detection of the reference beam and the light returning from the material under examination can generate a heterodyne signal which is and analyzed to determine selected characteristics of the material.
In another preferred embodiment, the field-based LSS system includes a dynamic light scattering system to measure dynamic properties 'of tissue structure such as cell nuclei within a small selected region by combining field-based LSS with temporal measurements. Motions of small objects can be detected by scattering of laser light through examination of its statistical properties such as, for example, autocorrelation function and power spectra.
A preferred embodiment of the invention involves a system for fiber optic delivery and a collection of light from tissue to provide real time diagnostic measurements either ih vitro or ih vivo. A preferred method of using the system involves analysis of the collected data by storing the data in memory of computer and comparing the collected data with reference data previously stored in electronic memory or database.
The foregoing and other objects, features and advantages of the invention will be apparent from the following more particular description of preferred embodiments of the invention, as illustrated in the accompanying drawings in which like reference characters refer to the same parts throughout the different views. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 illustrates a schematic diagram of a system using field-based light scattering spectroscopy in accordance with the present invention.
Figure 2 graphically illustrates a theoretical fit of the data points measured in accordance with the presentinvention.
Figure 3 illustrates a schematic diagram of a preferred embodiment of an angular light scattering system in accordance with the present invention.
Figures 4A and 4B graphically illustrate theoretical angular distributions in accordance with the present invention.
Figures SA and SB graphically illustrate path delay and intensity measurements related to the angular light scattering system in accordance with the present invention.
Figure 6 illustrates a block diagram of a preferred embodiment of a dynamic light scattering spectroscopy system in accordance with the present invention.
Figure 7 graphically illustrates. the correlation measurement related to the dynamic light scattering spectroscopy system in accordance with the present invention.
Figure 8 schematically illustrates a fiber optic system used in conjunction with the present inventions.
DETAILED DESCRIPTION OF THE INVENTION
Refernng to Figure l, a preferred embodiment of a field based LSS system 10 is illustrated, including a Michelson interferometer with two low-coherence light sources. In this particular embodiment, a Coherent MIR.A Tiaapphire laser 12 operating in femtosecond mode (150 fs) produces 800 nri1 light. The measured coherence length is about 30 ~,m. A portion of this light is handled by converter 15 such that it is split off and up-converted to 400 nm by means of a CSK
Optronics LBO second harmonic generation crystal. Note that two or more separate light sources can be used to provide the two or more wavelengths used for particular measurements. The converted light is then recombined with the original beam.
There is preferably substantial overlap (i.e., greater than 50% and preferably greater than 80°l°) between the two wavelength components. Reduced overlap will increase the scan time to illuminate a given surface are or volume. The combined beam 14 is then divided by a beam-sputter 16 into a probe beam 18 and a reference beam 20.
Alternatively, a single beam can be used with a variable wavelength source.
The probe beam 18 is focused onto the sample by means of a 12.7 mm focal length achromatic lens 22. The powers of the 400 and 800 nm components at the sample are 14 and 5.5 mW, respectively. The beam waists at the lens are 1.1 and 2.1 mm (FWI~VI), resulting in beam waists at the focal point of 5.9 and 6.2 ~uxn, and the Rayleigh lengths of 270 and 150 ~.m, respectively. As the coherence length is shorter than the Rayleigh lengths, the probe region can be approximately represented as a cylinder. The calculated beam overlap within this region is 98%, thus insuring that the same particles within the material 26 are illuminated simultaneously at both wavelengths.
The reference beam 20 is reflected from a mirror 24 moving along linear axis 27 at a constant speed, inducing Doppler shifts of 14.6 and 7.3 kHz at 400 and nm, respectively. It is then recombined with the probe light backscattered from the material 26 and transmitted to a dichroic mirror 30, which separates the 400 and 800 nm components and delivers them to separate auto-balanced photoreceivers 32,34 (New Focus 2007). Apertures 36, 38 are positioned in front of the detectors to limit the amount of collected light. In this particular embodiment the apertures have a radius of 0.9 mm. The distance from the collection lens 22 to the detector s is about 50 cm. A portion of the reference beam is separated with a second beamsplitter and sent to a second dichroic mirror 40. The two output components from the mirror are delivered to the reference ports 42, 44 of the photoreceivers. This serves to cancel power fluctuations at the two wavelengths.
The heterodyne signal at each wavelength, which results from the interference of the back-scattered probe beam and the appropriate Doppler-shifted reference beam, is detected by the photoreceiver. Each signal is measured using a Stanford Research 830 loclc-in amplifier array 46, 48. The magnitude of the heterodyne signal is displayed on a display 50 or oscilloscope and/or recorded by a computer 52.
In an example of a measurement of a sample in accordance with the invention, a cuvette with a thin layer of clear gelatin about 3 mW thicle is followed by a layer of polystyrene microspheres suspended in gelatin. The probe beam enters the sample through the first layer and is brought to a focus at the interface between the first and the second layers. The two layers have the same refractive index, which ensures that there is no reflective interface between them. Therefore, any heterodyne signal observed can be attributed solely to scattering from the microspheres.
The polystyrene microspheres range from 0.53 to 6.09 ~.m and can be obtained from Bangs Laboratories, Polysciences Inc. and Spherotech Inc. The volume concentration of microspheres within the second layer is maintained at 1.3%
for these samples. Based on an illumination volume of about 850 ~m3, as defined by the waist at the focus of the probe beam and its coherence length, the average number of microspheres illuminated ranges from 140 for the 0.535 ~,m microspheres to 0.1 for the 6.1 ~.m microspheres. The refractive index ratio of the microspheres to the gelatin is measured independently to be 1.19~0:01. To account for attenuation and absorption of the probe beam in the first layer of gelatin, the results are normalized using heterody~.ie signals from a cuvette in which the second layer of gelatin is replaced by a mirror.
The amplitude of the measured heterodyne signal depends primarily on three factors: the number of microspheres illuminated, their positions within the probe beam, and the amount of back scattering they produce. Errors that may arise from the first two factors include that they can potentially skew the measurement and mask the periodic structure of the light scattering spectrum. These contributions can be minimized or eliminated by aligning the 400 and 800 nm beam components and employing approximately equal beam waists at the focus, as described above.
This ensures that the beam components at both wavelengths illuminate the same region of the sample, and thus the same microspheres. Therefore, by taping the ratio of the _g_ signals at the two wavelengths, any dependence on the number and positions of the microspheres are normalized out.
In general, the heterodyne signal from the interference of the reflection from a group of scatterers with a reference beam can be expressed as:
~~2gr (P)ET (P)e"~T (a> ~P ~ xetowr ~ = [2Er (P)ET (P)~' ~T(a) ~P ~ (l a) with Er ( p ) the electric f eld of the reference beam incident on the detector, p the spatial coordinate on the detector's surface, ET ( p) is the amplitude component of the signal electric.field which is polarized in the same direction as the reference beam. Henceforth, all discussion of electric fields refer to those polarized in the plane of the reference beam's polarization. The symbol ~t~ denotes the heterodyne frequency, t represents time, and yr T(p) the associated phase shift of the signal field to the reference field. ~...~p indicates averaging over the detector area and L . .I denotes amplitude.
~[2Er~P~ET~P~~'vT(n)T~ ~_ ~~ [2E,.(P)E.i~t(P)e'W.'i,i(a)~ ~ (1b) J~° j-_1 t=1 Where yr 1,i ( p) denotes the amplitude of the scattered electric field associated with the ith trajectory which has undergone exactly j scattering events in the sample, and ~ ~,Z (p) the associated phase.
Each term on the right hand side of Eq. (1) can be rewritten as:
0 ~2Er~P)L',j,i(P)e"~sl~(P)~ -_aJ~ErrL''.l.i~'~''' ,F

where E,. and E~,i denote the root mean square averages of EY ( p) and E~,i ( p) over the detector areas, and ~r ~,i (p) the effective average phase shift. The phase and amplitude of the scattered field must be averaged over the detector since Iight incident on scatterers will be scattered in different amounts at different angles. The reference field is represented as having a Gaussian transverse spatial variation across the detector. The quantity a~,i is the spatial coherence factor. It measures the uniformity of y ~,Z (p) , as well as the correlation of E J,~ ( p) with Er ( p) . Its maximum value, 2, occurs when y~ ~,I (p) is constant over the detector surface, and E~,i ( p) is identical in profile to E,. (p) . Note that Er and E~,1 are simply the square roots of the average reference intensity and scattered intensity at the detector, respectively.
Note that a~,Z decreases very rapidly with increasing values of j. This indicates that a heterodyne measurement strongly favors the detection of singly scattered light over multiply scattered light, provided that the two are present in comparable amounts. Such is the case for light scattered near the surface of a turbid medium such as tissue. Thus, in our case the measured heterodyne signal can be approximated by the terms associated with the ~,i 's in Eq. (1b):
N
[2E,.(P)ET~P)e'~r(p)]pl NEr ~al°'El~ie~~n~.
t=i Note that since the ntunber of trajectories with only one scattering event equals the number of scatterers illuminated, N, the summation in Eq. (3) is terminated at N.
The amplitude of the heterodyne signal can be evaluated from Eq. (3) as:

H= [2E'r(P)E'~'(P)ea'l'T(p)~
,P
N _ N N _ _ _ _ E 'i~ (al'iEl'i)Z +i~ i~lal,ial,i~E1 iEl i~ cOS(~l,i - ~l,i') r i'~i _ N 2 '~ Er ~1(al,iEl,i) The cross terms in Eq. (4) are negligible if the yr l,i's are uncorrelated and the number of scatterers, N, in the probe region is large. The later requirement can be satisfied by averaging over a sufficient number of measurements taken at different points on the sample. The assumption that the yr l,i's are uncorrelated can be understood by noting that, in addition to the non-zero phase shift that the scattered light accumulates from the scatterers, there are additional random phase shifts due to the different round trip distances traversed by the light from scatterers at various depths. Given that the average spatial distance between neighboring scatterers in this measurement is at least l.~ ~,m or about 2.2 times the X00 nm wavelength, the scattered light contributions are unlikely to be correlated.
For a medium with microspheres of diameter D illuminated by light of wavelength ~,, we have an average ~ eld, E (~,, D) , given by the square root of the average intensity of light per scatterer reaching the detector, I (~,, D) .
The root mean square ensemble average of the heterodyne signal from Eq. (4) can then be written as:
_ N _ _ _ II (a,, D) ensernble-rats- Er ~ (al,i ~a'~ D, di )El i (~', D, di ))2 ehse~able - P392S
i0 _ ~I NEra (~,, D)Er (a,, D) _ Era (~,, D) I (a,, D) with d1 the displacement of the i'~' scatterer from the focal point. a (7~, D) is a defined quantity which accounts for the contributions of the a~,i (~,, D, dt )'s to the ensemble averaged heterodyne signal. This is a fundamental characterization of the scattering process. Note that x.11 quantities in Eq.(4) except E,. are functions of 7~, D
and di ; which is made explicit in Eq. (5).
In intensity-based LSS, N can be found by measuring the backscattered intensity at many wavelengths. More details regarding such measurements can be found in U.S. Application No. 08/948,734 filed on October 10, 1997. the entire contents of the referenced application being incorporated herein by reference.
The same can be done in field-based LSS. Unfortunately, where the heterodyne signals are measured at only two wavelengths, N cannot be accurately determined.
However, since the beam components at the two wavelengths are well aligned, and thus illuminate the same focal area, the same scatterers are simultaneously illuminated at both wavelengths in each measurement. To provide the ensemble average, malce M individual measurements of H, each at a different beam position and sum, then sum over all the measurements at each wavelength. The ratio, R, of the two sums then eliminates any dependence on N:
M
~,Hr2(~~D) N a(~~D) I(~~D) R = r=i M 2 a (~z ~ D) I(~a ~ D) ~Hr (°~a ~ D) r=~
Note that in addition to the scattered intensity, I (7~, D) , which determines the reflectance spectnun in a conventional LSS experiment, field-based LSS is also influenced by the wavelength-dependent spatial coherence factor, a (7~, D) .
This indicates that the field-based LSS spectrum differs from the intensity-based LSS
spectrum.
The ensemble averaged quantity a (~,, D) I (~,, D) can be calculated from Mie scattering theory for spherical objects. First express the contribution to the heterodyne signal from a single scatterer at a displacement dl from the focal point.
This is given by:
sta., D ~~ , B~ ) H(~, D, di ) _ C2Er (~, P)Ei (~, di ) k~, ~ , (7) P
with Et (7~, d~ ) the electric field strength at the scatterer and S(?~, D, ~i Bi , ) the amplitude function in the plane of polarization of the reference and incident beam.
(Note that, S(7~, D, ~i Bi , ) is a complex function.) In this notation, ~1 is the angle subtended from the vertical, defined as the direction of the probe beam's propagation, and ~~ is the angle in the plane normal to the vertical, with ~~
0 being along the direction of polarization. ~ ~ and ~t can be related to p and dZ by geometrical optics. k is the optical wavenumber and ~ the distance from the scatterer to the detector.
S(~,, D, ~i 9i , ) can, in turn, be expressed as:
S(~, D, ~t , ~e ) = sine (~~ )St (~, D, Bt ) + cost (~Z )Sz (~, D, Bt ), with Sl (7~, D, Bi , ) and Sa (7~, D, Bt , ) the amplitude functions given by Mie theory for in-plane and out-of plane scattering, relative to the scatterer.
From Eq. (4), the ensemble averaged quantity a (7~, D) I (~., D) equals the root mean squaxe ensemble average of the heterodyne signal associated with N--1. In other words, it can be expressed in terms of the single scatterer's heterodyne contribution:
a a _1j~j~a~ 1 a~~ I~~ "~~~~~)eruen:blerms ~~ all where Q is the number of elements in the ensemble average.
Using this formulation, the predictions for R for microspheres ranging from 0.1 to 7.0 ~,m in diameter can be determined. The result is then averaged over a 5%
variation in diameter, to account for the distribution of microsphere sizes in a given sample. This distribution is consistent with the size distribution given by the manufacturers.
Measured data are plotted in Figure 2. Each data point consists of 30 or 45 measurements taken at various sample positions. Averaging over scatterer position is determined by taking the root mean square of the measured heterodyne signal at each wavelength. The ration of the results, R, is then calculated and plotted.
The solid line in the plot is the theoretical fit. The fit was done using the procedure outlined above. In the fitting procedure, the aperture size of the detectors were free parameters. The best fit requires an aperture size of 0.7 mm in radius, instead of 0.9 mm, for the 400 nm measurements. This can be accounted for as a possible misalignment of the incoming beam at the detector, which can result in a smaller effective detection area.
The fit agrees well with the measurements for microspheres smaller than 5 microns in diameter. For the larger microsphere sizes the determination which is based on Mie scattering which assumes an incident plane wave field is not as accurate as the microsphere size becomes comparable to the beam waists (about microns).
While both field-based and intensity-based LSS determine the size and relative refractive index of scatterers by measuring variations in scattering across the spectrum, field-based LSS differs from intensity-based LSS in two important ways.
First, it is sensitive to phase front variations in the scattered wave, as well as to the intensity variations of the back-scattered light. Second, it permits greater localization of the region to be measured.
Single particle light scattering is characterized by the scattering amplitude, S(?~, D, ø, 8, ) , which is a complex function with a phase that varies with angular coordinates. Tlus phase is not measured in intensity-based LSS. In contrast, field-based LSS is very sensitive to variations in phase. To illustrate this, consider a measurement with a plane wave incident field. The average scattering intensity from a single scatterer measured at the detector can then be expressed as:
I(~,,D)= k j2 [,~S(~2,D,~,B)'Z~P

From Eq. (7), a similar field-based LSS experiment in which the incident and reference fields are both plane waves will give a heterodyne signal of the form:
z 4Er Etz . z D) = kz~z I [S(~~D~~~B)~pl . (11) Setting aside the average reference field intensity, EY a subtle but important difference between these two expressions can be seen. In Eq. (10), the magnitude of ~S'(~D,~,e) is taken before averaging over the detector area, whereas in Eq.
(I1), the sequence of operations is reversed. Thus, phase variations become important in field-based LSS, and the measured signal is proportional to the degree of phase-front uniformity of the scattered light. As Mie theory shows, the phase front uniformity degrades as a function of ~D) , and so for a scatterer of a given size, there is more variation across the spectrum with field-based LSS than with intensity-based LSS.
This, in turn, makes calculation of scatterer size easier and more sensitive.
The theory employs an ensemble-averaged spatial coherence factor, a (7~, D) , to quantify the degree of phase front uniformity. Note that 6(~,, D) is fundamentally related to f (7~, D, ~, 8) . By substituting Eqs. (10) and (11) into Eq.
(5) and setting N--1 (as the above calculation is for a single scatterer), provides:
~[s(~' D' ~' B)] ~ z a (7~, D) = 2 P (12) [I S('~~D~ ~~ B)~Z]P ' Note that Eq. (12) is applicable only in situations where both the reference and incident fields are uniform, and are appropriately modified for non-uniform input fields.
To clarify the physical significance of a (7~, D) , re-express Eqs. (10) and (11) in terms of the scattering cross section, a (~,, D) , and phase function, f(~.,D,~,e) [12]. The phase function is the normalized function which describes the angular intensity distribution of the scattered light. Since it is given by f (~~ D~ ~~ ~) = k26(~~ D) ~S(~~ D~ ~~ ~)~z express Eq. (10) as I (~~ D) = EZ ~'<~, D~[f (~~ ~~ ~~ e)JP (14) The corresponding equation for the heterodyne signal, Eq. (11), can then be written in terms of a (7~, D) as:
H(~~D)Z = EKE' ~C~~DJa'<~,DJ[f(~~D~~~B)JP (15) Tii intensity-based light scattering measurements, the scattering cross-section, 6 (~,, D) , and the anisotropy factor, g(7~, D) , are often used to characterize the scattering process. In contrast, specification of the scattering process in field-based light scattering requires a third parameter derived from S(~,, D, ~, B) the spatial coherence factor a (~,, D) .
In addition to this sensitivity to phase variations, field-based LSS has the advantage of being able to localize a much smaller region for study. This localization~can be achieved with low-coherence interferometric techniques, as exemplified by OCT. Spatial localization is limited by the coherence length, and a typical femtosecond laser began can localize a region for study on the order of tens of microns. This opens the possibility of selectively probing each individual scatterer, such as a single cell nucleus, even when it is surrounded by other scatterers.
The prospect of probing one scatterer at a time also simplifies the computations.
For a single scatterer, the cross-teens in Eq. (4) will not be present, thereby eliminating the need to average over numerous samples. We also note that with the high sensitivity afforded by heterodyne techniques, the full-spectral response of field-based LSS may provide a means to resolve features of a scatterer that are considerably smaller than an optical wavelength.
Figure 3 illustrates another preferred embodiment 60 of the invention in which a source of light 82, a beamsplitter 80 which splits the light along two paths, a first path extends through lens 74 onto the tissue or object being measured and the second path is directed through lens 66 reflects off mirror 64 and is redirected through the beamsplitter 80 along with light returning from the object being measured through lens 78 to detector system 76. In this embodiment both the lens 66 and mirror 64 are displaced over a distance 70 during the measurement.
Additionally, the lens 66 can be displaced in a second direction 68 to alter the direction of propagation. By scanning lens 66 and mirror 64 together as the angle of the reference field is scanned, the position of the reference field in the detector plan remains fixed. Figures 4A and 4B graphically illustrate the angular dependence of measurements made using the present invention. Figures SA and SB illustrate path delay and intensity measurement of a monolayer. The angular resolution of this measurement is about 0.45 mrad, and the longitudinal resolution is about 11.6 ~,m.
Figure 6 schematically illustrates a process in which measurements employing autocorrelation methods provide information and images of changes in cellular or subcellular structures in tissue. Figure 7 graphically illustrates the results of measurements of correlation in the intensity of light from a He-Ne laser scattered by a dilute suspension of beads. The autocorrelation time of 0.51 msec is consistent with multiple scattering of the light by 0.22 micron beads.
Illustrated in connection with Figure 8 is a fiber optic system 200 for light delivery and/or collection in conjunction with the light scattering spectroscopic systems and methods of the invention described previously. A light source provides a beam 202 that includes at least two wavelengths 7~1, ~,2 which are coupled to the proximal end of optical fiber 204. A beam splitter 206 incorporated into the fiber optic system delivers light components through fibers 208 and 210, and through lenses 216 and 214, respectively. A first light component is reflected by moving mirror 220 traveling in direction 220, and returns through fibers 210 and 212.
A
second light component is directed onto tissue 218, and light scattered by the tissue is returned through fibers 208 and 212. Dichroic mirror 230 separates the two wavelength ~,1 and ~,2 which are detected by detectors 240 and 242, respectively.
The heterodyne detection systems 250 and 252 are used to process the detected systems as described previously in connection with Figure 1. The systems described herein can be used in conjunction with standard endoscopics to provide diagnostic information retrieved from lumens or tissue within the human body ih ~ivo.
While this invention has been particularly shown and described with references to preferred embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the scope of the invention encompassed by the appended claims.

Claims (40)

What is claimed is:
1. A method of measuring a material comprising irradiating a region of interest in the material with spatially coherent light having a first beam with a first wavelength and a second beam with a second wavelength;
directing reference light having the first wavelength and the second wavelength along an optical path having, a variable path length;
detecting scattered light from the material in response to the irradiating light and detecting the reference light while varying the path length; and generating a heterodyne signal from the detected scattered light and the detected reference light.
2. The method of Claim 1 further comprising forming an image of the region of interest.
3. The method of Claim 1 further comprising measuring a size of a material within a region of tissue.
4. The method of Claim 1 wherein the first beam and the second beam irradiate a focal area within the region of interest.
5. The method of Claim 1 further comprising measuring the material at a plurality of first and second wavelengths.
6. The method of Claim 1 further comprising combining scattered light and the reference light and subsequently detecting the combined light.
7. The method of Claim 1 further comprising measuring a refractive index of a material within a region of tissue.
8. The method of Claim 1 further comprising recording data in electronic memory and comparing the data to reference data.
9. The method of Claim 1 further comprising using a fiber optic device to transmit light.
10. The method of Claim 1 further comprising using a low coherence light source.
11. The method of Claim 1 further comprising detecting backscattered light from a region of interest.
12. The method of Claim 1 further comprising adjusting a depth within the material being measured.
13. The method of Claim 1 further comprising aligning the first beam and the second beam to overlap at the region of interest.
14. An optical system for measuring a material comprising:
a light source and an optical system that irradiates a region of interest in a material with spatially coherent light having a first beam with a first wavelength and a second beam with a second wavelength;
a reference light beam having the first wavelength and the second wavelength along an optical path having a variable path length;
an actuator that adjusts the variable path length;
a detector system that detects scattered light from the material in response to the irradiating light and detects the reference light while varying the path length, the detector system generating a heterodyne signal from the detected scattered light and the detected reference light.
15. The system of Claim 14 further comprising a scanning assembly that scans the first beam and the second beam across the material such that an image of the region of interest.
16. The system of Claim 14 further comprising a data processor that computes a size of size tissue structure within the region of interest.
17. The system of Claim 14 wherein the first beam and the second beam irradiate a focal area within the region of interest.
18. The system of Claim 14 further comprising a light source emitting a plurality of first and second wavelengths.
19. The system of Claim 14 wherein the light source comprises a laser system that generates first and second wavelengths.
20. The system of Claim 14 further comprising a fiberoptic probe.
21. The system of Claim 14 wherein the light source comprises a wavelength tunable laser.
22. The system of Claim 14 further comprising a scanner that alters a beam path through the material.
23. The system of Claim 22 wherein the scanner alters an angle of the beam path relative to the material.
24. The system of Claim 14 further comprising a time correlation system.
25. The system of Claim 14 further comprising a fiber optic coupler and an interferometer.
26. The system of Claim 16 wherein the structure comprises a cell or a cell nucleus.
27. The system of Claim 14 further comprising a computer having a memory with stored reference data.
28. A method of analyzing a material comprising:
detecting scattered light from the material in response to irradiating light and detecting a reference light while varying a path length;
generating a heterodyne signal from the detected scattered light and the detected reference light; and comparing the heterodyne signal with reference data.
29. The method of Claim 28 further comprising forming an image of the region of interest.
30. The method of Claim 28 further comprising determining a size of a material within a region of tissue.
31. The method of Claim 28 wherein a first beam and a second beam irradiate a focal area within a region of interest of the material.
32. The method of Claim 28 further comprising detecting light from the material at a plurality of first and second wavelengths.
33. The method of Claim 28 further comprising combining scattered light and the reference light and subsequently detecting the combined light.
34. The method of Claim 28 further determining a refractive index of a material within a region of tissue.
35. The method of Claim 28 further comprising recording data in electronic memory and comparing the dat to reference data.
36. The method of Claim 28 further comprising using a fiber optic device to collect light.
37. The method of Claim 28 further comprising using a low coherence light source.
38. The method of Claim 28 further comprising detecting backscattered light from a region of interest.
39. The method of Claim 28 further comprising adjusting a depth within the material being measured.
40. The method of Claim 28 further comprising aligning the first beam and the second beam to overlap at the region of interest.
CA002377724A 2000-04-28 2001-04-27 Methods and systems using field-based light scattering spectroscopy Abandoned CA2377724A1 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US20018700P 2000-04-28 2000-04-28
US60/200,187 2000-04-28
PCT/US2001/040612 WO2001084124A2 (en) 2000-04-28 2001-04-27 Methods and systems using field-based light scattering spectroscopy

Publications (1)

Publication Number Publication Date
CA2377724A1 true CA2377724A1 (en) 2001-11-08

Family

ID=22740679

Family Applications (1)

Application Number Title Priority Date Filing Date
CA002377724A Abandoned CA2377724A1 (en) 2000-04-28 2001-04-27 Methods and systems using field-based light scattering spectroscopy

Country Status (8)

Country Link
US (1) US6847456B2 (en)
EP (1) EP1277040A2 (en)
JP (1) JP5134177B2 (en)
KR (1) KR20020027377A (en)
CN (1) CN1279345C (en)
AU (1) AU5980701A (en)
CA (1) CA2377724A1 (en)
WO (1) WO2001084124A2 (en)

Families Citing this family (45)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6611339B1 (en) * 2000-06-09 2003-08-26 Massachusetts Institute Of Technology Phase dispersive tomography
US7557929B2 (en) 2001-12-18 2009-07-07 Massachusetts Institute Of Technology Systems and methods for phase measurements
US6934035B2 (en) * 2001-12-18 2005-08-23 Massachusetts Institute Of Technology System and method for measuring optical distance
US7365858B2 (en) * 2001-12-18 2008-04-29 Massachusetts Institute Of Technology Systems and methods for phase measurements
FR2834567B1 (en) * 2002-01-07 2004-03-19 Centre Nat Rech Scient METHOD AND DEVICE FOR MICROSCOPIC VISUALIZATION WITH LOCAL PROBES OF A THREE-DIMENSIONAL OBJECT
ATE541202T1 (en) * 2002-01-24 2012-01-15 Gen Hospital Corp DEVICE AND METHOD FOR LOCATION AND REDUCTION OF NOISE OF SIGNALS IN LOW COHERENCE INTERFEROMETRY (LCI) AND OPTICAL COHERENCE TOMOGRAPHY (OCT) USING PARALLEL DETECTION OF SPECTRAL BANDS
US20040001662A1 (en) * 2002-06-27 2004-01-01 Biotechplex Corporation Method of and apparatus for measuring oscillatory motion
US6847449B2 (en) * 2002-11-27 2005-01-25 The United States Of America As Represented By The Secretary Of The Navy Method and apparatus for reducing speckle in optical coherence tomography images
US20060256342A1 (en) * 2003-03-28 2006-11-16 Wong Lid B Confocal microscope system for real-time simultaneous temporal measurements of metachronal wave period and ciliary beat frequency
US7102758B2 (en) 2003-05-06 2006-09-05 Duke University Fourier domain low-coherence interferometry for light scattering spectroscopy apparatus and method
US7006219B2 (en) * 2003-05-20 2006-02-28 Technology Innovations, Llc Biological imager
US6861972B2 (en) * 2003-07-28 2005-03-01 Ellistar Sensor Systems, Inc. Object detection apparatus and method
EP1655592A4 (en) * 2003-08-12 2007-08-29 Bussan Nanotech Res Inst Inc Detection device, optical path length measurement device, measurement instrument, optical member evaluation method, and temperature change detection method
US8634072B2 (en) * 2004-03-06 2014-01-21 Michael Trainer Methods and apparatus for determining characteristics of particles
US8705040B2 (en) * 2004-03-06 2014-04-22 Michael Trainer Methods and apparatus for determining particle characteristics by measuring scattered light
US10955327B2 (en) * 2004-03-06 2021-03-23 Michael Trainer Method and apparatus for determining particle characteristics utilizing a plurality of beam splitting functions and correction of scattered light
ATE421077T1 (en) * 2004-10-29 2009-01-15 Degudent Gmbh METHOD AND DEVICE FOR COLLECTING CONTOUR DATA AND/OR OPTICAL PROPERTIES OF A THREE-DIMENSIONAL SEMIT-TRANSPARENT OBJECT
US7385460B1 (en) 2004-11-17 2008-06-10 California Institute Of Technology Combined electrostatic and optical waveguide based microfluidic chip systems and methods
US7428048B1 (en) 2004-12-30 2008-09-23 Spectral Molecular Imaging Inc. Imaging elastic scattering spectroscopy
US20080119740A1 (en) * 2004-12-30 2008-05-22 Iddan Gavriel J Device, System, and Method for Optical In-Vivo Analysis
US20080204716A1 (en) * 2005-03-07 2008-08-28 Michael Trainer Methods and apparatus for determining characteristics of particles
EP2950065A1 (en) * 2005-10-11 2015-12-02 Duke University Method for fiber-based endoscopic angle-resolved low coherence interferometry
US8537366B2 (en) 2005-10-11 2013-09-17 Duke University Systems and methods for endoscopic angle-resolved low coherence interferometry
US7531774B2 (en) * 2006-06-05 2009-05-12 General Dynamics Advanced Information Systems, Inc. Measurement-diverse imaging and wavefront sensing with amplitude and phase estimation
CA2658481A1 (en) * 2006-07-21 2008-01-24 Oncoscope, Inc. Protective probe tip, particularly for use on a fiber-optic probe used in an endoscopic application
WO2008115089A1 (en) * 2007-03-16 2008-09-25 Germanov, Evgeny Pavlovich Diagnosis method using ultrasonic, acoustic and electromagnetic waves
US7781733B2 (en) * 2007-05-16 2010-08-24 International Business Machines Corporation In-situ high-resolution light-optical channel for optical viewing and surface processing in parallel with charged particle (FIB and SEM) techniques
US7920269B2 (en) * 2007-08-03 2011-04-05 Chung Yuan Christian University System and method for measuring interferences
WO2009036418A1 (en) * 2007-09-13 2009-03-19 Duke University Apparatuses, systems, and methods for low-coherence interferometry (lci)
CA2711643A1 (en) * 2008-01-08 2009-07-16 Oncoscope, Inc. Systems and methods for tissue examination, diagnostic, treatment, and/or monitoring
US7772569B2 (en) * 2008-04-01 2010-08-10 The Jackson Laboratory 3D biplane microscopy
US8662962B2 (en) * 2008-06-30 2014-03-04 3M Innovative Properties Company Sandpaper with non-slip coating layer and method of using
US9823127B2 (en) 2010-01-22 2017-11-21 Duke University Systems and methods for deep spectroscopic imaging of biological samples with use of an interferometer and spectrometer
AU2011207444A1 (en) 2010-01-22 2012-08-09 Duke University Multiple window processing schemes for spectroscopic optical coherence tomography (OCT) and fourier domain low coherence interferometry
CN101788474B (en) * 2010-02-10 2011-11-23 哈尔滨工业大学 Method for quantitatively testing threshold value of photo-induced light scattering exposure energy density of photorefractive materials
CA2903201C (en) 2013-03-14 2023-02-21 Research Development Foundation Apparatus and methods for optical coherence tomography and two-photon luminescence imaging
US9151671B2 (en) * 2013-08-26 2015-10-06 Thermo Electron Scientific Instruments Llc Motorized variable path length cell for spectroscopy
CN104122228B (en) * 2014-07-14 2016-06-15 江苏大学 A kind of Integrated Light interferes the micro imaging system and method analyzed with scattered information
US9633277B2 (en) 2014-09-12 2017-04-25 Research Development Foundation Apparatus and methods for identifying and evaluating bright spot indications observed through optical coherence tomography
JP5859154B1 (en) 2015-03-06 2016-02-10 リオン株式会社 Particle counter
WO2017136773A1 (en) * 2016-02-04 2017-08-10 Canon U.S.A., Inc. Diffusing wave spectroscopy with heterodyne detection
CN106680159B (en) * 2017-01-13 2023-05-23 浙江大学 Particle rapid detection method and system based on dynamic light scattering sample ensemble analysis
WO2019183291A1 (en) * 2018-03-20 2019-09-26 Cornell University Aberration diverse optical coherence tomography (oct) imaging to suppress optical scattering noise
CN110646377B (en) * 2019-09-27 2022-07-05 京东方科技集团股份有限公司 Blood detection system and detection method thereof
US11506478B2 (en) * 2019-11-27 2022-11-22 California Institute Of Technology Three dimensional (3D) imaging using optical coherence factor (OCF)

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4080264A (en) 1976-03-01 1978-03-21 Massachusetts Institute Of Technology Immunoassay by light scattering spectroscopy
US5481360A (en) * 1992-06-19 1996-01-02 Citizen Watch Co., Ltd. Optical device for measuring surface shape
US5410413A (en) 1993-08-18 1995-04-25 Petrometrix Ltd. Optical head probe using a gradient index lens and optical fibers
DE4411017C2 (en) 1994-03-30 1995-06-08 Alexander Dr Knuettel Optical stationary spectroscopic imaging in strongly scattering objects through special light focusing and signal detection of light of different wavelengths
US5627666A (en) 1994-07-27 1997-05-06 Board Of Regents Of The University Of Colorado Liquid crystal phase modulator using cholesteric circular polarizers
US5919140A (en) 1995-02-21 1999-07-06 Massachusetts Institute Of Technology Optical imaging using time gated scattered light
US5838485A (en) * 1996-08-20 1998-11-17 Zygo Corporation Superheterodyne interferometer and method for compensating the refractive index of air using electronic frequency multiplication
US6015969A (en) * 1996-09-16 2000-01-18 The Regents Of The University Of California Multiple-wavelength spectroscopic quantitation of light-absorbing species in scattering media
US5760901A (en) 1997-01-28 1998-06-02 Zetetic Institute Method and apparatus for confocal interference microscopy with background amplitude reduction and compensation
JP2001515382A (en) 1997-03-06 2001-09-18 マサチューセッツ インスティチュート オブ テクノロジー Equipment for optical scanning of living tissue
US5943129A (en) 1997-08-07 1999-08-24 Cambridge Research & Instrumentation Inc. Fluorescence imaging system
JP3999437B2 (en) 2000-03-10 2007-10-31 富士フイルム株式会社 Optical tomographic imaging system

Also Published As

Publication number Publication date
JP5134177B2 (en) 2013-01-30
AU5980701A (en) 2001-11-12
US20020101593A1 (en) 2002-08-01
EP1277040A2 (en) 2003-01-22
JP2003532112A (en) 2003-10-28
KR20020027377A (en) 2002-04-13
CN1451094A (en) 2003-10-22
WO2001084124A3 (en) 2002-06-20
WO2001084124A2 (en) 2001-11-08
US6847456B2 (en) 2005-01-25
CN1279345C (en) 2006-10-11

Similar Documents

Publication Publication Date Title
US6847456B2 (en) Methods and systems using field-based light scattering spectroscopy
US6611339B1 (en) Phase dispersive tomography
JP5555277B2 (en) System and method for angle resolved low coherence interferometry with endoscope
US6853457B2 (en) Optical amplification in coherence reflectometry
USRE46412E1 (en) Methods and systems for performing angle-resolved Fourier-domain optical coherence tomography
USRE42641E1 (en) Depth-resolved spectroscopic optical coherence tomography
US6900943B2 (en) Optical amplification in coherent optical frequency modulated continuous wave reflectometry
US20090219544A1 (en) Systems, methods and computer-accessible medium for providing spectral-domain optical coherence phase microscopy for cell and deep tissue imaging
Mitsui Dynamic range of optical reflectometry with spectral interferometry
Schmitt et al. Optical characterization of dense tissues using low-coherence interferometry
US6256102B1 (en) Dual-beam low-coherence interferometer with improved signal-to-noise ratio
JP2000046729A (en) Apparatus and method for high-speed measurement of optical topographic image by using wavelength dispersion
CN108572161A (en) Optical coherence tomography based on partial wave front interferometer
JP3688608B2 (en) Optical coherence tomographic image measuring device with spectroscopic function
Yang et al. Feasibility of field-based light scattering spectroscopy
EP1972926A1 (en) Assessment of parameters by optical measurement
Schmitt et al. Interferometric versus confocal techniques for imaging microstructures in turbid biological media
Xu et al. Polarization Imaging of Back-Scattered Terahertz Speckle Fields
Gao et al. Influence of multiple scattering on the image quality obtained with optical coherence microscopy
Poupardin et al. Study of coherent reflectometer for imaging internal structures of highly scattering media
Graf et al. Development of a pre-clinical Fourier domain low coherence interferometry system
Wax et al. Structure and organization of cellular components measured in neoplastic tissues using angular low coherence interferometry
Piao Flow velocity estimation in optical Doppler tomography and a preliminary study on radiation detection for hybrid optical coherence tomography/scintigraphy

Legal Events

Date Code Title Description
EEER Examination request
FZDE Discontinued