WO2008048263A1 - Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus - Google Patents
Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus Download PDFInfo
- Publication number
- WO2008048263A1 WO2008048263A1 PCT/US2006/040932 US2006040932W WO2008048263A1 WO 2008048263 A1 WO2008048263 A1 WO 2008048263A1 US 2006040932 W US2006040932 W US 2006040932W WO 2008048263 A1 WO2008048263 A1 WO 2008048263A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- blood flow
- magnetic field
- optical
- doppler
- detecting
- Prior art date
Links
- 238000012014 optical coherence tomography Methods 0.000 title claims abstract description 90
- 230000003287 optical effect Effects 0.000 title claims abstract description 63
- 238000003325 tomography Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 title claims description 40
- 108010054147 Hemoglobins Proteins 0.000 title claims description 36
- 102000001554 Hemoglobins Human genes 0.000 title claims description 36
- 238000012285 ultrasound imaging Methods 0.000 title description 3
- 230000005291 magnetic effect Effects 0.000 claims abstract description 96
- 230000017531 blood circulation Effects 0.000 claims abstract description 71
- 238000003384 imaging method Methods 0.000 claims abstract description 40
- 238000002604 ultrasonography Methods 0.000 claims abstract description 28
- 239000000523 sample Substances 0.000 claims description 51
- 210000004369 blood Anatomy 0.000 claims description 21
- 239000008280 blood Substances 0.000 claims description 21
- 230000033001 locomotion Effects 0.000 claims description 15
- 230000003595 spectral effect Effects 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 9
- 210000004204 blood vessel Anatomy 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 claims description 5
- 229910000859 α-Fe Inorganic materials 0.000 claims description 4
- 235000017284 Pometia pinnata Nutrition 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 claims description 2
- 230000005855 radiation Effects 0.000 claims description 2
- 210000005166 vasculature Anatomy 0.000 claims 2
- 240000007653 Pometia tomentosa Species 0.000 claims 1
- 239000012472 biological sample Substances 0.000 claims 1
- 230000036962 time dependent Effects 0.000 claims 1
- 210000003743 erythrocyte Anatomy 0.000 abstract description 23
- 230000007246 mechanism Effects 0.000 abstract description 2
- 210000001519 tissue Anatomy 0.000 description 20
- 239000002245 particle Substances 0.000 description 11
- 239000012530 fluid Substances 0.000 description 9
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 8
- 206010028980 Neoplasm Diseases 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000005298 paramagnetic effect Effects 0.000 description 7
- 201000011510 cancer Diseases 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000004907 flux Effects 0.000 description 5
- 239000013307 optical fiber Substances 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 239000002872 contrast media Substances 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 4
- 238000005534 hematocrit Methods 0.000 description 4
- 238000001727 in vivo Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000004005 microsphere Substances 0.000 description 4
- 230000010287 polarization Effects 0.000 description 4
- 238000002583 angiography Methods 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 3
- 239000011521 glass Substances 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 230000002792 vascular Effects 0.000 description 3
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N Iron oxide Chemical compound [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 description 2
- 208000006787 Port-Wine Stain Diseases 0.000 description 2
- 206010060862 Prostate cancer Diseases 0.000 description 2
- 208000000236 Prostatic Neoplasms Diseases 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000006213 oxygenation reaction Methods 0.000 description 2
- 230000010412 perfusion Effects 0.000 description 2
- 210000002307 prostate Anatomy 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- 206010067193 Naevus flammeus Diseases 0.000 description 1
- 240000009305 Pometia pinnata Species 0.000 description 1
- 208000000453 Skin Neoplasms Diseases 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000008081 blood perfusion Effects 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 238000001246 colloidal dispersion Methods 0.000 description 1
- 239000013068 control sample Substances 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 239000008367 deionised water Substances 0.000 description 1
- 229910021641 deionized water Inorganic materials 0.000 description 1
- 108010002255 deoxyhemoglobin Proteins 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000005292 diamagnetic effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000009552 doppler ultrasonography Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 208000002026 familial multiple nevi flammei Diseases 0.000 description 1
- 229940079405 ferumoxides Drugs 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000001631 haemodialysis Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000003709 heart valve Anatomy 0.000 description 1
- 230000000322 hemodialysis Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 239000004816 latex Substances 0.000 description 1
- 229920000126 latex Polymers 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000011553 magnetic fluid Substances 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 239000002105 nanoparticle Substances 0.000 description 1
- 239000002078 nanoshell Substances 0.000 description 1
- 230000017074 necrotic cell death Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000008807 pathological lesion Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 210000005259 peripheral blood Anatomy 0.000 description 1
- 239000011886 peripheral blood Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000002428 photodynamic therapy Methods 0.000 description 1
- 230000010399 physical interaction Effects 0.000 description 1
- 210000002381 plasma Anatomy 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 235000020043 port wine Nutrition 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000979 retarding effect Effects 0.000 description 1
- 230000002207 retinal effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 201000000849 skin cancer Diseases 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 230000026683 transduction Effects 0.000 description 1
- 238000010361 transduction Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 238000000827 velocimetry Methods 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/0515—Magnetic particle imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4416—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
Definitions
- the present invention relates in general to the art of medical diagnostic imaging and in particular to imaging blood flow using magneto-motive optical Doppler tomography (MM-ODT), Optical Coherence Tomography, or Ultrasound, which combines an externally applied temporally oscillating high-strength magnetic field with ODT, OCT, or Ultrasound to detect erythrocytes moving according to the field gradient
- MM-ODT magneto-motive optical Doppler tomography
- Optical Coherence Tomography or Ultrasound
- One common sensing technique involves the use of ultrasound.
- ultrasound is directed into the body of the patient and tiny particles such as red blood cells, which are suspended in the blood plasma, scatter the ultrasonic energy back towards the receiver or transducer.
- the transducer then converts the back-scattered ultrasonic energy into an electrical signal that is processed in some known manner to determine the presence of a flow and an estimate of the flow velocity.
- Magnetic resonance imaging is based on an imaging technique for magnetically exciting nuclear spins in a subject positioned in a static magnetic field by applying a radio-frequency (RF) signal of the Larmor frequency, and reconstructing an image
- Magnetic resonance angiography provides detailed angiographic images of the body in a non-invasive manner.
- MRA Magnetic resonance angiography
- conventional MRA which does not use contrast agents, magnetic resonance signal from flowing blood is optimized, while signal from stationary blood or tissue structures is suppressed.
- contrast-enhanced MRA a contrast agent is injected into the blood stream to achieve contrast between flowing blood and stationary tissue.
- EPI echo planar imaging
- fMRI Functional magnetic resonance imaging
- OCT optical coherence tomography
- Optical coherence tomography is a technology that allows for noninvasive, cross-sectional optical imaging of biological media with high spatial resolution and high sensitivity.
- OCT is an extension of low-coherence or white-light interferometry, in which a low temporal coherence light source is utilized to obtain precise localization of reflections internal to a probed structure along an optic axis. This technique is extended to enable scanning of the probe beam in the direction perpendicular to the optic axis, building up a two-dimensional reflectivity data set, used to create a cross-sectional gray-scale or false- color image of internal tissue backscatter.
- OCT uses the short temporal coherence properties of broadband light to extract structural information from heterogeneous samples such as biologic tissue.
- OCT has been applied to imaging of biological tissue in vitro and in vivo.
- Systems and methods for substantially increasing the resolution of OCT and for increasing the information content of OCT images through coherent signal processing of the OCT interferogram data have been developed to provide cellular resolution (i.e., in the order of 5 micrometers).
- numerous advancements in OCT have been reported including real-time imaging speeds.
- ODT Optical Doppler tomography
- OCT optical coherence tomography
- the principle employed in ODT is very similar to that used in radar, sonar and medical ultrasound.
- ODT uses a low coherence or broadband light source and optical interferometer to obtain high spatial resolution gating with a high speed scanning device such as a conventional rapid scanning optical delay line (RSOD) to perform fast ranging of microstructure and particle motion detection in biological tissues or other turbid media.
- RSOD rapid scanning optical delay line
- the Doppler frequency shift is proportional to the cosine of the angle between the probe beam and the scatterer's flow direction. When the two directions are perpendicular, the Doppler shift is zero. Because a priori knowledge of the Doppler angle is not available, and conventional intensity OCT imaging provides a low contrast image of microvasculature structure, detecting small vessels with slow flow rates is difficult. However, the Doppler angle can be estimated by combining Doppler shift and Doppler bandwidth measurements. (See for example, Piao DQ, Zhu Q, Quantifying Doppler Angle and Mapping
- 5,991,697 describe a method and apparatus for Optical Doppler Tomographic imaging of a fluid flow in a highly scattering medium comprising the steps of scanning a fluid flow sample with an optical source of at least partially coherent radiation through an interferometer, which is incorporated herein by reference).
- MM-OCT magneto-motive OCT
- fMRI Functional magnetic resonance imaging
- T2 differentiating relaxation times
- the present invention is a novel extension of MM-OCT to image hemoglobin in blood erythrocytes.
- the approach requires no exogenous contrast agent to detect blood flow and location.
- the invention is a method and apparatus of imaging a blood flow using
- Doppler optical coherence tomography which comprises an externally applied temporally oscillating high-strength magnetic field with a Doppler optical tomography system to detect hemoglobin moving according to the magnetic field gradient.
- the invention comprises a solenoid cone-shaped ferrite core with an extensively increased magnetic field strength at the tip of the core, hi a further embodiment, the invention comprises focusing the magnetic force on targeted samples.
- Another embodiment is an apparatus for imaging blood flow, comprising: (1) a magnetomotive optical Doppler tomography (MM-ODT) imaging system; (2) an oscillating magnetic field applied to the moving blood flow; (3) a light source (4) a rapid scanning optical delay line in a reference arm and aligned such that no phase modulation is generated when the group phase delay is scanned at 4 kHz; and (5) a hardware in-phase and a quadrature demodulator with high bandpass filters.
- MM-ODT magnetomotive optical Doppler tomography
- Another aspect of the invention is a method for imaging a blood flow, comprising applying a magnetic field to the blood flow, wherein the blood flow comprising a plurality of hemoglobin molecules and wherein the magnetic field interacts with the hemoglobin molecules to cause a change in the blood flow; and detecting the blood flow by detecting the change in the blood flow caused by the interaction with the hemoglobin molecules, wherein the change is detected using a magnetomotive optical Doppler tomography imaging system.
- Another embodiment of the method for imaging a blood flow comprises temporally oscillating the magnetic field.
- Another embodiment of the invention is an apparatus for imaging a blood flow comprising a magnetic field generator for applying a magnetic field to the blood flow, wherein the blood flow comprises hemoglobin molecules; and an ultrasound detection system for detecting the blood flow while it is in the presence of the magnetic field.
- FIG. 1 is a schematic diagram of the MM-ODT system.
- FIG. 2 is a schematic diagram of the probe beam, flow sample and solenoid coil.
- FIG. 3a and 3* are OCT and ODT images of a stationary turbid solution without an external magnetic field, respectively.
- FIG. 3c and FIG. 3d are OCT and ODT Attorney Docket: 6100-043 Customer No. 29335 images with a 50 Hz magnetic field, respectively.
- the white bar represents 200 ⁇ m, accordingly.
- FIGS. 4a-4d are M-mode ODT images of the diluted deoxygenated blood flow (18% hematocrit) without and with an external magnetic field.
- FIG. 4a and FIG. 46 are ODT images of 5 mm/s blood flow without and with a 5 Hz magnetic field, respectively.
- FIG. 4c and FIG. 4d are ODT images of 30 mm/s blood flow without and with a 50 Hz magnetic field, respectively.
- the Black bar indicates 200 ⁇ m, accordingly.
- FIG. Sa and FIG. Sb axe the Doppler frequency shift profiles without an external magnetic field, and with a 50 Hz magnetic field, respectively.
- FIG. 6 is a block diagram illustrating an exemplary phase sensitive OCT system.
- FIG. 7 is a graph of the optical path length change in the DP-OCT system.
- FIG. 8 is a schematic diagram of Magneto Motive Ultrasound for Blood.
- FIG. 9 is an ultrasound image of blood exposed to an oscillating magnetic field.
- FIG. 10 is a schematic diagram of Ultrasonography and OCT coupled with magnetic field generator.
- MM-ODT Doppler imaging of blood flow using an external oscillating magnetic field is described below.
- RBCs red blood cells
- MM-ODT allows imaging of blood flow and velocity.
- the controlled and increased Doppler frequency in MM-ODT provides an investigational tool to study in vivo blood transport, as shown in the article Hemoglobin Contrast in Magnetomotive Optical Doppler Tomography, Opt. Lett. 31, 778-780 (2006), herein incorporated by reference.
- microstructure of the blood flow and flow velocity information are all encoded in the interferogram of a Doppler OCT system. It should be readily apparent to those skilled in the optical arts, that different OCT systems and different OCT information
- Time domain OCT requires a mirror in the reference arm scanning at a constant velocity
- spectral domain OCT includes swept source OCT and Fourier domain OCT.
- An example of time-domain Doppler OCT is provided below.
- the OCT light source 12 comprises a super luminescent diode, which is used as the low coherence light source, hi one embodiment, light source 12 is centered at 1.3 ⁇ m with a bandwidth of 90 run.
- Light from source 12 is coupled into a single-mode optical fiber based interferometer 20 by circulator 14, where the interferometer 20 can provide 1 mW of optical power at the sample 60.
- Light is split into a reference arm 22 and sample arm 24 by a 2 x 2 splitter 16.
- a rapid- scanning optical delay (RSOD) line 26 is coupled to the reference arm 22.
- RSOD rapid- scanning optical delay
- the rapid-scanning optical delay line 26 is aligned such that no phase modulation is generated when the group phase delay is scanned at ⁇ 4 kHz.
- a collimated beam 36 is redirected to sample 60 by two galvanometers 30 that permit three-dimensional scanning.
- the galvanometers can be an X-scanner and a Y-scanner.
- the sample 60 can be internal or external to the body, where the probe beam is focused by an objective lens 32.
- the objective lens 32 yields a 10- ⁇ m spot at the focal point.
- Phase modulation can be generated using an electro-optical waveguide phase modulator, which can produce a carrier frequency ( ⁇ 1 MHz).
- the magnetic field generator 100 is in proximity to the sample 60.
- a dual-balanced photodetector 34 is coupled to the 2 x 2 splitter 16 and the circulator 14.
- the photodetector 34 of a 80 MHz bandwidth reduces the light source noise from the OCT interference signal.
- a hardware in-phase and quadrature demodulator 40 with high/bandpass filters 42 and low/bandpass filters 44 improves imaging speed.
- Doppler information was calculated with the Kasai autocorrelation velocity estimator.
- Labview software 50 National Instruments, Austin, TX
- the maximum frame rate of the MM-ODT system 10 was 16 frames per second for a 400 x 512 pixel sized image.
- the Doppler frequency shift can be determined with the use of a short time fast Fourier transform (STFFT).
- a phase-resolved technique can determine the Doppler frequency shift to decouple the
- OCT differential phase optical coherence tomography
- spectral domain phase-sensitive OCT can be used to determine the Doppler frequency shift, as readily determined by one skilled in the optical arts.
- FIG. 2 shows an example of the magnetic field generator 80 with a capillary tube 90.
- the magnetic field generator 80 includes a solenoid coil 82 (Ledex 4EF) with a cone-shaped ferrite core 84 at the center and driven by a current amplifier 86 supplying up to 960 W of power.
- the magnetic field generator 80 can be placed underneath the sample 60 during MM-ODT imaging.
- the sinusoidal current can vary the magnetic force applied to the capillary tube 90 in order to introduce movement of magnetic fluids, which include red blood cells that contain hemoglobin.
- the probe beam is oriented parallel to the gradient of the magnetic field's strength.
- the material parameter characterizing magnetic materials is the magnetic volume susceptibility, ⁇ .
- Magnetic volume susceptibility is dimensionless in SI units and is defined by the equation M — ⁇ H, where M is the magnetization at the point in question and H is the local density of the magnetic field strength.
- M is the magnetization at the point in question
- H is the local density of the magnetic field strength.
- Hemoglobin's high iron content, due to four Fe atoms in each hemoglobin molecule, and the large concentration of hemoglobin in human red blood cells give Hemoglobin magnetomotive effects in biological tissue.
- the magnetic volume susceptibility of the hemoglobin molecule consists of a paramagnetic component due to the electron spins of the four iron atoms. The paramagnetic susceptibility is given by the Curie Law,
- ⁇ 0 permeability of free space and has the value 4 ⁇ xlO 7 H/m
- N p is the volume density of paramagnetic iron atoms in hemoglobin
- the Bohr magneton, ⁇ B : 9.274xl0 "24 J/T
- Boltzmann's constant, k 1.38x 10 ⁇ 23 J/K
- T is the absolute Attorney Docket: 6100-043 Customer No. 29335 temperature (K).
- the calculated susceptibility of a RBC is about HxIO "6 assuming a 90% concentration of hemoglobin per RBC.
- the calculated susceptibility of a RBC is dependent on the oxygenation of the hemoglobin. The calculations can be adjusted accordingly, depending on the oxygenation of the RBC, which can be measured by known techniques.
- a RBC placed in a magnetic field gradient experiences forces and torques that tend to position and align it with respect to the field's direction.
- the magnetic force, in the direction of the probing light z, is given by
- V is the particle volume
- B is the magnitude of the magnetic flux density
- a ⁇ is the difference between the susceptibility of the particle and the surrounding medium.
- I R and Is are the back scattered intensities from the reference and sample arms, respectively
- f o is the fringe carrier frequency
- X 0 is the center wavelength of the light source
- z(t) is the RBC displacement.
- the displacement, z(t) is dominated a constant acceleration which can be, however, ignored in confined models (i.e. blood vessel or capillary tube) with assumptions that, first, the probing area is much smaller than magnetic field area, and that secondly probing time starts after steady states of inner pressures.
- Expansion of the right-hand side of Eq. (3) using Bessel functions gives
- Jk(m) is the Bessel function of the first kind of order k for argument m which is 4 ⁇ A/ ⁇ 0 .
- the amplitude of the k* sideband is proportional to j k (m).
- M-mode OCT/ODT images of a capillary glass tube filled with a stationary turbid solution with and without an external magnetic field as a control sample were recorded, as shown in FIGS. 3-4.
- a 750 ⁇ m-inner diameter glass capillary tube 200 was placed perpendicularly to the probing beam 70, as shown in FIG. 2.
- Fluids used for flow studies were injected through the tube at a constant flow rate controlled by a dual-syringe pump (Harvard Apparatus 11 Plus, Holliston, MA) with ⁇ 0.5% flow rate accuracy.
- FIGS. 3a and 3b show M-mode OCT and ODT images without any external magnetic field, respectively.
- the ODT image in FIG. 3b contains small random phase fluctuations due to ambient vibration through the optical path.
- FIGS. 3c and 3d show M-mode OCT and ODT images with a 50 Hz external magnetic field, respectively. No distinguishable Doppler shift could be observed in the ODT image FIG. 3d indicating no interaction between the external magnetic field and the moving microspheres.
- Deoxygenated blood was extracted from the vein of a human male's left arm, and diluted with saline. During preparation, blood was not exposed directly to air so as to remain deoxygenated. To simulate flow, blood was injected through the capillary tube 200 by a syringe pump at a relatively constant flow rate. As FIG. 4 shows, the oscillating Doppler frequency shift, resulting from RBC movement, could be observed at two different flow rates (5 and 30 mm/s ⁇ Because the flow direction was nearly perpendicular to the probing beam no significant Doppler frequency shift was distinguishable at the 5 mm/s flow rate FIG. 4a without any external magnetic field. In the case of the high blood flow rate of 30 mm/s, as shown at FIG. 4c, the Doppler frequency shift caused by the blood flow could be observed. And for maximum contrast enhancement, the probe beam can be directed parallel to the gradient of the magnetic field's strength. However, application of a 50 Hz magnetic
- Doppler frequency shift profiles were calculated from the ODT images by averaging 20 lines at a selected depth indicated by horizontal arrows, as shown in FIGS. 4a and 4b.
- FIG. 5a indicates no significant Doppler frequency shift over a 100 ms time period
- FIG. 5b displays ⁇ 200 to 300 Hz Doppler frequency shifts oscillating 20 times over 100 ms (200 Hz).
- the Doppler frequency shift indicates that RBCs physically move into and away from the incident light while passing through the external magnetic field depending on whether their magnetic properties are paramagnetic or diamagnetic, as shown in FIG. 4, and that the 200 Hz oscillation of the Doppler frequency shift correlates with the 50 Hz magnetic field, as shown in FIG. 5.
- Frequencies 4 times higher than that of the external magnetic field (fa n ) can be observed.
- the frequency of the force on paramagnetic targets was twice that of the magnetic flux density; therefore, a 50 Hz B field displaces the targets at 100 Hz.
- the fringe signal Eq.
- the invention is a new investigational tool to study in vivo blood transport and the first implementation of MM-ODT for improved Doppler imaging of blood flow using an external oscillating magnetic field introducing a mechanical movement of RBCs during blood flow by a temporally oscillating high-strength magnetic field.
- MM-ODT to allow imaging of tissue function in a manner similar to functional magnetic resonance images (f- MRI) of deoxygenated blood hi organs, when the sample arm of the MM-ODT system is coupled to a probe (not shown).
- probes are generally known in the arts, such as endoscopic probes, catheter probes, and the like.
- the MM-ODT can be used for Port- Wine vessel mapping and
- the MM-ODT can be used to detect blood vessel location and size for cancer and port-wine stains, since these conditions are characterized by blood vessel growth and increases in hemoglobin content. Accordingly, other blood vessel detection for tissue abnormality identification can be envisioned with this invention. Generally, the MM- ODT can be used wherever blood flow detection is necessary in operations, chemotherapy, hemodialysis, and the like.
- a spectral domain phase sensitive OCT system 100 can be used to image the blood flow and determine velocity information, when an oscillating magnetic field is applied to the blood flow, as shown in FIG. 6.
- Spectral domain phase sensitive OCT generally uses a broadband light source in a general interferometric setup, where the mirror in the reference typically does not move or does not require a rapid scanning delay line.
- a spectral interferogram is detected, in which each A-scan is entirely encoded.
- a Fourier transformation of the A-scan and velocity distribution of the blood flow can be extracted by known methods.
- the blood flow profile is the envelope of a calculated A-scan through Fourier transformation on the spectral interferogram.
- This OCT system 100 is only an example of one OCT imaging modality which can be used to image blood flow with a temporally oscillating magnetic field, and is not intended to suggest any limitation on the scope of OCT architectures applicable to the invention.
- the OCT system 100 includes a general-purpose computing device in the form of a computer 101 and includes a magnet control 114 and a magnet 116.
- Light energy is generated by a light source 117.
- the light source 117 can be a broadband laser light source coupled into optical fiber emitting light energy over a broad range of optical frequencies.
- the wavelength range can be from about 400 nanometers to about 1400 nanometers. Longer wavelengths (>600 run) can be used for deeper scanning.
- the light source emits light having a wavelength near the infrared spectrum to identify hemoglobin for OCT imaging, which places hemoglobin in motion and increases optical scattering of the hemoglobin.
- the light energy can be emitted over a multiplicity of optical wavelengths, frequencies, and pulse durations to achieve OCT imaging. As used
- optical fiber can refer to glass or plastic wire or fiber.
- Optical fiber is indicated on FIG. 6 as lines connecting the various blocks of the figures. Where light energy is described as “passing,” “traveling,” “returning,” “directed,” or similar movement, such movement can be via optical fiber.
- a fraction of the generated light energy passes from the light source 117 into an optical spectrum analyzer 118.
- the optical spectrum analyzer 118 measures optical frequency as the light energy is emitted from the light source 117 as a function of time.
- the ' optical spectrum analyzer 118 samples a portion of the light emitted by the light source 117.
- the optical spectrum analyzer 118 monitors the power spectral density of light entering the splitter 119.
- the remaining fraction of light energy from the light source 117 passes into a splitter 119.
- the splitter 119 can be a device with four ports, with Port 1 allowing light energy to enter the splitter 119.
- Ports 2 and 3 allow light energy to leave and re-enter the splitter 119 to the reference reflector 120 and OCT probe 122, respectively.
- Port 4 allows light energy to leave the splitter 119 to coupling lens 124.
- the splitter 119 couples the light into Port 1.
- the splitter 119 divides the light according to a pre-determined split ratio selected by a user.
- the split ratio can be 50/50 wherein half of the light energy entering the splitter 119 at Port 1 exits the splitter 119 through Port 2 and half exits the splitter 119 through Port 3.
- the split ratio can be 60/40 wherein 60% of the light energy passes through Port 2 and 40% of the light energy passes through Port 3.
- a fraction of the light energy (determined by the split ratio) that exits the splitter 119 through Port 2 travels to a reference reflector surface 120.
- the light energy is reflected from the reference reflector surface 120 back to the splitter 119 into Port 2.
- the reference reflector 120 can be a planar metallic mirror or a multilayer dielectric reflector with a specified spectral amplitude/phase reflectivity.
- the remaining fraction of light that entered splitter 119 through Port 1 exits splitter 119 through Port 3 and enters an OCT probe 122.
- the OCT probe 122 can be a turbine-type catheter as described in Patent Cooperation Treaty application PCT/US04/12773 filed 4/23/2004 which claims priority to U.S.
- the OCT probe 122 can be located within a subject 121 to allow light reflection off of subject's 121 blood flow.
- the light energy that entered OCT probe 122 is reflected off of the blood flow of subject 121. once an oscillating magnetic field has been temporally applied by magnet 116. The reflected light energy passes back through the OCT probe 122 into the splitter 119 via Port 3. The reflected light energy that is returned into Port 2 and Port 3 of the splitter 119 recombines and interferes according to a split ratio. The light recombines either constructively or destructively, depending on the difference of pathlengths. A series of constructive and destructive combinations of reflected light create an interferogram (a plot of detector response as a function of optical path length difference). Each reflecting layer from the subject 121 and the blood flow will generate an interferogram.
- the splitter 119 can recombine light energy that is returned through Port 2 and Port 3 so that the light energies interfere.
- the light energy is recombined in the reverse of the split ratio. For example, if a 60/40 split ratio, only 40% of the light energy returned through Port 2 and 60% of the light energy returned through Port 3 would be recombined.
- the recombined reflected light energy is directed out Port 4 of the splitter 119 into a coupling lens 137.
- the coupling lens 137 receives light from the output of the splitter 119 and sets the beam etendue (beam diameter and divergence) to match that of the optical spectrometer 138.
- the coupling lens 137 couples the light into an optical spectrometer 138.
- the optical spectrometer 138 can divide the recombined reflected light energy light into different optical frequencies and direct them to different points in space which are detected by a line scan camera 139.
- the line scan camera 139 performs light to electrical transduction resulting in digital light signal data 108.
- the digital light signal data 108 is transferred into the computer 101 via the OCT input interface 111. Interface between the line scan camera 139 and computer 101 can be a Universal Serial Bus (USB) 3 or the like.
- the digital light signal data 108 can be stored in the mass storage device 104 or system memory 112 and utilized by the image construction software 106 and the Labview image construction software 107.
- the image construction software 106 can generate an image of the blood flow of subject 121 from the light signal data 108, by receiving light signal data 108 generating amplitude and phase data.
- the amplitude and phase data (optical path length difference (c ⁇ ) or optical time-delay ( ⁇ )) can be separated into discrete channels and a plot of intensity vs. depth (or amplitude vs. depth) can be generated for each channel.
- Such plot is known as an A-scan, where the composition of all the A-scans can comprise one image.
- movement image construction software 107 generates an image of the movement of the hemoglobin from the light signal data 108.
- the movement image construction software 107 receives light
- additional information can be extracted from the light signal data to generate additional images.
- the light signal data can be further processed to generate a Stokes parameter polarimetric image when used in conjunction with polarization detectors and polarizing lenses to extract polarization data from the light signal 108, as readily known to one skilled in the art of optical coherence tomography.
- the differential phase OCT image data is shown in FIG. 7, indicating the optical path difference with the magnetic field. Filters can be added to reduce the noise of the signal generated by the magnetic field.
- the phase sensitive OCT system 100 can be configured for swept source OCT, which is a different type of spectral domain OCT.
- swept source OCT a tunable laser source replaces the broadband laser light source 117.
- the scanning rate can be at wavelengths of 800nm ⁇ 1500nm.
- the reference reflector surface 120 is in-line with sample path 120.
- the optical spectrometer 125 and line scan camera 126 are replaced with a general photodetector.
- Ultrasonography is the ultrasound-based diagnostic imaging technique used to visualize muscles and internal organs, their size, structures and any pathological lesions. "Ultrasound" applies to all acoustic energy with a frequency above human hearing (20,000 Hertz or 20 kilohertz). Typical diagnostic sonography scanners operate in the frequency range of 2 to 13 megahertz, hundreds of times greater than this limit. The choice of frequency is a trade-off between the image spatial resolution and penetration depth into the patient, with lower frequencies giving less resolution and greater imaging depth. Doppler ultrasonography uses the Doppler Effect to assess whether blood is moving towards or away from a probe, and its relative velocity.
- an ultrasound probe 212 is coupled with the magnetic field generator 100, as shown in FIG. 8.
- Ultrasound 210 is directed into the body of the patient by known techniques and the moving red blood cells backscatter the ultrasonic energy back towards the transducer of the ultrasound.
- the oscillating magnetic field generated by the magnetic field generator 100 increases the contrast of the ultrasonic energy 210 received from the red blood cells.
- the transducer then converts the back- scattered ultrasonic energy 210 into an electrical signal that is processed in some known manner to determine an estimate of the flow.
- An enhanced ultrasound image is produced, as displayed in FIG. 9.
- a rectal ultrasound probe is coupled with a magnet to evaluate the prostate gland for cancer.
- ultrasound is used for prostate cancer screening; however, the approach provides poor sensitivity and specificity.
- all cancers are known in the art to be highly vascular, so then the application of the magnetic field by generator 100 coupled with ultrasound enhances the contrast available from the endogenous RBCs in the prostate for cancer detection at an earlier stage. It is generally known in the art that cancers have enhanced metabolic properties compared to normal tissues, so then cancerous cells have higher oxygen content from hemoglobin and a greater concentration of deoxygenated hemoglobin compared to normal tissues.
- An exemplary ultrasound image for prostate cancer screening is shown in FIG. 9.
- MM-ODT, OCT, and ultrasound coupled with an oscillating magnetic field all can be used in clinical management of patients who need microvasculature monitoring, as shown in FIG. 10.
- the images could monitor and determine tissue profusion and viability before, during and after operation procedures.
- the images could be used to detect oxygenated and deoxygenated blood supply, detect blood vessel location, sizes of cancer, vascular tissue abnormality identification, and port-wine stain.
- the images could monitor photodynamic therapy or evaluate cranial injuries. While medical applications of the invention have been described, the embodiments of the invention are applicable to any circumstance where image contrast is needed for fluids comprising endogenous metallic compositions.
- MM-ODT, OCT, and ultrasound can be used in various combinations to detect vascular occurrences.
Abstract
A novel contrast mechanism for imaging blood flow using magneto-motive optical Doppler tomography (MM-ODT), Optical Coherence Tomography, and Ultrasound. MM-ODT, OCT, and ultrasound combined with an externally applied temporally oscillating high-strength magnetic field detects erythrocytes moving according to the field gradient.
Description
Attorney Docket: 6100-043 Customer No. 29335
TITLE
[0001] HEMOGLOBIN CONTRAST IN MAGNETO-MOTIVE OPTICAL
DOPPLER TOMOGRAPHY, OPTICAL COHERENCE TOMOGRAPHY, AND ULTRASOUND IMAGING METHODS AND APPARATUS
FUNDING
[0002] This invention was supported by funds from the National Institutes of Health
(AR47551, EB002495 and EB002021) and the Texas Advanced Technology Program. The U.S. Government may have certain rights in the invention.
BACKGROUND
[0003] The present invention relates in general to the art of medical diagnostic imaging and in particular to imaging blood flow using magneto-motive optical Doppler tomography (MM-ODT), Optical Coherence Tomography, or Ultrasound, which combines an externally applied temporally oscillating high-strength magnetic field with ODT, OCT, or Ultrasound to detect erythrocytes moving according to the field gradient
[0004] The accurate determination of location and flow velocity of moving particles in highly scattering media, such as blood flow, is important for medical diagnostics. While the measurements of blood flow in the coronary arteries is an important aspect in diagnosing coronary artery diseases. Numerous non-invasive approaches have been developed using techniques such as Doppler ultrasound, conventional angiography, laser Doppler flowmetry and magnetic resonance angiography.
[0005] One common sensing technique involves the use of ultrasound. Using this technique, ultrasound is directed into the body of the patient and tiny particles such as red blood cells, which are suspended in the blood plasma, scatter the ultrasonic energy back towards the receiver or transducer. The transducer then converts the back-scattered ultrasonic energy into an electrical signal that is processed in some known manner to determine the presence of a flow and an estimate of the flow velocity.
[0006] Magnetic resonance imaging (MRI) is based on an imaging technique for magnetically exciting nuclear spins in a subject positioned in a static magnetic field by applying a radio-frequency (RF) signal of the Larmor frequency, and reconstructing an image
9206 2
Attorney Docket: 6100-043 Customer No. 29335 using MR signals induced by the excitation. MRI is widely applied in clinical medicine because of its capability of clearly depicting the slightest tissue of human brain in vivo.
[0007] Magnetic resonance angiography (MRA) provides detailed angiographic images of the body in a non-invasive manner. In conventional MRA, which does not use contrast agents, magnetic resonance signal from flowing blood is optimized, while signal from stationary blood or tissue structures is suppressed. In contrast-enhanced MRA, a contrast agent is injected into the blood stream to achieve contrast between flowing blood and stationary tissue.
[0008] The commonly known echo planar imaging (EPI) is a rapid MRI technique, which is used to produce tomographic images at high acquisition rates, typically several images per second. Functional magnetic resonance imaging (fMRI) has been found useful in perfusion and/or diffusion studies and in dynamic-contrast studies, etc. However, images obtained in EPI experiments tend to be vulnerable to an artifact known as "ghosting" or "ghost images."
[0009] Optical coherence tomography (OCT) is a technology that allows for noninvasive, cross-sectional optical imaging of biological media with high spatial resolution and high sensitivity. OCT is an extension of low-coherence or white-light interferometry, in which a low temporal coherence light source is utilized to obtain precise localization of reflections internal to a probed structure along an optic axis. This technique is extended to enable scanning of the probe beam in the direction perpendicular to the optic axis, building up a two-dimensional reflectivity data set, used to create a cross-sectional gray-scale or false- color image of internal tissue backscatter.
[0010] OCT uses the short temporal coherence properties of broadband light to extract structural information from heterogeneous samples such as biologic tissue. OCT has been applied to imaging of biological tissue in vitro and in vivo. Systems and methods for substantially increasing the resolution of OCT and for increasing the information content of OCT images through coherent signal processing of the OCT interferogram data have been developed to provide cellular resolution (i.e., in the order of 5 micrometers). During the past decade, numerous advancements in OCT have been reported including real-time imaging speeds.
Attorney Docket: 6100-043 Customer No. 29335
[0011] In diagnostic procedures utilizing OCT, it would also be desirable to monitor the flow of blood and/or other fluids, for example, to detect peripheral blood perfusion, to measure patency in small vessels, and to evaluate tissue necrosis. Another significant application would be in retinal perfusion analysis. Accordingly, it would be advantageous to combine Doppler flow monitoring with the above micron-scale resolution OCT imaging in tissue.
[0012] Conventional OCT imaging primarily utilizes a single backscattering feature to display intensity images. Functional OCT techniques process the backscattered light to provide additional information on birefringence, and flow properties. (See for example, Kemp NJ, Park J, Zaatar HN, Rylander HG, Milner TE, High-sensitivity determination of birefringence in turbid media with enhanced polarization-sensitive optical coherence tomography, Journal of the Optical Society of America A: Optics Image Science and Vision 2005, 22(3):552-560; Dave DP, Akkin T, Milner TE, Polarization-maintaining fiber-based optical low-coherence reflectometer for characterization and ranging of birefringence, Optics Letters 2003, 28(19):1775-1777; Rylander CG, Dave DP, Akkin T, Milner TE, Diller KR, Welch M, Quantitative phase-contrast imaging of cells with phase-sensitive optical coherence microscopy, Optics Letters 2004, 29(13):1509-1511; de Boer JF, Milner TE, Ducros MG, Srinivas SM, Nelson JS, Polarization-sensitive optical coherence tomography, Handbook of Optical Coherence Tomography, New York: Marcel Dekker, Inc., 2002, pp 237-274.)
[0013] Since the ability to characterize fluid flow velocity using OCT was demonstrated by Wang et ah, several phase resolved, real-time optical Doppler tomography (ODT) approaches have been reported. (See for example, Chen ZP, Milner TE, Dave D, Nelson JS, Optical Doppler tomographic imaging of fluid flow velocity in highly scattering media, Optics Letters 1997, 22(l):64-66; Wang XJ, Milner TE, Nelson JS.
[0014] Optical Doppler tomography (ODT) combines Doppler velocimetry with optical coherence tomography (OCT) for noninvasive location and measurement of particle flow velocity in highly scattering media with micrometer-scale spatial resolution. The principle employed in ODT is very similar to that used in radar, sonar and medical ultrasound. ODT uses a low coherence or broadband light source and optical interferometer to obtain high spatial resolution gating with a high speed scanning device such as a conventional rapid scanning optical delay line (RSOD) to perform fast ranging of microstructure and particle motion detection in biological tissues or other turbid media.
9206 2
Attorney Docket: 6100-043 Customer No. 29335
[0015] To detect the Doppler frequency shift signal induced by the moving particles, several algorithms and hardware schemes have been developed for ODT. The most straightforward method to determine the frequency shift involves the use of a short time fast Fourier transform (STFFT). However, the sensitivity of this method is mainly dependent on the FFT time window, which limits axial scanning speed and spatial resolution when measuring slowly moving blood flow in small vessels that requires high velocity sensitivity. However, a phase-resolved technique can decouple the Doppler sensitivity and spatial resolution while maintaining high axial scanning speed.
[0016] In ODT, the Doppler frequency shift is proportional to the cosine of the angle between the probe beam and the scatterer's flow direction. When the two directions are perpendicular, the Doppler shift is zero. Because a priori knowledge of the Doppler angle is not available, and conventional intensity OCT imaging provides a low contrast image of microvasculature structure, detecting small vessels with slow flow rates is difficult. However, the Doppler angle can be estimated by combining Doppler shift and Doppler bandwidth measurements. (See for example, Piao DQ, Zhu Q, Quantifying Doppler Angle and Mapping
Flow Velocity by a Combination of Doppler-shift and Doppler-bandwidth Measurements in
Optical Doppler Tomography, Applied Optics, 2003, 42(25): 5158-5166, and US Patent No.
5,991,697 describe a method and apparatus for Optical Doppler Tomographic imaging of a fluid flow in a highly scattering medium comprising the steps of scanning a fluid flow sample with an optical source of at least partially coherent radiation through an interferometer, which is incorporated herein by reference).
[0017] The ability to locate precisely the microvasculature is important for diagnostics and treatments requiring characterization of blood flow. Recently, several efforts to increase blood flow contrast mechanisms have been reported including protein microspheres incorporating nanoparticles into their shells, plasmon-resonant gold nanoshells, and use of magnetically susceptible micrometer sized particles with an externally applied magnetic field. (See for example, Lee TM, Oldenburg AL, Sitafalwalla S, Marks DL, Luo W, Toublan FJJ, Suslick KS, Boppart SA, Engineered microsphere contrast agents for optical coherence tomography, Optics Letters, 2003, 28(17): 1546-1548; Loo C, Lin A, Hirsch L5 Lee MH, Barton J, Halas N, West J, Drezek R. Nanoshell-enabled photonics-based imaging and therapy of cancer. Technology in Cancer Research & Treatment, 2004; 3(1): 33-40; and
9206 2
Attorney Docket: 6100-043 Customer No 29335
Oldenburg AL, Gunther JR, Boppart SA, Imaging magnetically labeled cells with magnetomotive optical coherence tomography, Optics Letters, 2005, 30(7): 747-749.)
[0018] Wang, et al., "Characterization of Fluid Flow Velocity by Optical Doppler
Tomography," Optics Letters, Vol. 20, No. 11, Jun. 1, 1995, describes an Optical Doppler Tomography system and method which uses optical low coherence reflectrometry in combination with the Doppler effect to measure axial profiles of fluid flow velocity in a sample. A disadvantage of the Wang system is that it does not provide a method to determine direction of flow within the sample and also does not provide a method for generating a two- dimensional color image of the sample indicating the flow velocity and directions within the image.
[0019] The use of an externally applied field to move magnetically susceptible particles in tissue has been termed magneto-motive OCT (MM-OCT). Functional magnetic resonance imaging (fMRI) detects deoxyhemoglobin which is a paramagnetic molecule. However, the paramagnetic susceptibility of human tissue is very low compared to other biocompatible agents such as ferumoxides (nanometer sized iron oxide particles). Therefore, it was believed that, other than differentiating relaxation times (T2) between oxygenated and deoxygenated blood, the magnetic field strength required to produce a retarding force on blood flow was well above that of current imaging fields. (See also, for example, Schenck JF., Physical interactions of static magnetic fields with living tissues, Progress in Biophysics and Molecular Biology 2005, 87(2-3): 185-204; and Taylor DS, Coryell, CD., Magnetic susceptibility of iron in hemoglobin. J. Am. Chem. Soc. 1938, 60:1177-1181.)
[0020] The present invention is a novel extension of MM-OCT to image hemoglobin in blood erythrocytes. In one aspect of the invention, the approach requires no exogenous contrast agent to detect blood flow and location.
SUMMARY OF THE INVENTION
[0021] The invention is a method and apparatus of imaging a blood flow using
Doppler optical coherence tomography, which comprises an externally applied temporally oscillating high-strength magnetic field with a Doppler optical tomography system to detect hemoglobin moving according to the magnetic field gradient.
9206 2
Attorney Docket: 6100-043 Customer No. 29335
[0022] In another embodiment, the invention comprises a solenoid cone-shaped ferrite core with an extensively increased magnetic field strength at the tip of the core, hi a further embodiment, the invention comprises focusing the magnetic force on targeted samples.
[0023] Another embodiment is an apparatus for imaging blood flow, comprising: (1) a magnetomotive optical Doppler tomography (MM-ODT) imaging system; (2) an oscillating magnetic field applied to the moving blood flow; (3) a light source (4) a rapid scanning optical delay line in a reference arm and aligned such that no phase modulation is generated when the group phase delay is scanned at 4 kHz; and (5) a hardware in-phase and a quadrature demodulator with high bandpass filters.
[0024] Another aspect of the invention is a method for imaging a blood flow, comprising applying a magnetic field to the blood flow, wherein the blood flow comprising a plurality of hemoglobin molecules and wherein the magnetic field interacts with the hemoglobin molecules to cause a change in the blood flow; and detecting the blood flow by detecting the change in the blood flow caused by the interaction with the hemoglobin molecules, wherein the change is detected using a magnetomotive optical Doppler tomography imaging system.
[0025] Another embodiment of the method for imaging a blood flow comprises temporally oscillating the magnetic field.
[0026] Another embodiment of the invention is an apparatus for imaging a blood flow comprising a magnetic field generator for applying a magnetic field to the blood flow, wherein the blood flow comprises hemoglobin molecules; and an ultrasound detection system for detecting the blood flow while it is in the presence of the magnetic field.
BRIEF DESCRIPTION OF THE DRAWINGS
[0027] FIG. 1 is a schematic diagram of the MM-ODT system.
[0028] FIG. 2 is a schematic diagram of the probe beam, flow sample and solenoid coil.
[0029] FIG. 3a and 3* are OCT and ODT images of a stationary turbid solution without an external magnetic field, respectively. FIG. 3c and FIG. 3d are OCT and ODT
Attorney Docket: 6100-043 Customer No. 29335 images with a 50 Hz magnetic field, respectively. The white bar represents 200 μm, accordingly.
[0030] FIGS. 4a-4d are M-mode ODT images of the diluted deoxygenated blood flow (18% hematocrit) without and with an external magnetic field. FIG. 4a and FIG. 46 are ODT images of 5 mm/s blood flow without and with a 5 Hz magnetic field, respectively. FIG. 4c and FIG. 4d are ODT images of 30 mm/s blood flow without and with a 50 Hz magnetic field, respectively. The Black bar indicates 200 μm, accordingly.
[0031] FIG. Sa and FIG. Sb axe the Doppler frequency shift profiles without an external magnetic field, and with a 50 Hz magnetic field, respectively.
[0032] FIG. 6 is a block diagram illustrating an exemplary phase sensitive OCT system.
[0033] FIG. 7 is a graph of the optical path length change in the DP-OCT system.
[0034] FIG. 8 is a schematic diagram of Magneto Motive Ultrasound for Blood.
[0035] FIG. 9 is an ultrasound image of blood exposed to an oscillating magnetic field.
[0036] FIG. 10 is a schematic diagram of Ultrasonography and OCT coupled with magnetic field generator.
DETAILED DESCRIPTION OF THE INVENTION
[0037] Magneto-Motive Optical Doppler Tomography (MM-ODT) for improved
Doppler imaging of blood flow using an external oscillating magnetic field is described below. By introducing mechanical movement of red blood cells (RBCs) during blood flow by a temporally oscillating high-strength magnetic field, MM-ODT allows imaging of blood flow and velocity. The controlled and increased Doppler frequency in MM-ODT provides an investigational tool to study in vivo blood transport, as shown in the article Hemoglobin Contrast in Magnetomotive Optical Doppler Tomography, Opt. Lett. 31, 778-780 (2006), herein incorporated by reference.
[0038] The microstructure of the blood flow and flow velocity information are all encoded in the interferogram of a Doppler OCT system. It should be readily apparent to those skilled in the optical arts, that different OCT systems and different OCT information
9206_2 7
Attorney Docket: 6100-043 Customer No. 29335 can be used to determine the Doppler frequency shift. It is not intended to suggest any limitation as to the scope or functionality with different OCT architectures or optical information used with an oscillating magnetic field, such as time domain Doppler OCT and spectral domain Doppler OCT. Time domain OCT requires a mirror in the reference arm scanning at a constant velocity, while spectral domain OCT includes swept source OCT and Fourier domain OCT. An example of time-domain Doppler OCT is provided below.
[0039] A schematic of the MM-ODT apparatus 10 is shown in FIG. 1. The OCT light source 12 comprises a super luminescent diode, which is used as the low coherence light source, hi one embodiment, light source 12 is centered at 1.3 μm with a bandwidth of 90 run. Light from source 12 is coupled into a single-mode optical fiber based interferometer 20 by circulator 14, where the interferometer 20 can provide 1 mW of optical power at the sample 60. Light is split into a reference arm 22 and sample arm 24 by a 2 x 2 splitter 16. A rapid- scanning optical delay (RSOD) line 26 is coupled to the reference arm 22. In one embodiment, the rapid-scanning optical delay line 26 is aligned such that no phase modulation is generated when the group phase delay is scanned at ~ 4 kHz. In the sample arm 24, a collimated beam 36 is redirected to sample 60 by two galvanometers 30 that permit three-dimensional scanning. In one embodiment, the galvanometers can be an X-scanner and a Y-scanner. The sample 60 can be internal or external to the body, where the probe beam is focused by an objective lens 32. In one embodiment, the objective lens 32 yields a 10-μm spot at the focal point. Phase modulation can be generated using an electro-optical waveguide phase modulator, which can produce a carrier frequency (~1 MHz). And the magnetic field generator 100 is in proximity to the sample 60.
[0040] A dual-balanced photodetector 34 is coupled to the 2 x 2 splitter 16 and the circulator 14. The photodetector 34 of a 80 MHz bandwidth reduces the light source noise from the OCT interference signal. A hardware in-phase and quadrature demodulator 40 with high/bandpass filters 42 and low/bandpass filters 44 improves imaging speed. Doppler information was calculated with the Kasai autocorrelation velocity estimator. Labview software 50 (National Instruments, Austin, TX) is coupled to the MM-ODT system with a dual processor based multitasking scheme. The maximum frame rate of the MM-ODT system 10 was 16 frames per second for a 400 x 512 pixel sized image. The Doppler frequency shift can be determined with the use of a short time fast Fourier transform (STFFT). Alternatively, a phase-resolved technique can determine the Doppler frequency shift to decouple the
9206 2
Attorney Docket: 6100-043 Customer No. 29335
Doppler sensitivity and spatial resolution while maintaining high axial scanning speed. Alternatively, differential phase optical coherence tomography (OCT) or spectral domain phase-sensitive OCT can be used to determine the Doppler frequency shift, as readily determined by one skilled in the optical arts.
[0041] FIG. 2 shows an example of the magnetic field generator 80 with a capillary tube 90. The magnetic field generator 80 includes a solenoid coil 82 (Ledex 4EF) with a cone-shaped ferrite core 84 at the center and driven by a current amplifier 86 supplying up to 960 W of power. The magnetic field generator 80 can be placed underneath the sample 60 during MM-ODT imaging. The combination of the ferrite core 84 and solenoid coil 82 using a high power operation dramatically increases the magnetic field strength (Bmax=0.14 Tesla) at the tip of the core 82 and also focuses the magnetic force on the targeted samples 60. The sinusoidal current can vary the magnetic force applied to the capillary tube 90 in order to introduce movement of magnetic fluids, which include red blood cells that contain hemoglobin. In one embodiment, the probe beam is oriented parallel to the gradient of the magnetic field's strength.
[0042] The material parameter characterizing magnetic materials, including biological tissue, is the magnetic volume susceptibility, χ. Magnetic volume susceptibility is dimensionless in SI units and is defined by the equation M — χH, where M is the magnetization at the point in question and H is the local density of the magnetic field strength. Hemoglobin's high iron content, due to four Fe atoms in each hemoglobin molecule, and the large concentration of hemoglobin in human red blood cells give Hemoglobin magnetomotive effects in biological tissue. The magnetic volume susceptibility of the hemoglobin molecule consists of a paramagnetic component due to the electron spins of the four iron atoms. The paramagnetic susceptibility is given by the Curie Law,
where μ0 is permeability of free space and has the value 4πxlO7 H/m, Np is the volume density of paramagnetic iron atoms in hemoglobin, N^= 4.97x 1025 iron atoms/m^ μe/f is the effective number of Bohr magnetons per atom reported as 5.35, and the Bohr magneton, μB :=9.274xl0"24 J/T, and Boltzmann's constant, k=1.38x 10~23 J/K, and T is the absolute
Attorney Docket: 6100-043 Customer No. 29335 temperature (K). The calculated susceptibility of a RBC is about HxIO"6 assuming a 90% concentration of hemoglobin per RBC. The calculated susceptibility of a RBC is dependent on the oxygenation of the hemoglobin. The calculations can be adjusted accordingly, depending on the oxygenation of the RBC, which can be measured by known techniques.
[0043] A RBC placed in a magnetic field gradient experiences forces and torques that tend to position and align it with respect to the field's direction. The magnetic force, in the direction of the probing light z, is given by
where V is the particle volume, B is the magnitude of the magnetic flux density, and Aχ is the difference between the susceptibility of the particle and the surrounding medium. The displacement [z(t)] of an RBC driven by a time varying magnetic flux density can be included in the analytic OCT fringe expression, If,
where IR and Is are the back scattered intensities from the reference and sample arms, respectively, fo is the fringe carrier frequency, X0 is the center wavelength of the light source, and z(t) is the RBC displacement. Integration of all forces (magnetic, elastic, and viscous) on the RBC gives displacement, z(t)= Acos(4πfmt) where A is a constant in units of length and fm is the modulation frequency of the magnetic flux density. In free space, the displacement, z(t), is dominated a constant acceleration which can be, however, ignored in confined models (i.e. blood vessel or capillary tube) with assumptions that, first, the probing area is much smaller than magnetic field area, and that secondly probing time starts after steady states of inner pressures. Expansion of the right-hand side of Eq. (3) using Bessel functions gives
where Jk(m) is the Bessel function of the first kind of order k for argument m which is 4πA/λ0. The amplitude of the k* sideband is proportional to jk(m). In coherent detection, the fraction of
9206_2 10
Attorney Docket: 6100-043 Customer No. 29335 optical power transferred into each of the first order sidebands is (Jj (m)) , and the fraction of optical power that remains in the carrier is (J0(m))2.
Example 1
MM-ODT imaging of the Doppler shift of hemoglobin by applying an oscillating magnetic field to a moving blood sample
[0044] M-mode OCT/ODT images of a capillary glass tube filled with a stationary turbid solution with and without an external magnetic field as a control sample were recorded, as shown in FIGS. 3-4. A 750 μm-inner diameter glass capillary tube 200 was placed perpendicularly to the probing beam 70, as shown in FIG. 2. Fluids used for flow studies were injected through the tube at a constant flow rate controlled by a dual-syringe pump (Harvard Apparatus 11 Plus, Holliston, MA) with ±0.5% flow rate accuracy. The turbid solution was a mixture of deionized water and 0.5-gm latex microspheres (μs=5mm"1). The magnetic flux density and its frequency were approximately 0.14 T and 50 Hz, respectively. M-mode OCT/ODT images were acquired for 100 ms per frame. FIGS. 3a and 3b show M-mode OCT and ODT images without any external magnetic field, respectively. The ODT image in FIG. 3b contains small random phase fluctuations due to ambient vibration through the optical path. FIGS. 3c and 3d show M-mode OCT and ODT images with a 50 Hz external magnetic field, respectively. No distinguishable Doppler shift could be observed in the ODT image FIG. 3d indicating no interaction between the external magnetic field and the moving microspheres.
[0045] Deoxygenated blood was extracted from the vein of a human male's left arm, and diluted with saline. During preparation, blood was not exposed directly to air so as to remain deoxygenated. To simulate flow, blood was injected through the capillary tube 200 by a syringe pump at a relatively constant flow rate. As FIG. 4 shows, the oscillating Doppler frequency shift, resulting from RBC movement, could be observed at two different flow rates (5 and 30 mm/sλ Because the flow direction was nearly perpendicular to the probing beam no significant Doppler frequency shift was distinguishable at the 5 mm/s flow rate FIG. 4a without any external magnetic field. In the case of the high blood flow rate of 30 mm/s, as shown at FIG. 4c, the Doppler frequency shift caused by the blood flow could be observed. And for maximum contrast enhancement, the probe beam can be directed parallel to the gradient of the magnetic field's strength. However, application of a 50 Hz magnetic
9206_2 11
Attorney Docket: 6100-043 Customer No. 29335 field increased the Doppler contrast of blood at both the slow and fast flow rates as shown at FIGS. 4b and 4d. The high flow rate of 30 mm/s gives a higher contrast image than the low flow rate image, but the Doppler frequency shift of the former as a function of depth is less homogeneous than the latter, which is indicative of perturbation by blood flow. The same blood was diluted to 5 % hematocrit (HCT), but no RBC movement could be observed below 8 % HCT.
[0046] Doppler frequency shift profiles were calculated from the ODT images by averaging 20 lines at a selected depth indicated by horizontal arrows, as shown in FIGS. 4a and 4b. FIG. 5a indicates no significant Doppler frequency shift over a 100 ms time period, whereas FIG. 5b displays ±200 to 300 Hz Doppler frequency shifts oscillating 20 times over 100 ms (200 Hz).
[0047] The Doppler frequency shift indicates that RBCs physically move into and away from the incident light while passing through the external magnetic field depending on whether their magnetic properties are paramagnetic or diamagnetic, as shown in FIG. 4, and that the 200 Hz oscillation of the Doppler frequency shift correlates with the 50 Hz magnetic field, as shown in FIG. 5. Frequencies 4 times higher than that of the external magnetic field (fan) can be observed. According to Eq. 2, the frequency of the force on paramagnetic targets was twice that of the magnetic flux density; therefore, a 50 Hz B field displaces the targets at 100 Hz. Although the fringe signal (Eq. 4) contains harmonics at frequency (2^n), the modulation frequency, fm, was set so that the second-order sideband (4/Jn) was dominant as shown in FIG. 5b. The particle motion can not be described as a pure sinusoidal function even if the modulated magnetic field is sinusoidal, due to the numerous forces that contribute to the motion within the field such as gravity, concentration gradient, and colloidal dispersion.
[0048] The invention is a new investigational tool to study in vivo blood transport and the first implementation of MM-ODT for improved Doppler imaging of blood flow using an external oscillating magnetic field introducing a mechanical movement of RBCs during blood flow by a temporally oscillating high-strength magnetic field. MM-ODT to allow imaging of tissue function in a manner similar to functional magnetic resonance images (f- MRI) of deoxygenated blood hi organs, when the sample arm of the MM-ODT system is coupled to a probe (not shown). Such probes are generally known in the arts, such as endoscopic probes, catheter probes, and the like.
9206_2 12
Attorney Docket: 6100-043 Customer No. 29335
[0049] Alternatively, the MM-ODT can be used for Port- Wine vessel mapping and
Skin Cancer vessel mapping. The MM-ODT can be used to detect blood vessel location and size for cancer and port-wine stains, since these conditions are characterized by blood vessel growth and increases in hemoglobin content. Accordingly, other blood vessel detection for tissue abnormality identification can be envisioned with this invention. Generally, the MM- ODT can be used wherever blood flow detection is necessary in operations, chemotherapy, hemodialysis, and the like.
[0050] A spectral domain phase sensitive OCT system 100 can be used to image the blood flow and determine velocity information, when an oscillating magnetic field is applied to the blood flow, as shown in FIG. 6. Spectral domain phase sensitive OCT generally uses a broadband light source in a general interferometric setup, where the mirror in the reference typically does not move or does not require a rapid scanning delay line. A spectral interferogram is detected, in which each A-scan is entirely encoded. A Fourier transformation of the A-scan and velocity distribution of the blood flow can be extracted by known methods. In spectral domain Doppler OCT, the blood flow profile is the envelope of a calculated A-scan through Fourier transformation on the spectral interferogram. Velocity information is encoded in the phase of the interferogram. The Doppler velocity can be extracted by measuring the phase shift between two successive calculated A-scans in spectral domain Doppler OCT. [0051] This OCT system 100 is only an example of one OCT imaging modality which can be used to image blood flow with a temporally oscillating magnetic field, and is not intended to suggest any limitation on the scope of OCT architectures applicable to the invention. Generally, the OCT system 100 includes a general-purpose computing device in the form of a computer 101 and includes a magnet control 114 and a magnet 116.
[0052] Light energy is generated by a light source 117. The light source 117 can be a broadband laser light source coupled into optical fiber emitting light energy over a broad range of optical frequencies. The wavelength range can be from about 400 nanometers to about 1400 nanometers. Longer wavelengths (>600 run) can be used for deeper scanning. Preferably, the light source emits light having a wavelength near the infrared spectrum to identify hemoglobin for OCT imaging, which places hemoglobin in motion and increases optical scattering of the hemoglobin. The light energy can be emitted over a multiplicity of optical wavelengths, frequencies, and pulse durations to achieve OCT imaging. As used
9206_2 13
Attorney Docket: 6100-043 Customer No. 29335 herein, optical fiber can refer to glass or plastic wire or fiber. Optical fiber is indicated on FIG. 6 as lines connecting the various blocks of the figures. Where light energy is described as "passing," "traveling," "returning," "directed," or similar movement, such movement can be via optical fiber.
[0053] A fraction of the generated light energy passes from the light source 117 into an optical spectrum analyzer 118. The optical spectrum analyzer 118 measures optical frequency as the light energy is emitted from the light source 117 as a function of time. The' optical spectrum analyzer 118 samples a portion of the light emitted by the light source 117. The optical spectrum analyzer 118 monitors the power spectral density of light entering the splitter 119. The remaining fraction of light energy from the light source 117 passes into a splitter 119. The splitter 119 can be a device with four ports, with Port 1 allowing light energy to enter the splitter 119. Ports 2 and 3 allow light energy to leave and re-enter the splitter 119 to the reference reflector 120 and OCT probe 122, respectively. Port 4 allows light energy to leave the splitter 119 to coupling lens 124. The splitter 119 couples the light into Port 1. The splitter 119 divides the light according to a pre-determined split ratio selected by a user. For example, the split ratio can be 50/50 wherein half of the light energy entering the splitter 119 at Port 1 exits the splitter 119 through Port 2 and half exits the splitter 119 through Port 3. In another example, the split ratio can be 60/40 wherein 60% of the light energy passes through Port 2 and 40% of the light energy passes through Port 3.
[0054] A fraction of the light energy (determined by the split ratio) that exits the splitter 119 through Port 2 travels to a reference reflector surface 120. The light energy is reflected from the reference reflector surface 120 back to the splitter 119 into Port 2. The reference reflector 120 can be a planar metallic mirror or a multilayer dielectric reflector with a specified spectral amplitude/phase reflectivity. The remaining fraction of light that entered splitter 119 through Port 1 exits splitter 119 through Port 3 and enters an OCT probe 122. The OCT probe 122 can be a turbine-type catheter as described in Patent Cooperation Treaty application PCT/US04/12773 filed 4/23/2004 which claims priority to U.S. provisional application 60/466,215 filed 4/28/2003, each herein incorporated by reference for the methods, apparatuses and systems taught therein. The OCT probe 122 can be located within a subject 121 to allow light reflection off of subject's 121 blood flow.
9206 2 14
Attorney Docket: 6100-043 Customer No. 29335
[0055] The light energy that entered OCT probe 122 is reflected off of the blood flow of subject 121. once an oscillating magnetic field has been temporally applied by magnet 116. The reflected light energy passes back through the OCT probe 122 into the splitter 119 via Port 3. The reflected light energy that is returned into Port 2 and Port 3 of the splitter 119 recombines and interferes according to a split ratio. The light recombines either constructively or destructively, depending on the difference of pathlengths. A series of constructive and destructive combinations of reflected light create an interferogram (a plot of detector response as a function of optical path length difference). Each reflecting layer from the subject 121 and the blood flow will generate an interferogram. The splitter 119 can recombine light energy that is returned through Port 2 and Port 3 so that the light energies interfere. The light energy is recombined in the reverse of the split ratio. For example, if a 60/40 split ratio, only 40% of the light energy returned through Port 2 and 60% of the light energy returned through Port 3 would be recombined. The recombined reflected light energy is directed out Port 4 of the splitter 119 into a coupling lens 137. The coupling lens 137 receives light from the output of the splitter 119 and sets the beam etendue (beam diameter and divergence) to match that of the optical spectrometer 138. The coupling lens 137 couples the light into an optical spectrometer 138. The optical spectrometer 138 can divide the recombined reflected light energy light into different optical frequencies and direct them to different points in space which are detected by a line scan camera 139. The line scan camera 139 performs light to electrical transduction resulting in digital light signal data 108. The digital light signal data 108 is transferred into the computer 101 via the OCT input interface 111. Interface between the line scan camera 139 and computer 101 can be a Universal Serial Bus (USB)3 or the like. The digital light signal data 108 can be stored in the mass storage device 104 or system memory 112 and utilized by the image construction software 106 and the Labview image construction software 107.
[0056] The image construction software 106 can generate an image of the blood flow of subject 121 from the light signal data 108, by receiving light signal data 108 generating amplitude and phase data. The amplitude and phase data (optical path length difference (cτ) or optical time-delay (τ)) can be separated into discrete channels and a plot of intensity vs. depth (or amplitude vs. depth) can be generated for each channel. Such plot is known as an A-scan, where the composition of all the A-scans can comprise one image. And movement image construction software 107 generates an image of the movement of the hemoglobin from the light signal data 108. The movement image construction software 107 receives light
9206 2 15
Attorney Docket: 6100-043 Customer No. 29335 signal data 108 for at least two successive sweeps of the light source 117 or the light source performs a Fourier transform on the light signal data 108 generating amplitude and phase data.
[0057] Optionally, additional information can be extracted from the light signal data to generate additional images. The light signal data can be further processed to generate a Stokes parameter polarimetric image when used in conjunction with polarization detectors and polarizing lenses to extract polarization data from the light signal 108, as readily known to one skilled in the art of optical coherence tomography. The differential phase OCT image data is shown in FIG. 7, indicating the optical path difference with the magnetic field. Filters can be added to reduce the noise of the signal generated by the magnetic field.
[0058] Alternatively, the phase sensitive OCT system 100 can be configured for swept source OCT, which is a different type of spectral domain OCT. In swept source OCT, a tunable laser source replaces the broadband laser light source 117. The scanning rate can be at wavelengths of 800nm~1500nm. Also, the reference reflector surface 120 is in-line with sample path 120. The optical spectrometer 125 and line scan camera 126 are replaced with a general photodetector.
[0059] In another embodiment, an enhanced detection of cancer with ultrasound imaging 200 is provided. Ultrasonography is the ultrasound-based diagnostic imaging technique used to visualize muscles and internal organs, their size, structures and any pathological lesions. "Ultrasound" applies to all acoustic energy with a frequency above human hearing (20,000 Hertz or 20 kilohertz). Typical diagnostic sonography scanners operate in the frequency range of 2 to 13 megahertz, hundreds of times greater than this limit. The choice of frequency is a trade-off between the image spatial resolution and penetration depth into the patient, with lower frequencies giving less resolution and greater imaging depth. Doppler ultrasonography uses the Doppler Effect to assess whether blood is moving towards or away from a probe, and its relative velocity. By calculating the frequency shift (υD) of a particular sample volume, for example a jet of blood flow over a heart valve, its speed and direction can be determined and visualized. Ultrasonagraphy and Doppler Ultrasonagraphy can best be understood by S. A. Kana Introduction to physics in modern medicine, Taylor & Francis, (2003). The basic physics of the Doppler effect involving
9206 2 16
Attorney Docket: 6100-043 Customer No. 29335 acoustic and electromagnetic waves of OCT is similar and many of the signal processing techniques (hardware and software) used to estimate the Doppler shift is analogous.
[0060] In one embodiment of the invention, an ultrasound probe 212 is coupled with the magnetic field generator 100, as shown in FIG. 8. Ultrasound 210 is directed into the body of the patient by known techniques and the moving red blood cells backscatter the ultrasonic energy back towards the transducer of the ultrasound. The oscillating magnetic field generated by the magnetic field generator 100 increases the contrast of the ultrasonic energy 210 received from the red blood cells. The transducer then converts the back- scattered ultrasonic energy 210 into an electrical signal that is processed in some known manner to determine an estimate of the flow. An enhanced ultrasound image is produced, as displayed in FIG. 9.
[0061] hi one example, a rectal ultrasound probe is coupled with a magnet to evaluate the prostate gland for cancer. Currently, ultrasound is used for prostate cancer screening; however, the approach provides poor sensitivity and specificity. Yet, all cancers are known in the art to be highly vascular, so then the application of the magnetic field by generator 100 coupled with ultrasound enhances the contrast available from the endogenous RBCs in the prostate for cancer detection at an earlier stage. It is generally known in the art that cancers have enhanced metabolic properties compared to normal tissues, so then cancerous cells have higher oxygen content from hemoglobin and a greater concentration of deoxygenated hemoglobin compared to normal tissues. An exemplary ultrasound image for prostate cancer screening is shown in FIG. 9.
[0062] MM-ODT, OCT, and ultrasound coupled with an oscillating magnetic field all can be used in clinical management of patients who need microvasculature monitoring, as shown in FIG. 10. The images could monitor and determine tissue profusion and viability before, during and after operation procedures. For example, the images could be used to detect oxygenated and deoxygenated blood supply, detect blood vessel location, sizes of cancer, vascular tissue abnormality identification, and port-wine stain. Alternatively, the images could monitor photodynamic therapy or evaluate cranial injuries. While medical applications of the invention have been described, the embodiments of the invention are applicable to any circumstance where image contrast is needed for fluids comprising endogenous metallic compositions. MM-ODT, OCT, and ultrasound can be used in various combinations to detect vascular occurrences.
9206_2 17
Attorney Docket: 6100-043 Customer No. 29335
[0063] While the invention has been described in connection with various embodiments, it will be understood that the invention is capable of further modifications. This application is intended to cover any variations, uses or adaptations of the invention following, in general, the principles of the invention, and including such departures from the present disclosure as, within the known and customary practice within the art to which the invention pertains.
9206 2 18
Claims
1. A system for imaging blood vessels, comprising: a magnet for applying a magnetic field to the vessels to be imaged, wherein the applied magnetic field interacts with hemoglobin within the vasculature; and an optical coherence tomography apparatus for detecting magnetic movement of the hemoglobin within the vasculature.
2. The system of Claim 1, wherein the magnetic field is oscillating.
3. The system of Claim 1, wherein the magnet comprises a solenoid coil with a ferrite core.
4. The system of Claim 3, wherein the optical coherence tomography apparatus further comprises a sample arm, and wherein a probe is coupled to the sample arm of the optical coherence tomography apparatus and the magnet.
5. The system of Claim 4, wherein the optical coherence tomography apparatus is a magneto-motive optical Doppler tomography imaging system for detecting the blood flow affected by the magnetic field.
6. The system of Claim 5, wherein the magneto-motive optical Doppler tomography system comprises: a light source generating light energy; an interferometer coupled to the light energy, wherein the interferometer includes a reference and sample light paths coupled to a light splitter; a modulator coupled to the interferometer for modulating the optical path length difference in the reference arm and the sample arm; a scanner coupled to the sample arm for scanning a biological sample; a rapid scanning optical delay line coupled to the reference arm; a photodetector coupled to the interferometer for detecting backscattered radiation received by the interferometer from the scanner to detect interference fringes; and a processor for processing the reflected light energy from the reference arm and a magnetic resonance signal reflected off of the moving blood flow to produce a tomographic image and a tomographic flow velocity image.
9206_2 19 Attorney Docket: 6100-043 Customer No. 29335
7. The system of Claim 6, wherein the magneto-motive optical Doppler tomography system includes the scanner with at least two galvanometers to permit three dimensional scans.
8. The system of Claim 7, wherein the magneto-motive optical Doppler tomography system includes a dual-balanced photodetector.
9. The system of claim 8, wherein the magneto-motive optical Doppler tomography system includes a circulator coupled to the interferometer and the photodetector.
10. The system of claim 9, wherein the magneto-motive optical Doppler tomography system includes the photodetector coupled to a hardware in-phase and a quadrature demodulator having high bandpass filters.
11. A method for imaging a blood flow, comprising: applying a magnetic field to the blood flow, wherein the blood flow comprises a plurality of hemoglobin molecules and wherein the magnetic field interacts with the hemoglobin to cause a change in the blood flow; and detecting the blood flow by detecting the change in the blood flow caused by the interaction with the hemoglobin molecules, wherein the change is detected using a optical coherence tomography system.
12. The method of Claim 11, wherein the applying of the magnetic field comprises temporally oscillating the magnetic field.
13. The method of Claim 12, further comprising coupling the optical coherence tomography system and the magnetic field to a probe.
14. The method of Claim 13, wherein the detecting the blood flow by detecting the change in the blood flow caused by the interaction with the hemoglobin molecules, wherein the change is detected using a magnetomotive optical Doppler tomography imaging system.
9206 2 20 Attorney Docket 6100-043 Customer No. 29335
15. The method of Claim 14, wherein the magneto-motive optical Doppler tomography system comprises the method of providing light energy through an interferometer; phase modulating the light energy in the interferometer at a modulation frequency; continuously scanning a blood flow sample with the light energy through the interferometer, wherein the blood flow sample includes a blood flow therein and a structure in which the blood flow is defined; detecting the magnetic resonance signal reflected off the moving blood sample and the interference fringes of the light energy backscattered from moving blood sample; and data processing Doppler frequency changes of the 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 blood flow in and the structure of said scanned blood flow sample is formed.
16. The method of Claim 15, where phase modulating the light energy includes aligning a rapid scanning optical delay line, and generating no phase modulation when the group phase delay is scanned at 4 kHz.
17. The method of Claim 16, where detecting the interference fringes of light energy backscattered from the moving blood sample includes reducing the light source noise from the interference signal with a dual balanced photodetector.
18. The method of Claim 17, further including improving imaging speed with a hardware in- phase and a quadrature demodulator with at least one high-bandpass filter.
19. The method of Claim 18, where the data processing Doppler frequency changes includes calculating with a Kasai autocorrelation velocity estimator.
20. The method of Claim 19, where the continuously scanning of the blood flow comprises three dimensionally scanning the blood flow with a plurality of galvanometers.
9206_2 21 Attorney Docket: 6100-043 Customer No. 29335
21. The apparatus of Claim 4, wherein the optical coherence tomography system for detecting blood flow is a spectral domain phase sensitive optical coherence tomography system.
22. The apparatus of Claim 21 , wherein the phase sensitive optical coherence tomographic imaging modality is a swept source phase sensitive optical coherence tomography system.
23. The method of Claim 13, wherein the detecting of the change in the blood flow caused by the interaction with the hemoglobin molecules is detected using a spectral domain phase sensitive optical coherence tomography system.
24. The method of Claim 23, wherein the detecting of the change in the blood flow caused by the interaction with the hemoglobin molecules is detected using a swept source phase sensitive optical coherence tomography system.
25. An apparatus for imaging a blood flow comprising: a magnetic field generator for applying a magnetic field to the blood flow, wherein the blood flow comprises hemoglobin molecules; and an ultrasound detection system for detecting the blood flow while it is in the presence of the magnetic field.
26. The apparatus of Claim 25, wherein the ultrasound detection system comprises a probe for transmitting and receiving sound energy to and from a subject and the probe includes the magnetic field generator.
9206 2 22
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2006/040932 WO2008048263A1 (en) | 2006-10-18 | 2006-10-18 | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2006/040932 WO2008048263A1 (en) | 2006-10-18 | 2006-10-18 | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2008048263A1 true WO2008048263A1 (en) | 2008-04-24 |
Family
ID=39314304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2006/040932 WO2008048263A1 (en) | 2006-10-18 | 2006-10-18 | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2008048263A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104545872A (en) * | 2015-01-12 | 2015-04-29 | 南京理工大学 | Method and device for reconstructing three-dimensional micro blood flow distribution on basis of linearly dependent coefficients |
WO2023288073A1 (en) * | 2021-07-15 | 2023-01-19 | The Regents Of The University Of California | Universal photonics tomography |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030004412A1 (en) * | 1999-02-04 | 2003-01-02 | Izatt Joseph A. | Optical imaging device |
US20030220749A1 (en) * | 2002-04-09 | 2003-11-27 | Zhongping Chen | Phase-resolved functional optical coherence tomography: simultaneous imaging of the stokes vectors, structure, blood flow velocity, standard deviation and birefringence in biological samples |
US20050171433A1 (en) * | 2004-01-08 | 2005-08-04 | Boppart Stephen A. | Multi-functional plasmon-resonant contrast agents for optical coherence tomography |
US20060025659A1 (en) * | 2004-08-02 | 2006-02-02 | Masashi Kiguchi | Optical measuring device for substances in vivo |
US20060256343A1 (en) * | 2005-01-20 | 2006-11-16 | Michael Choma | Methods, systems and computer program products for characterizing structures based on interferometric phase data |
-
2006
- 2006-10-18 WO PCT/US2006/040932 patent/WO2008048263A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030004412A1 (en) * | 1999-02-04 | 2003-01-02 | Izatt Joseph A. | Optical imaging device |
US20030220749A1 (en) * | 2002-04-09 | 2003-11-27 | Zhongping Chen | Phase-resolved functional optical coherence tomography: simultaneous imaging of the stokes vectors, structure, blood flow velocity, standard deviation and birefringence in biological samples |
US20050171433A1 (en) * | 2004-01-08 | 2005-08-04 | Boppart Stephen A. | Multi-functional plasmon-resonant contrast agents for optical coherence tomography |
US20060025659A1 (en) * | 2004-08-02 | 2006-02-02 | Masashi Kiguchi | Optical measuring device for substances in vivo |
US20060256343A1 (en) * | 2005-01-20 | 2006-11-16 | Michael Choma | Methods, systems and computer program products for characterizing structures based on interferometric phase data |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104545872A (en) * | 2015-01-12 | 2015-04-29 | 南京理工大学 | Method and device for reconstructing three-dimensional micro blood flow distribution on basis of linearly dependent coefficients |
CN104545872B (en) * | 2015-01-12 | 2017-04-19 | 南京理工大学 | Method and device for reconstructing three-dimensional micro blood flow distribution on basis of linearly dependent coefficients |
WO2023288073A1 (en) * | 2021-07-15 | 2023-01-19 | The Regents Of The University Of California | Universal photonics tomography |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8036732B2 (en) | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus | |
US9532718B2 (en) | Hemoglobin contrast in ultrasound and optical coherence tomography for diagnosing diseased tissue, cancers, and the like | |
US9687153B2 (en) | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus | |
Daly et al. | ‘Go with the flow’: a review of methods and advancements in blood flow imaging | |
US7751057B2 (en) | Magnetomotive optical coherence tomography | |
US6549801B1 (en) | Phase-resolved optical coherence tomography and optical doppler tomography for imaging fluid flow in tissue with fast scanning speed and high velocity sensitivity | |
Hosseinaee et al. | Towards non-contact photoacoustic imaging | |
Yousefi et al. | Eigendecomposition-based clutter filtering technique for optical microangiography | |
Liang et al. | Dynamic optical coherence elastography: a review | |
Liu et al. | Advances in doppler oct | |
Chen et al. | Doppler optical coherence tomography | |
Wijesinghe et al. | In vivo imaging of melanoma-implanted magnetic nanoparticles using contrast-enhanced magneto-motive optical Doppler tomography | |
Bonesi et al. | Study of flow dynamics in complex vessels using Doppler optical coherence tomography | |
Mackle et al. | Fibre optic intravascular measurements of blood flow: A review | |
Meglinski et al. | Diffusing wave spectroscopy: application for blood diagnostics | |
Oldenburg et al. | Magnetic and plasmonic contrast agents in optical coherence tomography | |
WO2008048263A1 (en) | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus | |
Subhash | Biophotonics Modalities for High-Resolution Imaging of Microcirculatory Tissue Beds Using Endogenous Contrast: A Review on Present Scenario and Prospects | |
Chang et al. | High-resolution optical Doppler tomography for in vitro and in vivo fluid flow dynamics | |
US9486140B2 (en) | Coherence gated doppler motion sensor and medical applications | |
Chen et al. | Doppler optical coherence tomography and its application in measurement of cerebral blood flow | |
Gelikonov et al. | Multimodal OCT for malignancy imaging | |
Yang et al. | Principles of doppler oct | |
Zhao et al. | Review of noninvasive methods for skin blood flow imaging in microcirculation | |
An et al. | Label-free optical micro-angiography for functional imaging of microcirculations within tissue beds in vivo |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 06826303 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC, EPO FORM 1205A OF 24.07.09 |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 06826303 Country of ref document: EP Kind code of ref document: A1 |