Method and apparatus for optical Doppler tomographic imaging of fluid flow velocity in highly scattering media

• US 5,991,697 A
• Filed: 12/31/1996
• Issued: 11/23/1999
• Est. Priority Date: 12/31/1996
• Status: Expired due to Term

• Current Assignee: California Institute of Technology
• Original Assignee: Regents of the University of California (University of California)
• Inventors: Nelson, John Stuart, Milner, Thomas Edward, Chen, Zhongping
• Primary Examiner(s): WACHSMAN, HAL D
• Application Number:
  US08/775,279
• Time in Patent Office: 1,057 Days
• Field of Search: 702/49, 702/19, 702/21, 702/22, 702/28, 702/30, 702/32, 702/40, 702/45, 702/50, 702/75, 702/76, 702/87, 702/88, 702/100, 702/159, 702/189, 365/528.17, 382/107, 382/128, 382/131, 382/133, 382/134, 600/425, 600/427, 600/441, 600/453-457, 73/861.25, 73/861, 73/861.02, 73/861.03, 73/863, 73/863.01-863.03, 73/601-603, 73/605, 73/606, 73/618, 367/7, 367/10, 356/337, 356/345, 356/346, 356/348, 356/349, 356/351, 356/352, 356/359-361, 356/39, 356/40, 356/301, 356/73.1, 356/27, 356/28, 356/28.5, 250/363.04, 250/550, 324/304-306, 324/310, 324/312
• US Class Current: 702/49
• CPC Class Codes: A61B 5/0261 using optical means, e.g. infrared light
• CPC Class Codes: A61B 5/412 Detecting or monitoring sepsis
• CPC Class Codes: G01P 5/26 by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave

Abstract

Optical Doppler tomography permits imaging of fluid flow velocity in highly scattering media. The tomography system combines Doppler velocimetry with high spatial resolution of partially coherent optical interferometry to measure fluid flow velocity at discrete spatial locations. Noninvasive in vivo imaging of blood flow dynamics and tissue structures with high spatial resolutions of the order of 2 to 10 microns is achieved in biological systems. The backscattered interference signals derived from the interferometer may be analyzed either through power spectrum determination to obtain the position and velocity of each particle in the fluid flow sample at each pixel, or the interference spectral density may be analyzed at each frequency in the spectrum to obtain the positions and velocities of the particles in a cross-section to which the interference spectral density corresponds. The realized resolutions of optical Doppler tomography allows noninvasive in vivo imaging of both blood microcirculation and tissue structure surrounding the vessel which has significance for biomedical research and clinical applications.

Claims

1. A method for tomographic imaging of a fluid flow in a highly scattering medium comprising:

providing a source of at least partially coherent radiation through an interferometer;

phase modulating said radiation in said interferometer at a modulation frequency;

continuously scanning a fluid flow sample with said source of at least partially coherent radiation through said interferometer, said fluid flow sample having a fluid flow therein and a structure in which said fluid flow is defined;

detecting interference fringes of said radiation backscattered from said sample into said interferometer; and

data processing Doppler frequency changes of said detected backscattered interference fringes with respect to said modulation frequency at each pixel of a scanned image to continuously measure the interference fringe intensities to obtain time dependent power spectra for each pixel location in a data window in a continuous scan from which a tomographic image of the fluid flow in and the structure of said scanned fluid flow sample is formed.

2. The method of claim 1 where scanning said fluid flow sample with a source of at least partially coherent radiation comprises scanning said fluid flow sample with a source of low coherence radiation.

3. The method of claim 1 where scanning said fluid flow sample comprises scanning said fluid flow sample with said source of radiation through a fiberoptic Michelson interferometer.

4. The method of claim 1 wherein said interferometer has a reference beam and a sample beam and scanning said fluid flow sample with said source of radiation through said interferometer further comprises matching polarity and equalize optical path length differences of said reference and said sample beams to optimize fringe contrast in said interferometer.

5. A method for tomographic imaging of a fluid flow in a highly scattering medium comprising:

providing a source of at least partially coherent radiation through an interferometer;

phase modulating said radiation in said interferometer at a modulation frequency;

scanning a fluid flow 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; and

data processing Doppler frequency changes of said detected backscattered interference fringes with respect to said modulation frequency at each pixel of a scanned image to form a tomographic image of said scanned fluid flow sample,

wherein said interferometer has a sample and a reference beam and wherein said reference beam is attenuated to reduce photon bunching noise and increase signal-to-noise ratio of interference fringes between said reference and a backscattered sample beam within said interferometer.

6. A method for tomographic imaging of a fluid flow in a highly scattering medium comprising:

providing a source of at least partially coherent radiation through an interferometer;

phase modulating said radiation in said interferometer at a modulation frequency;

scanning a fluid flow 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; and

data processing Doppler frequency changes of said detected backscattered interference fringes with respect to said modulation frequency at each pixel of a scanned image to form a tomographic image of said scanned fluid flow sample,

wherein said interferometer has a sample optical path and a reference optical path with a reference mirror disposed in said reference optical path, said radiation having a beam waist in said sample path which includes said fluid flow sample, said method further comprising dynamically maintaining a zero optical path length difference between said beam waist in said fluid flow sample and said reference mirror notwithstanding axial movement of said beam waist in said fluid flow sample by means of compensating axial movement of said reference mirror in said reference optical path as said fluid flow sample is scanned.

7. The method of claim 6 where scanning said fluid flow sample includes making a plurality of incremental axial scanning movements, and where maintaining said zero optical path length difference between said beam waist in said fluid flow sample and said reference mirror comprises compensating for each said incremental axial scanning movement with a corresponding axial movement in said reference mirror to maintain said optical path length changes in said sample and reference beams zero as said fluid flow sample is scanned.

8. An apparatus for optical Doppler tomographic imaging of fluid flow velocity in highly scattering medium comprising:

a source of at least partially coherent radiation;

an interferometer coupled to said source of radiation, said interferometer having a reference arm for a reference beam and a sample arm for a sample beam;

a modulator coupled to said interferometer for modulating optical path length difference in said reference and sample arms of said interferometer at a predetermined phase modulation frequency;

a scanner coupled to said interferometer for scanning a fluid flow sample with said sample beam of said interferometer;

a sensor coupled to said interferometer for detecting backscattered radiation received by said interferometer from said scanner to detect interference fringes within said interferometer; and

a data processor coupled to said sensor for processing signals from said sensor corresponding to said interference fringes established by said backscattered radiation in said interferometer and for controlling said scanner to generate tomographic images, wherein said data processor generates a power spectrum of interference fringe intensity at each pixel defined within said scan, and wherein said data processor simultaneously generates a tomographic structural image and a tomographic flow velocity image from said power spectrum at each pixel in said image.

9. The apparatus of claim 8 further comprising a polarization element coupled to said interferometer to match polarity of said reference and sample beams to optimize fringe contrast and equalize optical path length difference.

10. An apparatus for optical Doppler tomographic imaging of fluid flow velocity in highly scattering medium comprising:

a source of at least partially coherent radiation;

an interferometer coupled to said source of radiation, said interferometer having a reference arm for a reference beam and a sample arm for a sample beam;

a modulator coupled to said interferometer for modulating optical path length difference in said reference and sample arms of said interferometer at a predetermined phase modulation frequency;

a scanner coupled to said interferometer for scanning a fluid flow sample with said sample beam of said interferometer;

a sensor coupled to said interferometer for detecting backscattered radiation received by said interferometer from said scanner to detect interference fringes within said interferometer; and

a data processor coupled to said sensor for processing signals from said sensor corresponding to said interference fringes established by said backscattered radiation in said interferometer and for controlling said scanner to generate tomographic images,

where said reference arm of said interferometer includes a reference mirror, and where said sample beam has a beam waist, said apparatus further comprising means for maintaining a zero optical path length difference between said sample and reference beams as said beam waist of said sample beam is axially moved in said sample by compensating axial movement of said reference mirror.

11. An apparatus for optical Doppler tomographic imaging of fluid flow velocity in highly scattering medium comprising:

a source of at least partially coherent radiation;

an interferometer coupled to said source of radiation, said interferometer having a reference arm for a reference beam and a sample arm for a sample beam;

a modulator coupled to said interferometer for modulating optical path length difference in said reference and sample arms of said interferometer at a predetermined phase modulation frequency;

a scanner coupled to said interferometer for scanning a fluid flow sample with said sample beam of said interferometer;

a sensor coupled to said interferometer for detecting backscattered radiation received by said interferometer from said scanner to detect interference fringes within said interferometer; and

a data processor coupled to said sensor for processing signals from said sensor corresponding to said interference fringes established by said backscattered radiation in said interferometer and for controlling said scanner to generate tomographic images,

wherein said data processor generates a power spectrum of interference fringe intensity at each pixel defined within said scan, and

wherein said data processor simultaneously generates a tomographic structural image and a tomographic flow velocity image from said power spectrum at each pixel in said image.

12. A method for tomographic imaging of a fluid flow in a highly scattering medium comprising:

scanning a fluid flow sample with a source of at least partially coherent radiation through an interferometer to generate backscattered radiation in said interferometer;

phase modulating said radiation in said interferometer at a modulation frequency;

measuring spectral periods and phase change rates of modulation of an interference spectral density of said backscattered radiation in said interferometer; and

generating a tomographic image of said scanned fluid flow sample from said spectral periods and phase change rates to simultaneously obtain positions and velocities of objects moving within said fluid flow sample.

13. The method of claim 12 where scanning said fluid flow sample comprises scanning said fluid flow sample with said source of radiation through a fiberoptic Michelson interferometer.

14. The method of claim 12 wherein said interferometer has a reference beam and a sample beam and scanning said fluid flow sample with said source of radiation through said interferometer further comprises matching polarity and equalizing optical path length differences of said reference and said sample beams to optimize fringe contrast in said interferometer.

15. The method of claim 12 wherein said interferometer has a sample and a reference beam and wherein said reference beam is attenuated to reduce photon bunching noise and increase signal-to-noise ratio of interference fringes between said reference and a backscattered sample beam within said interferometer.

16. The method of claim 12 wherein said interferometer has a sample and reference beam and wherein scanning said fluid flow sample with modulated radiation comprises modulating the optical path lengths of said radiation in said reference and sample beams of said interferometer.

17. The method of claim 12 wherein said modulation of said interference spectral density of said radiation is characterized by a plurality of frequencies and wherein measuring said frequencies and phase change rates of said modulation of said interference spectral density of said radiation in said interferometer is performed by measuring said spectral periods and phase change rates in an optical spectrum analyzer coupled to said interferometer to receive interference spectral density of said backscattered radiation in said interferometer.

18. The method of claim 12 wherein generating a tomographic image of said scanned fluid flow sample from said spectral periods and phase change rates to simultaneously obtain positiond and velocities of objects moving within said fluid flow sample, is obtained from said interference spectral density by obtaining a position of each object, O_i, in said fluid flow sample by measurement of a corresponding spectral period Γ

(Δ

_l)=C/2Δ

_l, of a frequency in said interference spectral density, where c is the speed of light, and Δ

_l is the path length difference of the O_i object at t=0, where t is the time of measurement, and by measurement of a corresponding rate of phase change, dΦ

_i /dt=4π

ν

V_i /c, of said spectral period in said inteference spectral density, where V_i is the velocity of the O_i object at t=0, ν

is the frequency of said modulation, and c is the velocity of light.

19. An apparatus for optical Doppler tomographic imaging of fluid flow velocity in highly scattering medium comprising:

a source of at least partially coherent radiation;

an interferometer coupled to said source of radiation, said interferometer having a reference arm for a reference beam and a sample arm for a sample beam, and an output;

a modulator coupled to said interferometer for modulating optical path lengths in said reference and sample arms of said interferometer at a predetermined phase modulation frequency;

a probe coupled to said interferometer for exposing a fluid flow sample with said sample beam of said interferometer;

an optical spectrum analyzer coupled to said output of said interferometer for detecting backscattered radiation received by said interferometer from said probe to detect and analyze a plurality of modulated optical spectral interference densities produced by said interferometer; and

a data processor coupled to said optical spectrum analyzer for processing signals from said probe for spectral scanning said fluid flow sample corresponding to said plurality of modulated optical spectral interference densities established by said backscattered radiation in said interferometer to generate tomographic images from said plurality of interference spectral densities.

20. The apparatus of claim 19 further comprising a polarization element coupled to said interferometer to match polarity and equalize optical path lengths of said reference and sample beams to optimize fringe contrast.

21. An apparatus for optical Doppler tomographic imaging of fluid flow velocity in highly scattering medium comprising:

a source of at least partially coherent radiation;

an interferometer coupled to said source of radiation, said interferometer having a reference arm for a reference beam and a sample arm for a sample beam, and an output;

a modulator coupled to said interferometer for modulating optical path lengths in said reference and sample arms of said interferometer at a predetermined phase modulation frequency;

a probe coupled to said interferometer for exposing a fluid flow sample with said sample beam of said interferometer;

an optical spectrum analyzer coupled to said output of said interferometer for detecting backscattered radiation received by said interferometer from said probe to detect and analyze modulated spectral interference densities produced by said interferometer; and

a data processor coupled to said optical spectrum analyzer for processing signals from said probe corresponding to said interference spectral densities established by said backscattered radiation in said interferometer to generate tomographic images,

wherein said optical spectrum analyzer measures spectral periods and phase change rates of said modulation of said interference spectral densities of said radiation in said interferometer.

22. The apparatus of claim 21 where said data processor generates a tomographic image of said fluid flow sample from said spectral periods and phase change rates to simultaneously obtain positions and velocities of objects moving within said fluid flow sample from said interference spectral densities by obtaining position of each object, O_i, in said fluid flow sample by measurement of a spectral period Γ

(Δ

_l)=c/2Δ

_l, of a frequency in said interference spectral densities, where c is the speed of light, and Δ

_l is the path length difference of the O_i object at t=0, where t is the time of measurement, and by measurement of a corresponding rate of phase change, dΦ

_j /dt=4π

ν

V_i /c, of said frequency in said interference spectral density, where V_i is the velocity of the O_i object at t=0, ν

is the frequency of said modulation, and c is the speed of light.