Phase-resolved optical coherence tomography and optical doppler tomography for imaging fluid flow in tissue with fast scanning speed and high velocity sensitivity
First Claim
1. A method for phase-resolved optical tomography capable of simultaneously imaging fluid flow and morphology in a sample with fast scanning speed and high velocity sensitivity comprising:
- providing a source of at least partially coherent radiation through an interferometer, said at least partially coherent radiation characterized by a phase;
phase modulating said radiation in said interferometer at a modulation frequency;
scanning said sample with said source of at least partially coherent radiation through said interferometer in a sequence of pixel line scans, said sample having a fluid flow therein;
changing said phase of said at least partially coherent radiation in response to said fluid flow at each pixel of each pixel line scan;
detecting interference fringes of said radiation backscattered from said sample into said interferometer;
processing said detected interference fringes to determine ODT phase signals of said detected backscattered interference fringes at each pixel of said pixel line scans to determine said corresponding phase at each pixel in a data window;
comparing said phase between corresponding pixels in two line scans to generate a difference between said phase at said two corresponding pixels in two line scans whereby speckle is substantially reduced; and
generating a tomographic image of the fluid flow in said sample from said phase difference at each pixel.
2 Assignments
0 Petitions
Accused Products
Abstract
The invention is a fast-scanning ODT system that uses phase information derived from a Hilbert transformation to increase the sensitivity of flow velocity measurements while maintaining high spatial resolution. The significant increases in scanning speed and velocity sensitivity realized by the invention make it possible to image in vivo blood flow in human skin. The method of the invention overcomes the inherent limitations of the prior art ODT by using a phase change between sequential line scans for velocity image reconstruction. The ODT signal phase or phase shifts at each pixel can be determined from the complex function, {tilde over (Γ)}ODT(t), which is determined through analytic continuation of the measured interference fringes function, ΓODT(t), by use of a Hilbert transformation, by electronic phase demodulation, by optical means, or a fast Fourier transformation. The phase change in each pixel between axial-line scans is then used to calculate the Doppler frequency shift. Sequential measurements of a single line scan, measurements of sequential line scans or measurements of line scans in sequential frames may be used. Because the time between line scans is much longer than the pixel time window, very small Doppler shifts can be detected with this technique. In addition, spatial resolution and velocity sensitivity are decoupled. Furthermore, because two axial-line scans are compared at the same location, speckle modulations in the fringe signal cancel each other and, therefore, will not affect the phase-difference calculation.
-
Citations
42 Claims
-
1. A method for phase-resolved optical tomography capable of simultaneously imaging fluid flow and morphology in a sample with fast scanning speed and high velocity sensitivity comprising:
-
providing a source of at least partially coherent radiation through an interferometer, said at least partially coherent radiation characterized by a phase;
phase modulating said radiation in said interferometer at a modulation frequency;
scanning said sample with said source of at least partially coherent radiation through said interferometer in a sequence of pixel line scans, said sample having a fluid flow therein;
changing said phase of said at least partially coherent radiation in response to said fluid flow at each pixel of each pixel line scan;
detecting interference fringes of said radiation backscattered from said sample into said interferometer;
processing said detected interference fringes to determine ODT phase signals of said detected backscattered interference fringes at each pixel of said pixel line scans to determine said corresponding phase at each pixel in a data window;
comparing said phase between corresponding pixels in two line scans to generate a difference between said phase at said two corresponding pixels in two line scans whereby speckle is substantially reduced; and
generating a tomographic image of the fluid flow in said sample from said phase difference at each pixel. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16)
where P denotes the Cauchy principle value and A(t) and ψ
(t) are the amplitude and the phase of {tilde over (Γ
)}ODT(t) respectively.
-
-
6. The method of claim 5 where processing said detected interference fringes to determine ODT phase signals of said detected backscattered interference fringes through the use of a Hilbert transformation comprises determining said phase shift, Δ
- ψ
, in said computer at each pixel between sequential line scans by the Doppler frequency shift at each pixel between, ω
=Δ
ω
/T where T is the time interval between successive line scans.
- ψ
-
7. The method of claim 1 where processing said detected interference fringes to determine ODT phase signals of said detected backscattered interference fringes comprises determining said ODT phase signals in a computer or digital signal processor by:
-
applying a first fast Fourier transformation in said computer to said detected interference fringes to produce a transformed signal;
passing said transformed signal through a digital filter in said computer to reduce noise to produce a filtered signal;
multiplying said filtered signal by the Heaviside function, H(ω
) in said computer to produce a multiplied signal; and
applying a second fast Fourier transformation in said computer to said multiplied signal to producer {tilde over (Γ
)}ODT(t) ; and
extracting in said computer said ODT phase signals from {tilde over (Γ
)}ODT(t).
-
-
8. The method of claim 1 wherein said interferometer has a reference arm, a sample arm and a detection arm and further comprising prior to detecting interference fringes of said radiation backscattered from said sample into said interferometer:
-
providing circularly polarized reflected radiation from said reference arm;
providing linearly polarized reflected radiation from said sample arm;
combining said circularly polarized reflected radiation from said reference arm and said linearly polarized reflected radiation from said sample arm in a coupler; and
splitting radiation from said coupler into two polarized beams in said detection arm with a polarizing beam splitter with an axis of polarization set at 45°
to said linearly polarized reflected radiation from said sample arm,wherein detecting interference fringes comprises detecting said two polarized beams in two corresponding detectors to generate two corresponding signals comprising in turn the real and imaginary parts of {tilde over (Γ
)}ODT(t) so that said phase is directly measured.
-
-
9. The method of claim 1 where generating said tomographic image of the fluid flow in said sample comprises generating a structural image of said sample.
-
10. The method of claim 1 where generating said tomographic image of the fluid flow in said sample comprises generating a velocity image of fluid in said sample.
-
11. The method of claim 1 further comprising determining a phase shift Δ
- f=Δ
ψ
/T from an average Doppler frequency shift, <
Δ
f>
, between sequential line scans using the equation;
where T is the time interval between sequential scans, and n is the number of sequential scans averaged.
- f=Δ
-
12. The method of claim 11 where determining a phase shift from an average Doppler frequency shift between sequential line scans comprises averaging at least three sequential line scans to increase the signal-to-noise ratio in images.
-
13. The method of claim 1 further comprising determining a standard deviation of a Doppler power spectrum and generating said tomographic image of the fluid flow in said sample from said standard deviation of said Doppler power spectrum.
-
14. The method of claim 13 where determining said standard deviation of a Doppler power spectrum comprises determining said standard deviation from the equation:
-
where T is the time interval between sequential scans, P(ω
) is said Doppler power spectrum, and {overscore (ω
)} is the centroid value of a Doppler frequency shift.
-
-
15. The method of claim 1 where generating said corresponding phase at each pixel in a data window comprises extracting said phase shift by electronic demodulation of said detected interference fringes.
-
16. The method of claim 1 where said scanning is characterized by a scanning rate, and adding π
- to said phase difference because tan(φ
)=tan(φ
+π
), when a Doppler frequency difference related to said phase difference is greater than half of said scanning rate so that said Doppler frequency difference can exceed said scanning rate.
- to said phase difference because tan(φ
-
17. A method for phase-resolved optical tomography comprising:
-
providing a source of at least partially coherent radiation through an interferometer, said at least partially coherent radiation characterized by a phase;
phase modulating said radiation in said interferometer at a modulation frequency;
scanning said sample with said source of at least partially coherent radiation through said interferometer in a sequence of pixel line scans, said sample having a fluid flow therein;
changing said phase of said at least partially coherent radiation in response to said fluid flow at each pixel of each pixel line scan;
detecting interference fringes of said radiation backscattered from said sample into said interferometer;
processing said detected interference fringes to determine ODT phase signals of said detected backscattered interference fringes at each pixel of said pixel line scans to generate a corresponding rate of change of said phase at each pixel in a data window through the use of a Hilbert transformation at each pixel, and generating a tomographic image of the fluid flow in said sample from said rate of change of said phase at each pixel.
-
-
18. A method for phase-resolved optical tomography comprising:
-
providing a source of at least partially coherent radiation through an interferometer with a reference arm, a sample arm and a detection arm, said at least partially coherent radiation characterized by a phase;
phase modulating said radiation in said interferometer at a modulation frequency;
scanning said sample with said source of at least partially coherent radiation through said interferometer in a sequence of pixel line scans, said sample having a fluid flow therein so that said phase of said at least partially coherent radiation is changed in response to said fluid flow at each pixel of each pixel line scan;
providing circularly polarized reflected radiation from said reference arm;
providing linearly polarized reflected radiation from said sample arm;
combining said circularly polarized reflected radiation from said reference arm and said linearly polarized reflected radiation from said sample arm in a coupler; and
splitting radiation from said coupler into two polarized beams in said detection arm with a polarizing beam splitter with an axis of polarization set at 45°
to said linearly polarized reflected radiation from said sample arm,detecting interference fringes of said radiation backscattered from said sample into said interferometer in said two polarized beams in two corresponding detectors to generate two corresponding signals comprising in turn the real and imaginary parts of {tilde over (Γ
)}ODT(t) so that said phase is directly measured; and
generating a tomographic image of the fluid flow in said sample from said phase at each pixel.
-
-
19. An apparatus for phase-resolved optical tomography capable of simultaneously imaging fluid flow and morphology in sample with fast scanning speed and high velocity sensitivity comprising:
-
an interferometer;
a source of at least partially coherent radiation through coupled to said interferometer, said at least partially coherent radiation characterized by a phase;
phase modulator coupled to said source to modulate said radiation in said interferometer at a modulation frequency;
a scanner for scanning said sample with said source of at least partially coherent radiation through said interferometer in a sequence of pixel line scans, said sample having a fluid flow therein so that said phase of said at least partially coherent radiation is changed in response to said fluid flow at each pixel of each pixel line scan;
a detector to detect interference fringes of said radiation backscattered from said sample into said interferometer;
a processor to determine said corresponding phase at each pixel of said pixel line scans from said ODT phase signals of said detected backscattered interference fringes and to compare said phase between corresponding pixels in two line scans to generate a difference between said phase at said two corresponding pixels in two line scans whereby speckle is substantially reduced; and
a display to generate a tomographic image of the fluid flow in said sample from said difference at each pixel. - View Dependent Claims (20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32)
where P denotes the Cauchy principle value and A(t) and ψ
(t) are the amplitude and the phase of {tilde over (Γ
)}ODT(t) respectively.
-
-
24. The apparatus of claim 19 where said display for generating said tomographic image of the fluid flow in said sample comprises generates a velocity image of said sample.
-
25. The apparatus of claim 19 where said processor determines a phase shift Δ
- f=Δ
ψ
/T from an average Doppler frequency shift, <
Δ
f>
, between sequential line scans using the equation;
where T is the time interval between sequential scans, and n is the number of sequential scans averaged.
- f=Δ
-
26. The apparatus of claim 25 where said processor determines a phase shift from an average Doppler frequency shift between sequential line scans comprises averaging at least three sequential line scans to increase the signal-to-noise ratio in images.
-
27. The apparatus of claim 19 where said processor determines a standard deviation of a Doppler power spectrum and said display generates said tomographic image of the fluid flow in said sample from said standard deviation of said Doppler power spectrum.
-
28. The apparatus of claim 27 where said processor determines said standard deviation of a Doppler power spectrum comprises determining said standard deviation from the equation:
-
where T is the time interval between sequential scans, P(ω
) is said Doppler power spectrum, and {overscore (ω
)} is the centroid value of a Doppler frequency shift.
-
-
29. The apparatus of claim 19 where said processor determines said corresponding phase at each pixel in a data window by extracting said phase by electronic demodulation of said detected interference fringes.
-
30. The apparatus of claim 19 where said scanner is characterized by a scanning rate, and wherein said processor adds π
- to said phase difference when a Doppler frequency difference related to said phase difference is greater than half of said scanning rate, because tan(φ
)=tan(φ
+π
), so that said Doppler frequency difference can exceed said scanning rate.
- to said phase difference when a Doppler frequency difference related to said phase difference is greater than half of said scanning rate, because tan(φ
-
31. The apparatus of claim 19 where said processor determines said ODT phase signals of said detected backscattered interference fringes through the use of a Hilbert transformation by determining said phase shift, Δ
- ψ
, in said computer at each pixel between sequential line scans by the Doppler frequency shift at each pixel between, ω
=Δ
ψ
/T where T is the time interval between successive line scans.
- ψ
-
32. The apparatus of claim 19 where said display for generating said tomographic image of the fluid flow in said sample comprises generates a structural image of said sample.
-
33. An apparatus for phase-resolved optical tomography of a sample comprising:
-
an interferometer a source of at least partially coherent radiation coupled to said interferometer, said at least partially coherent radiation characterized by a phase;
a modulator for phase modulating said radiation in said interferometer at a modulation frequency;
a scanner for scanning said sample with said source of at least partially coherent radiation through said interferometer in a sequence of pixel line scans, said sample having a fluid flow therein in which said phase of said at least partially coherent radiation is changed in response to said fluid flow at each pixel of each pixel line scan;
a detector for detecting interference fringes of said radiation backscattered from said sample into said interferometer;
a processor for processing said detected interference fringes to determine ODT phase signals of said detected backscattered interference fringes at each pixel of said pixel line scans to generate a corresponding rate of change of said phase at each pixel in a data window through the use of a Hilbert transformation at each pixel, and a display for generating a tomographic image of the fluid flow in said sample from said rate of change of said phase at each pixel.
-
-
34. An apparatus for phase-resolved optical tomography in a sample comprising:
-
an interferometer with a reference arm, a sample arm and a detection arm;
a source of at least partially coherent radiation coupled to said interferometer, said at least partially coherent radiation characterized by a phase;
a modulator for phase modulating said radiation in said interferometer at a modulation frequency;
a scanner for scanning said sample with said source of at least partially coherent radiation through said interferometer in a sequence of pixel line scans, said sample having a fluid flow therein so that said phase of said at least partially coherent radiation is changed in response to said fluid flow at each pixel of each pixel line scan;
a circular polarizer for circularly polarizing reflected radiation from said reference arm;
a linear polarizer for linearly polarized reflected radiation from said sample arm;
a coupler for combining said circularly polarized reflected radiation from said reference arm and said linearly polarized reflected radiation from said sample arm; and
a beam splitter for splitting radiation from said coupler into two polarized beams in said detection arm with a polarizing beam splitter with an axis of polarization set at 45°
to said linearly polarized reflected radiation from said sample arm,two corresponding detectors for detecting interference fringes of said radiation backscattered from said sample into said interferometer in said two polarized beams to generate two corresponding signals comprising in turn the real and imaginary parts of {tilde over (Γ
)}ODT(t) so that said phase is directly measured; and
a display for generating a tomographic image of the fluid flow in said sample from said phase at each pixel.
-
-
35. A method for optical tomography comprising:
-
providing a source of at least partially coherent radiation through an interferometer, said at least partially coherent radiation characterized by a phase;
providing a rapid-scanning optical delay line and an electro-optic modulator coupled to said rapid-scanning optical delay line;
phase modulating said radiation in said interferometer by means of said rapid-scanning optical delay line and electro-optic modulator in combination at a modulation frequency with a stable phase modulation;
scanning said sample with said source of at least partially coherent radiation through said interferometer;
detecting interference fringes of said radiation backscattered from said sample into said interferometer. - View Dependent Claims (36, 37, 38)
changing said phase of said at least partially coherent radiation in response to said fluid flow at each pixel of each pixel line scan;
processing said detected interference fringes to determine ODT phase signals of said detected backscattered interference fringes at each pixel of said pixel line scans to determine said corresponding phase at each pixel in a data window;
comparing said phase between corresponding pixels in two line scans to generate a difference between said phase at said two corresponding pixels in two line scans whereby speckle is substantially reduced; and
generating a tomographic image of the fluid flow in said sample from said phase difference at each pixel.
-
-
39. An apparatus for optical tomography comprising:
-
an interferometer;
a source of at least partially coherent radiation coupled to said interferometer, said at least partially coherent radiation characterized by a phase;
a rapid-scanning optical delay line;
an electro-optic modulator coupled to said rapid-scanning optical delay line to modulate said radiation in said interferometer at a modulation frequency with a stable phase modulation;
a scanner for scanning said sample with said source of at least partially coherent radiation through said interferometer; and
a detector to detect interference fringes of said radiation backscattered from said sample into said interferometer. - View Dependent Claims (40, 41, 42)
a processor to determine said corresponding phase at each pixel of said pixel line scans from said ODT phase signals of said detected backscattered interference fringes and to compare said phase between corresponding pixels in two line scans to generate a difference between said phase at said two corresponding pixels in two line scans whereby speckle is substantially reduced; and
a display to generate a tomographic image of the fluid flow in said sample from said difference at each pixel.
-
Specification