The present invention relates to apparatus, method, logic arrangement and storage medium for dramatically increasing the sensitivity in the detection of optical coherence tomography and low coherence interferometry signals by detecting a parallel set of spectral bands, each band being a unique combination of optical frequencies.
Two methods currently exist to implement depth ranging in turbid media The first method is known as Low Coherence Interferometry (“LCI”). This method uses a scanning system to vary the reference arm length and acquire the interference signal at a detector and demodulating the fringe pattern to obtain the coherence envelope of the source cross correlation function. Optical coherence tomography (“OCT”) is a means for obtaining a two-dimensional image using LCI. OCT is described by Swanson et al. in U.S. Pat. No. 5,321,501. Multiple variations on OCT have been patented, but many suffer from less than optimal signal to noise ratio (“SNR”), resulting in non-optimal resolution, low imaging frame rates, and poor depth of penetration. Power usage is a factor in such imaging techniques. For example in ophthalmic uses, only a certain number of milliwatts of power is tolerable before thermal damage can occur. Thus, boosting power is not feasible to increase SNR in such environments. It would be desirable to have a method of raising the SNR without appreciably increasing power requirements.
A second method for depth ranging in turbid media is known in the literature as spectral radar. In spectral radar the real part of the cross spectral density of sample and reference arm light is measured with a spectrometer. Depth profile information is encoded on the cross-spectral density modulation. Prior designs for spectral radar is primarily found in the literature.
The use of spectral radar concepts to increase the signal to noise ratio of LCI and OCT have been described earlier. However, in this description, only the real part of the complex spectral density is measured and the method uses a large number of detector elements (about 2,000) to reach scan ranges on the order of a millimeter. It would be desirable to have a method that would allow for an arbitrary number of detector elements. Secondly, the previously described method uses a single charge coupled device (“CCD”) to acquire the data. Since the charge storage capacity is limited, it requires a reduction of the reference arm power to approximately the same level as the sample arm power, giving rise to auto correlation noise on the sample arm light. In addition, since no carrier is generated, the 1/f noise will dominate the noise in this system. Thirdly, even with the short integration times of state of the art CCD technology, phase instabilities in the interferometer reduce fringe visibility of the cross spectral density modulation.
The present invention can increase the SNR of LCI and OCT by splitting the LCI broad bandwidth source into a number “N” of spectral bands. In one exemplary embodiment, the N spectral bands are individually detected and processed to provide an increase in the SNR by a factor of N. This increase in SNR enables LCI or OCT imaging by a factor of N times faster, or alternatively allows imaging at the same speed with a source that has N times lower power. As a result, the present invention overcomes two of the most important shortcomings of conventional LCI and OCT, namely, source availability and scan speed. The factor N may reach more than 1,000, and allows construction of OCT and LCI systems that can be more than three orders of magnitude improved from OCT and LCI technology currently in practice.
The present invention improves current data acquisition speeds and availability of sources for OCT. Shot noise is due to the statistical fluctuations of the current that are due to the quantized or discrete electric charges. The reduction of shot noise allows for much lower source powers or much higher acquisition rates. Limitations in current data acquisition rates (approximately 4 frames/sec) are imposed by available source power and availability of fast mechanisms for scanning delay. An increase in the sensitivity of the detection by a factor of 8 would allow real time imaging at a speed of about 30 frames per second. An increase of the sensitivity by a factor of about 1,000-2,000 would allow for the use of sources with much lower powers and higher spectral bandwidths which are readily available, cheaper to produce, and can generate higher resolution LCI or OCT scans.
For ophthalmic applications of OCT, the efficient detection preferably allows for a significant increase of acquisition speed. The limitation in ophthalmic applications is the power that is allowed to enter the eye according to the ANSI standards (approximately 700 microwatts at 830 nm). Current data acquisition speed in ophthalmic applications is approximately 100-500 A-lines per second. The power efficient detection would allow for A-line acquisition rates on the order of about 100,000 A-lines per second, or video rate imaging at about 3,000 A-lines per image.
The gain in SNR is achieved because the shot noise has a white noise spectrum. An intensity present at the detector at frequency ω (or wavelength λ) contributes only to the signal at frequency ω, but the shot noise is generated at all frequencies. By narrowing the optical band width per detector, the shot noise contribution at each frequency can be reduced, while the signal component remains the same.
In summary, the present invention improves a performance of LCI and OCT, and as a result can be used in developing LCI and OCT diagnostic technologies for medical and non-medical applications.
Other features and advantages of the present invention will become apparent upon reading the following detailed description of embodiments of the invention, when taken in conjunction with the appended claims.
The invention is illustrated in the drawings in which like reference characters designate the same or similar parts throughout the figures of which:
Overview
Certain exemplary embodiments of the present invention include a hybrid method that implements aspects of LCI and OCT where the reference arm is scanned, and spectral radar, which does not require reference arm scanning.
In one embodiment, the signal in the detection arm of an OCT system is split into more than one spectral band before detection. Each spectral band is detected by a separate photodetector and amplified. For each spectral band, the signal can be band pass filtered around the signal band by analog electronics and digitized, or, alternatively, the signal may be digitized and band pass filtered in software. As a consequence, the shot noise contribution to the signal can be reduced by a factor equal to the number of spectral bands, while output of the signal remains the same. The reduction of the shot noise increases the dynamic range and sensitivity of the system.
In another exemplary embodiment of the present invention, an apparatus is provided for spectral radar that does not require reference arm scanning. For many detectors, no ranging or reference arm scanning is needed, and the method may be similar to the method which can be employed for a spectral radar except that phase information of the cross spectral density is preferably preserved.
In other exemplary embodiments, the present invention describes an arrangement for spectral radar that eliminate phase instability in the interferometer, obtaining the complex spectral density and eliminating auto correlation noise on the sample arm light, relative intensity noise, and 1/f noise.
Theory
Time Domain Versus Spectral Domain OCT
Nearly all conventional OCT systems are based on Time Domain scanning. In such conventional systems, the length of the reference arm in a Michelson interferometer is rapidly scanned over a distance corresponding to the imaging depth range. An alternative procedure to scanning the reference arm, is one that measures the cross-spectral density at the detection arm of the Michelson interferometer using a spectrometer. In Spectral Domain OCT, no mechanical (e.g., motionless) scanning of the reference arm is required, while an apparatus for generating a phase shift can be used. Only recently was it recognized that a significant signal to noise gain can be achieved by direct measurement of the cross-spectral density.
Certain exemplary embodiments of the present invention provide a detection principle based on Spectral Radar concepts (further referred to as Spectral Domain OCT) or a hybrid method between Spectral Domain and Time Domain OCT that can be more sensitive than current state of the art Time Domain OCT, allowing a substantial increase in the acquisition speed to resolution ratio.
Principle of Shot Noise Reduction in Spectral Domain OCT
The best signal to noise performance of Time Domain OCT systems is obtained when the noise is shot noise limited. Shot noise can be reduced significantly by replacing the single element detector with a multi-element array detector. When the detection arm light is spectrally dispersed on the array detector, each element of the array detects a small wavelength fraction of the spectral width of the source. The shot noise is preferably reduced by a factor equal to the number of elements of the array. The principle of the signal to noise improvement is based on the white noise characteristic of shot noise and the observation that only electromagnetic waves of the same wavelength produce interference fringes.
The shot noise power density Nshot (f) (in units [W/Hz], [A2/Hz] or [V2/Hz]) is proportional to the current (or equivalently the optical power times the quantum efficiency) generated in the detector. For a monochromatic beam of wavelength λ1 entering the interferometer, the fringe frequency or carrier f at the detector is determined by the velocity ν of the mirror, f1=2ν/λ1. The shot noise is proportional to the power (or spectral density S(ω)) at wavelength λ1. A second wavelength λ2 is preferably coupled into the interferometer. A second fringe frequency or carrier at frequency f2=2ν/λ2 is simultaneously present. The shot noise at this second frequency is preferably the sum of the shot noise generated by the optical power at wavelength λ1 and λ2. Also, at frequency f1 the shot noise is the sum of the shot noise generated by the optical power at wavelength λ1 and λ2. Thus, at both frequencies a cross-shot noise term is generated by the simultaneous presence of both wavelengths at the detector. By spectrally dispersing each wavelength to a separate detector, the cross shot noise term can be eliminated. In this way, Spectral Domain OCT offers a significant improvement of signal to noise ratio over Time Domain OCT systems.
Signal to Noise Analysis of Time Domain Versus Spectral Domain OCT.
Signal
Analysis of the Signal to Noise Ratio (SNR) in Time Domain OCT has been described in related publications. The interference fringe peak amplitude in time domain OCT is given by
Ipeak=√{square root over (PrefPsample)}, (1)
with Pref, Psample the reference and sample arm power in Watts, respectively. In terms of electrical power at the detector, the signal in units [A2] is defined as
S=η2e2PrefPsample/Eν2, (2)
with η the quantum efficiency, e the charge quantum and Eν=hc/λ the photon energy. The reference and sample arm powers are given by the respective reflected spectral densities,
Pref,sample=∫Sref,sample(ω)dω. (3)
Assuming that the reference and sample spectral densities are equal to the source spectral density S(ω), where the sample arm spectral density is attenuated by a large factor, i.e., Sref(ω)=S(ω), Ssample(ω)=αS(ω) with α<<1, and inserting the above expression of reference and sample arm into the original definition of the signal gives,
S=η2e2α[∫S(ω)dω]2/Eν2. (4)
Thermal, Shot Noise and Relative Intensity Noise Contributions
Three contributions to the total noise of OCT signals are: thermal noise, shot noise and relative intensity noise. Thermal noise is generated by the feedback resistor, shot noise is related to the finite nature of the charge quantum resulting in statistical fluctuations on the current, and relative intensity noise is related to the temporal fluctuations due to chaotic character of classical light sources. These three contributions to the noise density in units [A2/Hz] are given by,
k is Boltzmann's constant, T the temperature in Kelvin, Rfb the value of the feedback resistor, and τcoh the coherence time of the source. Coherence time is related to the full spectral width at half maximum Δλ of a Gaussian source by the following relation, τcoh=√{square root over (2ln2/π)}λ02/(cΔλ). Shot noise limited detection is achieved when the second term in Eq. (5) dominates the other noise contributions.
Signal to Noise Ratio (SNR)
The signal to noise ratio (SNR) is given by
with BW the signal bandwidth, and parameters S and Nnoise(f) as described above.
Space and Frequency Domain Description of the OCT Signal
The OCT signal is most easily described in the space domain. For a single object in the sample arm, the interference term of the OCT signal is proportional to the real part of the Fourier transform of the source spectrum S(ω),
I(Δz)∝Re∫exp(ikΔz)S(k)dk, (7)
with Δz the path length difference between sample and reference arm and k the wave vector. As a function of time, the OCT signal is given by,
I(t)∝Re∫exp(2iωtν/c)S(ω)dω, (8)
with v the reference arm mirror velocity. The frequency spectrum of the signal is given by a Fourier transform of the signal in the time domain, resulting in a complex function. The absolute value of this function is equal to the spectral density,
|I(f)|=|∫I(t)e2iπftdt|=S(πfc/ν), (9)
which shows that the signal bandwidth is directly proportional to the source spectral width and scales linearly with the reference arm mirror velocity, i.e., imaging speed. Eq. (9) also directly relates the absolute value of the frequency spectrum, |I(f)|, to the signal S (Eq. (4)).
Eq (9) also demonstrates that each angular frequency of the light source or equivalently each wavelength of the source is represented at its own frequency in the measured interferometric signal. The depth profile information I(t) can be obtained from the complex cross spectral density I(f) by a Fourier transform.
The complex cross spectral density can also be obtained by splitting the signal I(t) in several spectral bands using a dispersive or interferometric element. At each detector, only part of the complex cross spectral density is determined. Combining the cross spectral densities of each detector, the full spectral density of the signal is retrieved.
Thus, the same information can be obtained by separating spectral components to individual detectors. Combining the signal of all detectors in software or hardware would result in the same signal as obtained with a single detector.
Signal to Noise Gain with Spectral Domain OCT
In the detection arm, the spectrum can be split into two equal halves, where two detectors each detect one half of the spectrum. According to Eq (9), the frequency spectra at detectors 1 and 2 are given by |I1(f)|=S(πfc/ν) for f<f0, I1(f)=0 for f>f0 and I2(f)=0 for f<f0, |I2(f)|=S(πfc/ν) for f>f0, respectively. The frequency spectrum as would be acquired by a single detect6r in time domain OCT is given by the sum of I1(f) and I2(f); I(f)=I1(f)+I2(f). Thus, the signal S after combining the spectra is equal, however I1(f)=0 for f>f0 and I2(f)=0 for f<f0, the bandwidth BW per detector can be reduced by a factor of 2.
The noise is determined by the sum of the shot noise contributions at detectors one and two. From Eqs. (5) and (6), the shot noise per detector is proportional to the reference arm power at the detector times the bandwidth for the detector. Since the spectrum was split in equal halves, the reference power at detectors 1 and 2 is, respectively,
Pref1=0.5Pref, Pref2=0.5 Pref. (10)
The sum of the shot noise contribution for the two detectors is,
NnoiseSD∝Pref1×0.5BW+Pref2×0.5BW=0.5PrefBW, (11)
which may compared with the shot noise of a single detector in time domain OCT,
NnoiseTD∝PrefBW. (12)
Thus, by spectrally dispersing the detection arm light over two separate detectors, the signal remains the same, while the noise is reduced by a factor of 2, resulting in a net
SNR gain by a factor of 2.
Extending the above analysis, it can be demonstrated that the shot noise contribution is reduced by a factor equal to the number of detectors. The sum of shot noises for N detector elements, where each detector element receives one Nth of the total reference power, is given by,
The signal is the same as in Time Domain OCT, and the SNR ratio for Spectral Domain OCT is given by,
Thus Spectral Domain OCT enables a SNR improvement over Time Domain OCT of a hundred to a thousand fold, depending on the number of detector elements N. Using a charge coupled array or an integrating device as a detector, such as, but not limited to, a line scan camera, the ratio N/BW is replaced by the integration time τi of the array, which results in,
Advantages
The exemplary embodiment of the present invention reduce shot noise and other forms of noise which allows for much lower source powers, or much higher acquisition rates than current systems. The increased detection sensitivity allows for real time imaging. Such imaging speed can help practitioners where motion artifacts are a continuing problem, such as in gastrointestinal, ophthalmic and arterial imaging environments. By increasing the frame rate while maintaining or improving the signal to noise ratio such artifacts can be minimized. The present invention also enable one to screen large areas of tissues with OCT and allows clinical viable screening protocols using this method.
Sources
The source arm 203 contains at least light source 202 that is used to illuminate the interferometer with low-coherence light. The source temporal coherence length is preferably shorter than a few microns (a preferred range is about 0.5 μm-30 μm). Examples of sources include, but are not limited to, semiconductor optical amplifier, superluminescent diodes, light-emitting diodes, solid-state femtosecond sources, amplified spontaneous emission, continuum sources, thermal sources, combinations thereof and the like. Other appropriate sources known to those skilled in the art may be used. While light is referred to herein as the source, it is intended that other electromagnetic radiation ranges may be suitable for use, depending on the circumstances.
Interferometer
The sample arm 208 collects light reflected from the tissue sample 130 and is combined with the light from the reference arm 206 to form interference fringes. The reference arm 206 returns light back to be combined with the source arm 203. The reference arm can also be transmissive with no reflection. This action of beam splitting/recombining may be performed using a beam splitter 204 (Michelson), or circulator(s) (Mach-Zehnder) or other means known to those skilled in the art for separating a beam into multiple paths and recombining these multiple beams in a manner that interference between the beams may be detected. The splitting may be accomplished in free space or by using a splitter 204 having passive fiber optic or waveguide components.
Sample Arm
For LCI applications, the sample arm may be terminated by an optical probe comprising a cleaved (angled, flat, or polished) optical fiber or free space beam. A lens (such as, but not limited to, aspherical, gradient index, spherical, diffractive, ball, drum or the like) may be used to focus the beam on or within the sample. Beam directing elements (such as, but not limited to, mirror, prism, diffractive optical element or the like) may also be contained within the probe to direct the focused beam to a desired position on the sample. For OCT applications, the position of the beam may be changed on the sample as a function of time, allowing reconstruction of a two-dimensional image. Altering the position of the focused beam on the sample may be accomplished by a scanning mirror (such as, but not limited to, a galvanometer, piezoelectric actuator or the like), electrooptic actuator, or moving the optical fiber (for example, rotating the optical fiber, or linearly translating the optical fiber). The sample arm probe may be a fiber optic probe that has an internally moving element where the motion is initiated at a proximal end of the probe and the motion is conveyed by a motion transducing arrangement (such as, but not limited to, wire, guidewire, speedometer cable, spring, optical fiber and the like) to the distal end. The fiber optic probe may be enclosed in a stationary sheath which is optically transparent where the light exits the probe at the distal end.
Reference arm delay
A mechanism 270 in the reference arm 206 allows for scanning the group delay of the reference arm 206. This group delay can be produced by any of a number of techniques known to those skilled in the art, such as, but not limited to, stretching an optical fiber, free space translational scanning using a piezoelectric transducer, or via a grating based pulse shaping optical delay line. Preferably, the delay is introduced by a non-mechanical or motionless arrangement. By “non-mechanical” it is meant that no mechanically moving parts are utilized. The absence of mechanically moving parts is believed to reduce the known deficiencies of using mechanical devices to introduce delay. As opposed to traditional LCI or OCT systems described in the literature, the reference arm 206 in the present invention does not necessarily need to scan over the full ranging depth in the sample, and preferably scans over at least a fraction of the ranging depth equal to one over the number of detectors (1/N). This scanning feature is fundamentally different from known delay scanning schemes used in conventional known LCI and OCT systems. The reference arm 206 optionally has a phase modulator mechanism (described more fully herein), such as, but not limited to, an acoustooptic modulator, electrooptic phase modulator or the like, for generating a carrier frequency. In order to reduce the scan range of the reference arm 206, the spectrum is preferably split into a plurality of spectral bands according to a method that will be explained below.
Detection
Referring to
Detector signals can be amplified by Trans Impedance Amplifiers (“TIAs”) 116, band pass filters 124 (digitally or using analog circuitry) and digitized by A/D converters and stored in a computer 122 for further processing. Each detector 114 is preferably configured to be shot noise limited. Shot noise limited detection is preferably achieved by adjusting the intensity of light returned from the reference arm 108 so that the shot noise dominates over the thermal noise of the resistor in the TIA 116 and is higher than the relative intensity noise (“RIN”). Each detector 114 is balanced for such dual noise reduction.
In one embodiment of the present invention, the number of detectors 114, N can be in the range of 2-10,000 or more. A preferred range of N is about 8-10,000 detectors. In one preferred embodiment, eight detectors 114 (or a number in that area) can provide real time, or close to real time, imaging.
Alternatively, another way for detection includes an integrating one-dimensional or two-dimensional detector 114 array which is capable of obtaining images at a rate preferably greater than 1/f noise (f=frequency) (see
This system could be implemented using a single detector 114 with dual-balanced detection enabled by either interleaving dual balanced rows of the array detector or by placing two similar array detectors adjacent to one another. If two array detectors 114 and 115 are used, the values are subtracted from one another to achieve dual balance detection. If more than two array detectors are used the signals can be selectively subtracted and complex spectral density can be obtained.
The spectral intensity as a function of wavelength is preferably constant. However, if it is not, the spectrum can be shaped in the reference, sample and/or source arms to make it constant. Spectral shapers are known in the art.
Processing
The signal of each detector 114 is band pass filtered around the signal frequency, such as by FFT's. The signal of all detectors 114 can be combined as explained hereinabove to obtain the complex cross spectral density in the frequency domain. By Fourier transform, the complex cross spectral density can be converted to a depth profile in the tissue. Several methods to process the complex spectral density to obtain depth profile information are known to those skilled in the art, such as, but not limited to, by obtaining at least two signals with a pi/2 phase shift in the reference arm and then reconstructing the complex spectral density by some linear combination of the two signals.
Following detection analog processing includes a trans impedance amplifier, band pass filter, and digitization of the signal. This signal may then be converted to reflectivity as a function of depth by the Fourier transform operation. Digital processing includes digitization, digital band pass filtering in either the frequency domain or time domain (FIR or IIR filter) and inverse Fourier transformation to recover the tissue reflectivity as a function of depth.
System Integration
Processing of the multiple signals may be performed using an imaging or diagnostic console which performs basic operations including, mathematical image reconstruction, display, data storage. Alternatively, another embodiment, shown in
Scan Range of the Reference Arm.
The ranging depth in the sample 130 is determined by the resolution with which the cross spectral density can be determined. In a method using a single detector the spectral resolution of the complex spectral density is determined by the scan range of the reference arm. The larger the scan range, the higher the spectral resolution and the larger the ranging depth in the sample. In a system with a spectral separating unit and multiple detectors, the resolution of the cross spectral density is a combination of reference arm scan range and spectral separating characteristics.
Any suitable wavelength band shape may be used for separating. For arbitrary spectral band shapes, the scan range of the reference arm 18 is determined by the delay that is needed to completely resolve the spectral components in each band.
For instance, in one preferred embodiment, as depicted in
Embodiments of the Wavelength Separating Filter
Several techniques are known to separate or disperse the spectrum. One method uses a grating and a micro lens array to focus spectral components onto individual detectors. A second method uses prisms instead of a grating. A third method uses a grating and an addressable mirror array (such as, but not limited to, a “MEMS” mirror or digital light processing (“DLP”) apparatus or the like) to direct spectral components to individual detectors. A fourth method uses a linear array of optical filters prior to the array of individual detectors. A fifth method uses waveguides etched into a material or manufactured from fiber optic components to generate a pattern with the desired filter action. As an example, in
Relative Intensity Noise
One of the noise terms that is present at the detectors is relative intensity noise (“RIN”) or Bose-Einstein noise. RIN noise likely becomes dominant over shot noise for spectral widths less than a few nanometers. For many detector configurations, the spectral width at each detector may likely be smaller than a few nanometers, and the relative intensity noise can dominate the overall system noise. Thus, balanced detection, can preferably be implemented to eliminate the RIN. Several methods known in the art exist to implement balanced detection. One such method will be discussed below in further detail. For example, but not by way of limitation, as shown in
An example of such balanced detection is shown in
Signal Processing to Reconstruct the Signal after Spectral Separating and Detection.
Two cases will be discussed below as nonlimiting illustrations of exemplary embodiments of the present invention, firstly the case of continuous spectral bands (blocks), and secondly the comb-like spectral bands as depicted in
Case A: Continuous Spectral Bands.
The detection arm light is split into N spectral blocks, where each spectral block contains the intensity between two optical frequencies,
The signal for the full spectral width is obtained by an FFT of the signal in each band, an optional compensation of dispersion and other corrections to the phase and amplitude of each Fourier component to optimize the signal and to correct the spectral density for side lobe reduction, addition of the complex FFT spectra, and inverse FFT on the added complex FFT spectrum, optionally with data reduction before the inverse FFT, to obtain the optionally demodulated function R(t), which is the interferometric response for a depth scan with the full source spectrum.
Case B1: Comb like spectral bands and the reconstruction of the full depth range in the sample arm from reduced reference arm scans.
The following description provided below describes the principle of reconstruction of the full depth range in the sample arm from reduced reference arm scans according to the present invention. The procedure shall be explained in the case of separating the spectrum in two spectral bands. The exemplary method can be expanded for separating into many spectral bands.
The signal at the detector for a single detector system is defined by R(t). The depth range in the sample is given by the measurement time T of a single A-line (depth profile) times the group velocity generated by the reference arm delay line, zrange=νgT
The smallest resolvable frequency after an FFT is given by 1/T, which gives a smallest resolvable angular frequency Δω=2π/T. The filter as depicted in
B1(t) and B2(t) are the signals in band one and two respectively. The signal in spectral bands one and two after Fourier transform are given by B1(ω)=R(ω)cos2(ωT/4) and B2(ω)=R(ω)sin2(ωT/4).
This product in the Fourier domain can also be written as a convolution in the time domain. Assuming the signals periodic with time T, the signals B1(t) and B2(t) are given by B1(t)=R(t)+R(t+T/2) and B2(t)=R(t)−R(t+T/2).
Using the above equations, the signal R(t) from t=0 to t=T can be reconstructed from the signals B1(t) and B2(t) recorded from t=0 to t=T/2 by writing,
R(t)=B1(t)+B2(t) and R(t+T/2)=B1(t)−B2(t) for 0<t<T/2. For higher N>2, the identical procedure is performed such that R(t) is reconstructed from B1 to BN.
This demonstrates that the signals B1(t) and B2(t) only need to be recorded over half the depth range zrange. Thus, the depth ranging in the reference arm can be reduced by a factor of 2, while the ranging depth in the sample remains the same. If the signal is split into more spectral bands, like shown in
An exemplary flow diagram of the procedure described above is shown in
Case B2. Limit of large number of spectral bands
In the limit of a large number of spectral bands, N≧
the optical path length change in the reference arm approaches that of a wavelength, λ. In this limit, only a phase change across one wavelength is needed for reconstructing the entire axial scan over length L. In this case, the reference arm path delay may be accomplished by using any of the aforementioned ways for scanning the reference arm delay. Other preferred methods according to the present invention include insertion of an electrooptic modulator, acoustooptic modulator or phase control rapidly scanning optical delay line (“RSOD”) in the reference arm path to impart the path length delay of one wavelength. Also in this case, the wavelength separating unit does not separate the wavelengths into a comb pattern, but separates the spectrum into unique optical frequencies, with each frequency detected by a single detector.
Case C. Fourier domain reconstruction for arbitrary wavelength patterns
In contrast to the reconstruction of the LCI or OCT signal in the time or space domains, the signal may be reconstructed in the Fourier domain by adding the complex spectral components for each wavelength band to compose the Fourier transform of the LCI or OCT signal. Alterations of the phase for each Fourier component may be preferred in certain selected circumstances to correct for minimization of reference arm delay length.
Reconstruction of the Image or One Dimensional Axial Scan
Following reconstruction of the LCI or OCT signal in the real domain, the axial reflectivity may be determined by demodulating the reconstructed LCI or OCT signal. An arrangement for demodulation can include multiplication by a sinusoid and low pass filtering, envelope demodulation using envelope detection, square law demodulation and low pass filtering, quadrature demodulation followed by FIR, IIR filtering, or low pass filtering. In addition, the reconstruction of Stokes vectors (polarization) and flow from these LCI or OCT signals is known to those skilled in the art. Following reconstruction and demodulation, the data may be displayed in one or two-dimensional format (image) for interpretation and ultimately diagnosis of a tissue condition or defect in a medium. If the LCI or OCT signal is reconstructed in the Fourier domain, such reconstructed signal in the Fourier domain can be demodulated in the Fourier domain by shifting the Fourier spectrum and performing an inverse Fourier transform. As a result, the complex signal in the real domain (quadrature signal) is then reconstructed into axial reflectivity information by computing the amplitude of the real portion of the quadrature signal. The complex component is used for computing polarization or flow information. Alternatively, if the signal is reconstructed in the Fourier domain, it can be directly inverse Fourier transformed into the real domain and undergo the aforementioned processing described for the reconstructed real domain signals.
After spectral compounding of the source light in the first 50/50 splitter and the 80/20 splitter, light enters a modified Michelson interferometer. A configuration that implements balanced detection is shown. The sample arm goes to the probe (e.g., a slit lamp). Reference arm light is transmitted through two acousto-optic modulators with a difference frequency of 10 kHz to generate a constant carrier frequency that is independent of wavelength. The balanced detection outputs go to separate spectral detection units.
Spectral Detection Unit
Referring to
A scan cameras with N detector elements is used as spectral detection unit 128 (see
Line scan rates of 20 kHz can be achieved, allowing demodulation of a 10 kHz carrier to extract the complex cross-spectral density. Data is digitized and transferred to computer memory. Demodulation of the signal is done in software. Scan rates of 10,000 depth profiles per second or more can be achieved.
Dual Balanced Detection
Dual balanced detection is preferably used by the present invention, which is preferably utilized for the following reasons. Firstly, most light sources generate 1/f noise (f=frequency) at relatively low frequencies. Balanced detection will eliminate 1/f source noise. Secondly, an interference term of the sample arm light with itself (auto-correlation term) is present on top of the true signal term, which is the interference between sample and reference arm. This auto-correlation term can be eliminated by a differential technique. Balanced detection may eliminate this auto-correlation term from the measured signal. Thirdly, RIN can be reduced.
Data Acquisition and Processing Unit
The data rate at 20,000 spectral profiles per second, with 2000 detector elements and 8-10 bit resolution (the dynamic range of most line scan cameras) is 40-80 MB/sec. Maximum sustainable data transfer speed over the PCI bus is 100 MB/sec. In a computer with two independent PCI bridges to computer system memory, approximately 200 MB/sec of data can be transferred for real time processing of data from two line scan cameras simultaneously. Implementation of dual balanced detection in analog by subtracting line scan camera signals before digitization may reduce the data rate by a factor of 2. High-speed data acquisition boards are available at resolutions of 12-14 bits and speeds up to 100 Msamples/sec. A single 2048 point fast Fourier transform on a 2.5 GHz Pentium 4 processor takes 50 μsec. These numbers show that real-time processing of Spectral Domain OCT data at 20,000 spectral profiles/sec is within reach of current data acquisition and processing power of dual processor PC's . The data collected by the spectrometer can be sampled with equal wavelength increments. Fourier transform however, links z and k space (or t and w). Because of the nonlinear relation between k and λ the spectrum from the spectrometer should be interpolated to create evenly spaced samples in k domain. To achieve the optimal point spread function, dispersion in the sample and reference arm of the interferometer should be balanced. We have shown that dispersion imbalance can be corrected for by digital processing, allowing for correct compensation of dispersion for individual eye lengths.
Phase Tracking
The present invention also provides apparatus and methods for phase tracking in spectral domain (“SD”) OCT.
Fully Parallel SD OCT
One of the features of fully parallel SD OCT is spectral dispersion of the detection arm light onto a multi-element array such as but not limiting to an integrating device (e.g., CCD) and measurement of the real or complex spectral density at high speeds. The detection arm beam is separated by a spectral separating unit (e.g., grating) and focused onto the array. With respect to previous Spectral Domain OCT designs known in the art, two differences are apparent that will be discussed below: 1) implementation of balanced detection, and, 2) implementation of phase tracking.
Spectrometer design The depth range in SD OCT is inversely proportional to the spectral resolution. Using the complex spectral density, ranging depth z is given by,
Dual balanced detection: Dual balanced detection is advantageous for at least three reasons. First, most light sources generate 1/f noise at relatively low frequencies (tens of kHz range). In time domain (“TD”) OCT systems 1/f noise is not a problem because the signal carrier is in general in the MHz range where 1/f noise is not significant. In SD OCT, balanced detection may likely eliminate 1/f source noise. Second, an interference of the sample arm light with itself (auto-correlation term) is present on top of the true signal. This auto-correlation term can be eliminated by a differential technique. Balanced detection can be used to eliminate this auto-correlation term from the measured signal. Third, balanced detection may reduce relative intensity or Bose Einstein noise.
Phase Tracking: Phase tracking is preferable to eliminate phase instabilities in the interferometer. Phase instabilities can cause individual interferometric fringes to shift in location. If detection is slow relative to the shifting of the fringes, the resulting averaging results in an artifactual decrease in the measured fringe amplitude. Fast detection arrays can capture the cross spectral density at a rate of 20 to 40 kHz, resulting in integration times of 50 to 25 μsec, respectively. Phase instabilities arising on a time f lame shorter than the integration time of the array should be compensated.
Phase locking circuitry is common in electronics, and is frequently used in radar and ultrasound. Active phase tracking can be implemented by modulating the interferometer path length difference at 10 MHz with an electro-optic phase modulator in the reference arm over a fraction of the wavelength. By demodulating the intensity measured by one detector at the output of the interferometer at the frequency of the path length modulation, an error signal can be generated indicating in which direction the phase modulator should shift to lock onto a fringe amplitude maximum. By adding an offset to the phase modulator as determined by the error signal, the phase tracker actively locks onto a fringe maximum. The phase modulator can only modulate the path length difference over a few wavelengths. The processing unit can determine if the phase modulator has reached its range limit, and jump by a full wave in phase to maintain lock on a different fringe maximum. This approach exploits the fact that phase should be controlled only modulo 2π. In addition, the processing drives a slower component (e.g., the Rapid Scanning Optical Delay line) to extend the path length range of the phase modulator/RSOD combination over several millimeters. Phase locking can be performed on a fringe maximum, minimum, or zero crossing, based on the type of mixing performed in the demodulation circuit.
The present invention can also use autoranging technology, including processing algorithms, as disclosed in copending U.S. application Ser. No. 10/136,813, filed Apr. 30, 2002, entitled METHOD AND APPARATUS FOR IMPROVING IMAGE CLARITY AND SENSITIVITY IN OPTICAL COHERENCE TOMOGRAPHY USING DYNAMIC FEEDBACK TO CONTROL FOCAL PROPERTIES AND COHERENCE GATING, and commonly assigned to the assignee of the present invention, the disclosure of which is incorporated herein.
The autoranging mechanism may, in one exemplary embodiment, comprise a processor unit for (a) obtaining a first scan line; (b) locating a surface location “S” of a sample; (c) locating an optimal scan range “R” of the sample;(d) modifying a reference arm delay waveform to provide an output; (e) outputting the output to a reference arm; (f) determining whether the image is complete; and (g) moving to the next scan line if the image is not complete or remapping the image using the surface S data and the waveform data stored in the memory storage device if the image is complete.
If the light returned from the sample is of low amplitude, phase locking may be unstable due to the presence of noise. In another embodiment, a separate, preferably monochromatic, light source is input into the interferometer. The separate source wavelength may overlap with the broad bandwidth OCT or LCI source spectrum or may be centered at a different wavelength than the OCT or LCI source spectrum. The separate source is preferably of higher power and may be combined with the source arm (using wavelength division multiplexer, grating, prism, filter or the like) travel to the reference and sample arms and return back to the beam recombining element. The returned separate source light can then separated from the OCT or LCI light following transmission back through the beam recombining element (i.e. beam splitter output). A separation arrangement can perform spectral separation by a dispersing element, such as a dichroic mirror, filter, grating, prism, wavelength division multiplexer or the like. The separate source will be detected separately from the OCT or LCI broad bandwidth light using one or more detectors. The higher power provided by this separate source can enable detection of a higher amplitude interference pattern, and provide an improved input to the phase tracker, thus enabling more stable phase tracking.
Mixer Implementation: The intensity I(t) at the detector at a given moment within a single oscillation of the fringe pattern is given by
I(t)=cos [φ(t)]
where the phase φ gives the position in the fringe. For φ=0, the signal is at a fringe maximum, for φ=π, the signal is at a fringe minimum. At an arbitrary moment t, the phase φ(t) is given by,
φ(t)=α+β sin(ωt)
where α describes the position within a single oscillation of the fringe pattern, and β*sin(ωt) is the phase modulation introduced by the phase modulator, with β the amplitude of the phase modulation, and ω the frequency of the phase modulation signal. The intensity at the photodetector I(t) can be mixed with a carrier at frequency ω and 2ω, resulting in the mixer signal MixerC(t), MixerS(t), Mixer2ωC(t) and Mixer2ωS(t),
MixerC(t)=cos(ωt)*cos(α+β sin(ωt)); MixerS(t)=sin(ωt)*cos(α+β sin(ωt));
Mixer2ωC(t)=cos(2ωt)*cos(α+β sin(ωt));
Mixer2ωS(t)=sin(2ωt)*cos(α+β sin(ωωt))
The time average over a single oscillation of the carrier frequency ω of MixerC, MixerS, Mixer2ωC and Mixer2ωS is given by,
where J1(β) and J2(β) are a Bessel functions of the first kind; its value depends on β, the amplitude of the phase modulation. Thus, the signal
The present invention provides a method for tracking phase in an imaging system, as shown in
The method further may comprise that steps (a)-(f) are performed in parallel with other imaging processes. The adjustment of phase “φ” is defined as A(x2-x1), where “A” is a constant. Furthermore, optionally, step d) may further comprise the substeps of d1) determining whether A(x2-x1) is within range of the phase modulator; and d2) changing φ by an amount equal to A(x2-x1) if A(x2-x1) is within the range or changing φ by an amount equal to A(x2-x1)—m2π if A(x2-x1) is outside of the range, where M is an integer greater than 1. The method may optionally further comprise a substep d3) remeasuring signal x1.
Data Acquisition and Processing Unit
In general, the data collected by the spectrometer are sampled with equal wavelength increments. Fourier transform, however, links z and k space (or t and w). Because of the non-linear relation between k and λ the acquired spectrum is interpolated to create evenly spaced samples in the k domain. Alternatively, the light could be dispersed in such a way on the detection array that the light is samples in equal intervals in k space, such that the interpolation becomes obsolete. Alternatively, the detection array spacing could be designed to sample the light evenly spread in the k domain, such that the interpolation becomes obsolete. To achieve the optimal point spread function, dispersion in the sample and reference arm of the interferometer should preferably be balanced. Dispersion imbalance can be corrected by digital processing.
The present invention provides a probe for locating atherosclerotic plaque in a blood vessel, comprising: an interferometer; a spectral separating unit which splits signal received from the interferometer into a plurality of optical frequencies; and a detector arrangement capable of detecting at least a portion of the optical frequencies received from the spectral separating unit.
The present invention further provides an apparatus for delivering a therapeutic agent, comprising: a probe disposed in the housing and comprising: an interferometer, a spectral separating unit which splits signal received from the interferometer into a plurality of optical frequencies, a detector arrangement capable of detecting at least a portion of the optical frequencies received from the spectral separating unit; and a conduit cooperating with the probe, and comprising a proximal end for receiving the therapeutic agent and a distal end for delivering the therapeutic agent at a predetermined location, the location being determined by imaging the environment in proximity to the distal end using the probe.
An exemplary embodiment of the present invention will be further described below in connection with the following example, which is set forth for purposes of illustration only.
The method according to the present invention was verified in the laboratory by the following experiment.
In the existing OCT system, the shot noise power spectrum as determined from the spectral density due to the reference arm optical power was measured. Then ⅔ of the spectrum from the reference arm was blocked, and experimentally it was verified that the shot noise power spectrum was reduced by a factor of three, thus demonstrating that the shot noise is reduced by a factor of 3 if the spectrum is split in three spectral bands (see
An object with low reflectivity was inserted in the sample arm. Using the full spectral width of the source, the power spectrum of the interference between sample and reference arm light was determined in the lower half of the spectral density. Then the upper part of the source spectrum was blocked in the reference arm, and it was verified that the lower ⅓ of the power spectrum of the interference between sample and reference arm light had the same magnitude as in the previous measurement (see
This demonstrates that when the light in the detection arm is split in two spectral bands, the spectral density of the interference between sample and reference arm light within the spectral bandwidth of a single detector is unchanged. Combined with the measurement that showed a reduction in the shot noise power spectrum, the conclusion is that a reduction of shot noise can be realized by splitting the detection arm light in separate spectral bands.
Experimental Verification of the Noise Reduction.
To demonstrate the noise reduction in Spectral Domain OCT, an OCT system was used, including a Rapid Scanning Optical Delay line (RSOD) was used in the reference arm, enabling portions of the spectrum to be blocked. Detector signals were digitized at 2.5 Msamples/sec, allowing digital processing of the fringe information. First, the thermal noise density of the detector was measured as a function of frequency by blocking all light onto the detector. Second, the shot noise density of the reference arm power was measured with only the reference arm power incident on the detector. Third, both the sample and reference arm light were incident on the detector. The sample was a single scattering surface mounted in a model eye and 512 depth profiles were acquired in 2 seconds. The power density I(f)2 was measured, which is proportional to the spectral density squared (see Eq. (9)). Then we blocked half of the spectrum in the reference and measured again the shot noise density of the reference arm by blocking the sample arm, and the power density I(f)2 when both sample and reference arm light were incident on the detector. Shot noise and power densities were corrected for thermal noise by subtraction. Thermal noise was at least a factor of 3 smaller than the lowest shot noise level.
As is evident from
The next experiment further demonstrated that by dispersing the spectrum in the detection arm over several detectors, and by selectively band pass filtering the signals of each detector, the SNR is increased. The detection arm light was dispersed over 4 detectors by a diffraction grating as shown in
First, the thermal noise density of all four detectors was measured. Second, the shot noise density of the reference arm light in each detector channel 116 was measured. Third, both the sample and reference arm light were incident on the detector 114. The sample 130 was a single scattering surface mounted in a model eye and 512 depth profiles were acquired in 2 seconds. The power density I(f)2 in each detector channel 114 was measured. Then, the signals of the four detectors 114 were summed, and the combined power density I(f)2 was determined. The results are shown in
In
These experiments clearly demonstrate that spectrally dispersing the light in the detection arm can offer a significant SNR advantage.
Although only a few exemplary embodiments of this invention have been described in detail above, those skilled in the art will readily appreciate that many modifications are possible in the exemplary embodiments without materially departing from the novel teachings and advantages of this invention. Accordingly, all such modifications are intended to be included within the scope of this invention as defined in the following claims. It should further be noted that any patents, applications and publications referred to herein are incorporated by reference in their entirety.
This application claims benefit of copending U.S. provisional patent application No. 60/351,904, filed Jan. 24, 2002, entitled APPARATUS AND METHOD FOR RANGING AND SHOT NOISE REDUCTION OF LOW COHERENCE INTERFEROMETRY (LCI) AND OPTICAL COHERENCE TOMOGRAPHY (OCT) SIGNALS BY PARALLEL DETECTION OF SPECTRAL BANDS, and copending U.S. application Ser. No. 10/136,813, filed Apr. 30, 2002, entitled METHOD AND APPARATUS FOR IMPROVING IMAGE CLARITY AND SENSITIVITY IN OPTICAL COHERENCE TOMOGRAPHY USING DYNAMIC FEEDBACK TO CONTROL FOCAL PROPERTIES AND COHERENCE GATING, both commonly assigned to the assignee of the present application, the disclosures of which are incorporated by reference in its entirety herein.
The invention was made with the U.S. Government support under Grant Number DAMD17-02-2-0006 awarded by the U.S. Department of the Army. Thus, the U.S. Government has certain rights in the invention.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/US03/02349 | 1/24/2003 | WO | 00 | 7/9/2004 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO03/062802 | 7/31/2003 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
2339754 | Brace | Jan 1944 | A |
3601480 | Randall | Aug 1971 | A |
3856000 | Chikama | Dec 1974 | A |
3941121 | Olinger | Mar 1976 | A |
4030827 | Delhaye et al. | Jun 1977 | A |
4141362 | Wurster | Feb 1979 | A |
4295738 | Meltz et al. | Oct 1981 | A |
4300816 | Snitzer et al. | Nov 1981 | A |
4479499 | Alfano | Oct 1984 | A |
4533247 | Epworth | Aug 1985 | A |
4585349 | Gross et al. | Apr 1986 | A |
4601036 | Faxvog et al. | Jul 1986 | A |
4607622 | Fritch et al. | Aug 1986 | A |
4631498 | Cutler | Dec 1986 | A |
4770492 | Levin et al. | Sep 1988 | A |
4868834 | Fox et al. | Sep 1989 | A |
4925302 | Cutler | May 1990 | A |
4928005 | Lefevre et al. | May 1990 | A |
4965441 | Picard | Oct 1990 | A |
4993834 | Carlhoff et al. | Feb 1991 | A |
5039193 | Snow et al. | Aug 1991 | A |
5040889 | Keane | Aug 1991 | A |
5045936 | Lobb et al. | Sep 1991 | A |
5046501 | Crilly | Sep 1991 | A |
5065331 | Vachon et al. | Nov 1991 | A |
5120953 | Harris | Jun 1992 | A |
5127730 | Brelje et al. | Jul 1992 | A |
5197470 | Helfer et al. | Mar 1993 | A |
5202745 | Sorin et al. | Apr 1993 | A |
5248876 | Kerstens et al. | Sep 1993 | A |
5262644 | Maguire | Nov 1993 | A |
5291885 | Taniji et al. | Mar 1994 | A |
5293872 | Alfano et al. | Mar 1994 | A |
5293873 | Fang | Mar 1994 | A |
5304810 | Amos | Apr 1994 | A |
5305759 | Kaneko et al. | Apr 1994 | A |
5317389 | Hochberg et al. | May 1994 | A |
5321501 | Swanson et al. | Jun 1994 | A |
5353790 | Jacques et al. | Oct 1994 | A |
5383467 | Auer et al. | Jan 1995 | A |
5411016 | Kume et al. | May 1995 | A |
5419323 | Kittrell et al. | May 1995 | A |
5439000 | Gunderson et al. | Aug 1995 | A |
5441053 | Lodder et al. | Aug 1995 | A |
5450203 | Penkethman | Sep 1995 | A |
5454807 | Lennox et al. | Oct 1995 | A |
5459325 | Hueton et al. | Oct 1995 | A |
5459570 | Swanson et al. | Oct 1995 | A |
5465147 | Swanson | Nov 1995 | A |
5486701 | Norton et al. | Jan 1996 | A |
5491524 | Hellmuth et al. | Feb 1996 | A |
5491552 | Knuttel | Feb 1996 | A |
5526338 | Hasman et al. | Jun 1996 | A |
5562100 | Kittrell et al. | Oct 1996 | A |
5565986 | Knüttel | Oct 1996 | A |
5583342 | Ichie | Dec 1996 | A |
5590660 | MacAulay et al. | Jan 1997 | A |
5600486 | Gal et al. | Feb 1997 | A |
5601087 | Gunderson et al. | Feb 1997 | A |
5621830 | Lucey et al. | Apr 1997 | A |
5623336 | Raab et al. | Apr 1997 | A |
5697373 | Richards-Kortum et al. | Dec 1997 | A |
5698397 | Zarling et al. | Dec 1997 | A |
5710630 | Essenpreis et al. | Jan 1998 | A |
5719399 | Alfano et al. | Feb 1998 | A |
5735276 | Lemelson | Apr 1998 | A |
5740808 | Panescu et al. | Apr 1998 | A |
5748598 | Swanson et al. | May 1998 | A |
5784352 | Swanson et al. | Jul 1998 | A |
5785651 | Kuhn et al. | Jul 1998 | A |
5795295 | Hellmuth et al. | Aug 1998 | A |
5801826 | Williams | Sep 1998 | A |
5803082 | Stapleton et al. | Sep 1998 | A |
5807261 | Benaron et al. | Sep 1998 | A |
5817144 | Gregory | Oct 1998 | A |
5840023 | Oraevsky et al. | Nov 1998 | A |
5842995 | Mahadevan-Jansen et al. | Dec 1998 | A |
5843000 | Nishioka et al. | Dec 1998 | A |
5843052 | Benja-Athon | Dec 1998 | A |
5847827 | Fercher | Dec 1998 | A |
5862273 | Pelletier | Jan 1999 | A |
5865754 | Sevick-Muraca et al. | Feb 1999 | A |
5867268 | Gelikonov et al. | Feb 1999 | A |
5871449 | Brown | Feb 1999 | A |
5872879 | Hamm | Feb 1999 | A |
5877856 | Fercher | Mar 1999 | A |
5887009 | Mandella et al. | Mar 1999 | A |
5892583 | Li | Apr 1999 | A |
5920373 | Bille | Jul 1999 | A |
5920390 | Farahi et al. | Jul 1999 | A |
5921926 | Rolland et al. | Jul 1999 | A |
5949929 | Hamm | Sep 1999 | A |
5951482 | Winston et al. | Sep 1999 | A |
5956355 | Swanson et al. | Sep 1999 | A |
5968064 | Selmon et al. | Oct 1999 | A |
5983125 | Alfano et al. | Nov 1999 | A |
5987346 | Benaron et al. | Nov 1999 | A |
5991697 | Nelson et al. | Nov 1999 | A |
5994690 | Kulkarni et al. | Nov 1999 | A |
6002480 | Izatt et al. | Dec 1999 | A |
6004314 | Wei et al. | Dec 1999 | A |
6006128 | Izatt et al. | Dec 1999 | A |
6010449 | Selmon et al. | Jan 2000 | A |
6014214 | Li | Jan 2000 | A |
6033721 | Nassuphis | Mar 2000 | A |
6044288 | Wake et al. | Mar 2000 | A |
6048742 | Weyburne et al. | Apr 2000 | A |
6053613 | Wei et al. | Apr 2000 | A |
6069698 | Ozawa et al. | May 2000 | A |
6091496 | Hill | Jul 2000 | A |
6091984 | Perelman et al. | Jul 2000 | A |
6111645 | Tearney et al. | Aug 2000 | A |
6117128 | Gregory | Sep 2000 | A |
6120516 | Selmon et al. | Sep 2000 | A |
6134003 | Tearney et al. | Oct 2000 | A |
6134010 | Zavislan | Oct 2000 | A |
6134033 | Bergano et al. | Oct 2000 | A |
6141577 | Rolland et al. | Oct 2000 | A |
6151522 | Alfano et al. | Nov 2000 | A |
6159445 | Klaveness et al. | Dec 2000 | A |
6160826 | Swanson et al. | Dec 2000 | A |
6161031 | Hochmann et al. | Dec 2000 | A |
6166373 | Mao | Dec 2000 | A |
6174291 | McMahon et al. | Jan 2001 | B1 |
6175669 | Colston et al. | Jan 2001 | B1 |
6185271 | Kinsinger | Feb 2001 | B1 |
6191862 | Swanson et al. | Feb 2001 | B1 |
6193676 | Winston et al. | Feb 2001 | B1 |
6198956 | Dunne | Mar 2001 | B1 |
6201989 | Whitehead et al. | Mar 2001 | B1 |
6208415 | De Boer et al. | Mar 2001 | B1 |
6208887 | Clarke | Mar 2001 | B1 |
6263234 | Engelhardt et al. | Jul 2001 | B1 |
6264610 | Zhu | Jul 2001 | B1 |
6272376 | Marcu et al. | Aug 2001 | B1 |
6274871 | Dukor et al. | Aug 2001 | B1 |
6282011 | Tearney et al. | Aug 2001 | B1 |
6308092 | Hoyns | Oct 2001 | B1 |
6324419 | Guzelsu et al. | Nov 2001 | B1 |
6341036 | Tearney et al. | Jan 2002 | B1 |
6353693 | Kano et al. | Mar 2002 | B1 |
6377349 | Fercher | Apr 2002 | B1 |
6384915 | Everett et al. | May 2002 | B1 |
6393312 | Hoyns | May 2002 | B1 |
6394964 | Sievert, Jr. et al. | May 2002 | B1 |
6421164 | Tearney et al. | Jul 2002 | B2 |
6445944 | Ostrovsky | Sep 2002 | B1 |
6463313 | Winston et al. | Oct 2002 | B1 |
6469846 | Ebizuka et al. | Oct 2002 | B2 |
6485413 | Boppart et al. | Nov 2002 | B1 |
6485482 | Belef | Nov 2002 | B1 |
6501551 | Tearney et al. | Dec 2002 | B1 |
6501878 | Hughes et al. | Dec 2002 | B2 |
6549801 | Chen et al. | Apr 2003 | B1 |
6552796 | Magnin et al. | Apr 2003 | B2 |
6556305 | Aziz et al. | Apr 2003 | B1 |
6556853 | Cabib et al. | Apr 2003 | B1 |
6558324 | Von Behren et al. | May 2003 | B1 |
6564087 | Pitris et al. | May 2003 | B1 |
6564089 | Izatt et al. | May 2003 | B2 |
6615071 | Casscells, III et al. | Sep 2003 | B1 |
6622732 | Constantz | Sep 2003 | B2 |
6680780 | Fee | Jan 2004 | B1 |
6685885 | Nolte et al. | Feb 2004 | B2 |
6687007 | Meigs | Feb 2004 | B1 |
6687010 | Horii et al. | Feb 2004 | B1 |
6741355 | Drabarek | May 2004 | B2 |
6790175 | Furusawa et al. | Sep 2004 | B1 |
6806963 | Wälti et al. | Oct 2004 | B1 |
6816743 | Moreno et al. | Nov 2004 | B2 |
6980299 | de Boer | Dec 2005 | B1 |
7006231 | Ostrovsky et al. | Feb 2006 | B2 |
7231243 | Tearney et al. | Jun 2007 | B2 |
20010047137 | Moreno et al. | Nov 2001 | A1 |
20020016533 | Marchitto et al. | Feb 2002 | A1 |
20020064341 | Fauver et al. | May 2002 | A1 |
20020076152 | Hughes et al. | Jun 2002 | A1 |
20020122246 | Tearney et al. | Sep 2002 | A1 |
20020161357 | Rox et al. | Oct 2002 | A1 |
20020163622 | Magnin et al. | Nov 2002 | A1 |
20020172485 | Keaton et al. | Nov 2002 | A1 |
20020188204 | McNamara et al. | Dec 2002 | A1 |
20020196446 | Roth et al. | Dec 2002 | A1 |
20020198457 | Tearney et al. | Dec 2002 | A1 |
20030023153 | Izatt et al. | Jan 2003 | A1 |
20030026735 | Nolte et al. | Feb 2003 | A1 |
20030135101 | Webler | Jul 2003 | A1 |
20030164952 | Deichmann et al. | Sep 2003 | A1 |
20030171691 | Casscells, III et al. | Sep 2003 | A1 |
20030236443 | Cespedes et al. | Dec 2003 | A1 |
20040086245 | Farroni et al. | May 2004 | A1 |
20040100631 | Bashkansky et al. | May 2004 | A1 |
20040100681 | Bjarklev et al. | May 2004 | A1 |
20040150829 | Koch et al. | Aug 2004 | A1 |
20040166593 | Nolte et al. | Aug 2004 | A1 |
20050018201 | De Boer | Jan 2005 | A1 |
20050075547 | Wang | Apr 2005 | A1 |
20050083534 | Riza et al. | Apr 2005 | A1 |
Number | Date | Country |
---|---|---|
4309056 | Sep 1994 | DE |
19542955 | May 1997 | DE |
10351319 | Jun 2005 | DE |
0110201 | Jun 1984 | EP |
0251062 | Jan 1988 | EP |
0590268 | Apr 1994 | EP |
1426799 | Jun 2004 | EP |
1257778 | Dec 1971 | GB |
2030313 | Apr 1980 | GB |
2209221 | May 1989 | GB |
4135550 | May 1992 | JP |
4135551 | May 1992 | JP |
9219930 | Nov 1992 | WO |
9303672 | Mar 1993 | WO |
9533971 | Dec 1995 | WO |
9732182 | Sep 1997 | WO |
9801074 | Jan 1998 | WO |
9814132 | Apr 1998 | WO |
9835203 | Aug 1998 | WO |
9846123 | Oct 1998 | WO |
9838907 | Nov 1998 | WO |
9848838 | Nov 1998 | WO |
9944089 | Sep 1999 | WO |
9957507 | Nov 1999 | WO |
0058766 | Oct 2000 | WO |
0138820 | May 2001 | WO |
0142735 | Jun 2001 | WO |
0236015 | May 2002 | WO |
0238040 | May 2002 | WO |
02054027 | Jul 2002 | WO |
03020119 | Mar 2003 | WO |
03062802 | Jul 2003 | WO |
20040066824 | Aug 2004 | WO |
2004088361 | Oct 2004 | WO |
2004105598 | Dec 2004 | WO |
20050000115 | Jan 2005 | WO |
2005054780 | Jun 2005 | WO |
2006014392 | Feb 2006 | WO |
2006130797 | Dec 2006 | WO |
Number | Date | Country | |
---|---|---|---|
20050018201 A1 | Jan 2005 | US |
Number | Date | Country | |
---|---|---|---|
60351904 | Jan 2002 | US |