PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Opt Express. Author manuscript; available in PMC 2009 August 7.
Published in final edited form as:
PMCID: PMC2723161
NIHMSID: NIHMS129051

Velocity-resolved 3D retinal microvessel imaging using single-pass flow imaging spectral domain optical coherence tomography[diamond with plus]

Abstract

We demonstrate in vivo velocity-resolved, volumetric bidirectional blood flow imaging in human retina using single-pass flow imaging spectral domain optical coherence tomography (SPFI-SDOCT). This technique uses previously described methods for separating moving and non-moving scatterers within a depth by using a modified Hilbert transform. Additionally, a moving spatial frequency window is applied, creating a stack of depth-resolved images of moving scatterers, each representing a finite velocity range. The resulting velocity reconstruction is validated with and strongly correlated to velocities measured with conventional Doppler OCT in flow phantoms. In vivo velocity-resolved flow mapping is acquired in healthy human retina and demonstrate the measurement of vessel size, peak velocity, and total foveal blood flow with OCT.

1. Introduction

Spectral domain optical coherence tomography (SDOCT) has demonstrated clinical potential for in vivo high-resolution and high-speed imaging of biological structures [1, 2]. Advances in Doppler SDOCT (DOCT) have demonstrated several image acquisition schemes that enable real-time, high-resolution, volumetric display of blood flow maps [38]. Recently, several methods have been introduced which focus on utilizing spatial frequency modulations across lateral scans to resolve moving scatterers within a B-scan [912]. Techniques have also been presented which use reference delay modulations to specify a velocity measurement range [13]; or uses temporal frequency shifts to improve velocity resolution in low signal conditions [14]. The theoretical bases behind these techniques are closely related to those used in complex conjugate resolved imaging [1523]. A subset of these techniques [9, 10, 12] have demonstrated improved flow detection over conventional DOCT, resolving in vivo microvascular networks. Compared with DOCT, however, this subset of techniques has sacrificed velocity resolution for improved flow sensitivity, providing only flow direction without speed.

Spatial frequency flow detection techniques [9, 10, 12] can be considered optical analogs to power Doppler (PD) ultrasonography [2426]. Developed as a method of improving the sensitivity of Doppler ultrasound, the analog to DOCT, PD reports the power of the Doppler signal within specified frequency windows instead of the mean frequency shift. The advantage arises from the representation of the power spectrum of random phase noise. Since the noise in the power spectrum is uniformly low, random phase variations can be filtered out by raising the sensitivity threshold above the noise floor. The Doppler signal in PD is represented as an integral of the power spectrum, which improves the sensitivity and detection range of moving scatterers at the expense of eliminating velocity information. PD is relatively insensitive to Doppler angle and phase wrapping since these factors only modify the distribution of the Doppler power spectrum, but the total integrated power remains constant. The Doppler signal in PD is separated from non-moving components by filtering out all power spectrum components without Doppler shifts, thus only imaging moving scatterers. The resulting PD signal is related to the number of moving scatterers producing the corresponding Doppler shifts. Algorithmically, PD is equivalent to the single-pass flow imaging SDOCT (SPFI-SDOCT) technique we recently demonstrated [12].

Spatial frequency filtering techniques [9, 10, 12], while able to provide 3D flow maps with improved acquisition speed and sensitivity compared with DOCT, lack velocity-resolved blood flow information provided by techniques such as DOCT [4] and laser Doppler velocimetry (LDV) [27]. In ophthalmology, noninvasive quantification of blood circulation in tissues can facilitate the description of retinal vascular changes prior to and during ocular and systemic disease. Here, we demonstrate an improvement to SPFI-SDOCT which allows for velocity-resolved flow mapping using a similar modified Hilbert transform algorithm. The spatial frequencies of moving scatterers detected by SPFI are windowed such that only a small spatial frequency range, corresponding to a finite velocity range, is reconstructed. The window is then scanned across all spatial frequencies, forming a stack of B-scans each corresponding to a different velocity. The stack can then be summed across the entire velocity range to reconstruct a single velocity-resolved B-scan. OCT volume datasets can be processed to yield three-dimensional, velocity-resolved flow maps featuring the flow detection sensitivity improvement of SPFI, and the velocity resolution of conventional DOCT. Flow velocities calculated using this technique is validated against measurements taken with DOCT in flow phantoms. Velocity-resolved SPFI is then used to reconstruct flow maps of in vivo human retina estimating, for the first time to our knowledge, total foveal blood flow using OCT.

2. Theory

The theoretical basis for SPFI-SDOCT [12], and similar techniques [9, 10, 13, 17], have been developed previously. Briefly, for a recorded interferometric B-scan, I [k,x], the spatial Fourier transform along the lateral scan dimension (Fig. 1(a)) yields the spatial frequency content

Fig. 1
Flow chart of velocity-resolved SPFI-SDOCT signal processing. (a) Lateral Fourier transform of raw SDOCT spectral interferogram B-scan yields (b) spatial frequency of stationary scatterers centered around DC, and spatial frequency of moving scatterers ...
FTxu[I[k,x]]=R[k,u]+R¯[k,u]+V±[k,u+fD,±]+V±¯[k,ufD,±].
(1)

Here, both moving, V±[k,u + fD], and non-moving, R′[u,k], scatterers in the sample are imaged (Fig. 1(b)), where fD,± represents the Doppler frequency shift associated with the axial components of scatterer motion. All scatterer motion and associated flow refer to the axial components of motion.

In SPFI-SDOCT, we recognize that the spatial frequencies of moving and non-moving scatterers do not overlap at spatial frequencies above the non-moving scatterer bandwidth (Fig. 1(b)). An analytic signal for the spectral interferogram can be obtained by applying a Heaviside function (H[ufT]), frequency-shifted outside of the structural bandwidth (Fig. 1(c)), and then inverse Fourier transforming the result (Fig. 1(d)). Application of this modified Hilbert transform (HT*) enables bidirectional flow imaging by windowing Eq. (1) to yield

FTxu[I[k,x]]H[ufT]=α(V+[k,u+fD,+]+V¯[k,ufD,])
(2)

where α represents the fractional portion of bidirectional flow with Doppler frequencies outside of the spatial frequency bandwidth of non-moving scatterers (Fig. 1(c)). This threshold frequency (fT) defines the minimum detectable velocity in SPFI-SDOCT and is related to the spatial correlation of sequential A-scans [12]. The spatial oversampling and velocity resolution can be related by combining the velocity-related Doppler frequency shift with the spatial frequency resolution, thus yielding a velocity resolution

dv=Lλ02nDτwocosθD
(3)

where λ0 is the center wavelength, τ is the integration time, and θD is the Doppler angle between the scanning beam and the direction of scatterer motion. n is the index of refraction, L is the lateral scan length, D is the number of A-scans acquired across the lateral scan, and wo is the scanning beam spot size. Eq. (3) shows that velocity resolution in SPFI increases and the maximum detectable velocity decreases linearly with increased spatial oversampling.

Since velocity depends linearly on frequency, application of a spatial frequency window will necessarily filter out all velocities not traveling at velocities described by

Δv=Lλ0W[u]nDτcosθD.
(4)

Here, W [u] represents the spatial frequency FWHM bandwidth of the Gaussian window. There is, however, a tradeoff between the spatial frequency window width and the resulting spatial resolution of the velocity-resolved intensity image. For a Gaussian window in frequency-space, the resulting spatially resolved vessel map will be convolved with a Gaussian in the B-scan dimension with FWHM bandwidth W [x] = 4ln 2/(πW [u]). Thus, to avoid loss of lateral resolution in the velocity-resolved vessel maps, the spatial frequency window bandwidth, W [u], needs to be constrained such that the associated lateral blurring function, W [x], does not exceed the scanning beam spot size, wo. The velocity window limit can be represented as the bandwidth W [u] ≥ 4ln 2/(πwo). Given the discrete sampling parameters of SPFI-SDOCT, the minimum W [u] before loss of resolution is a factor of 8ln 2/π greater than the spatial frequency sampling rate. Combining this result with Eq. (3) gives the minimum resolvable velocity resolution as

dv=4ln2Lλ0πnDτwocosθD.
(5)

Finally, a datacube, I [x, z, v], for each B-scan can be formed by shifting the spatial frequency window across the bandwidth of moving scatterers (Fig. 1(c)), inverse transforming the result back to functions of I[k, x] (Fig. 1(d)), and then performing a spectral inverse Fourier transform (Fig. 1(e)). Each v-slice represents a velocity range given by Eq. (4), and the datacube can be summed across the v-dimension to create a single velocity-resolved B-scan. Similar to conventional DOCT, velocity wrapping occurs as the spatial frequency content of moving scatterers wraps across the Nyquist sampling upper limit. In this case, the scatterers are mapped to the opposite image half-plane and are therefore represented as moving in the opposite direction. Phase unwrapping techniques, similar to those used for DOCT [28], can be applied for a singly wrapped velocity profile. Higher spatial frequencies that wrap across Nyquist and have components greater than the threshold frequency, fT, will necessarily lose their velocity components due to SPFI windowing. This fundamentally limits the resolvable spatial frequency range to twice the Nyquist frequency.

3. Methods

Velocity-resolved SPFI-SDOCT was implemented on a high-speed SDOCT retinal imaging system (Fig. 2) employing a light source with a central wavelength of 859nm and a FWHM bandwidth of 99nm (Superlum, Ltd.). The sample arm was a modified slit lamp equipped with scanning galvanometers and relay optics for retinal imaging of subjects. The retinal scanner optics were designed for a 15–20μm transverse resolution, as limited by the optics of the eye, across a 12×12mm field. The reference arm was dispersion compensated using a water cell and matched optics, and interferometric signals were captured using a 2048 pixel line-scan camera (e2v, Ltd.). Custom software (Bioptigen, Inc.) performed data acquisition, archiving, and real-time processing and display of image magnitude. Using a 700μW sample beam, the SNR measured near DC was 110dB with an axial resolution of 4.72μm in tissue and a 6dB falloff at 0.8mm. DC removal, k-space resampling, and flow imaging using the modified Hilbert transform algorithm [12] were computed during post-processing using Matlab (MathWorks, Inc.). Vessels and structure were visualized using Amira (Visage Imaging, Inc.) and OSA ISP (Kitware, Inc.).

Fig. 2
Velocity-resolved SPFI-SDOCT retinal system. The sample arm is a modified slit lamp equipped with scanning galvanometers and relay optics for convenient patient retinal imaging. SPFI windowing and velocity-map reconstruction were calculated in post-processing ...

Velocity-resolved bidirectional flow imaging was validated and compared with conventional DOCT on a flow phantom. Two glass micro-capillary tubes (1.5mm outer diameter, 0.6mm inner diameter) were connected using silastic tubing to a syringe pump (Harvard Apparatus) and pumped with 1% liposyn at 10μL/min, 20μL/min, 30μL/min, and 40μL/min. The micro-capillaries were then positioned adjacent to each other on an angled stage such that fluid in the tubes flowed in opposite directions in a B-scan cross-section (Fig. 3(a)), simulating bidirectional flow. B-scans of the phantom were acquired across a 2mm scan range with 2500 A-scans/frame for SPFI and 1000 A-scans/frame with 4 sequential A-scans at each lateral position for DOCT. The SPFI dataset was laterally oversampled compared to DOCT because velocity resolution increases as a function of spot-size overlap on the sample (Eq. (3)). Both datasets were acquired with an integration time of 50μs. At these sampling parameters, the total imaging time for SPFI is a factor of 1.6 times faster than that of DOCT. SPFI parameters were chosen to demonstrate velocity resolution at the lower limit of the detection range for a given integration time as a comparison with DOCT. Since the velocity resolution of DOCT is limited by the phase noise of the system instead of spatial sampling parameters, as is the case with SPFI, an appropriate lateral spacing used to minimize scanner jitter. The number of sequential A-scans used is indicative of common DOCT sampling parameters.

Fig. 3
Validation of SPFI measured flow velocities. (a) Two micro-capillaries were connected and oriented such that 1% liposyn flowed in opposite directions in a B-scan cross-section. (b),(c) Velocity profiles measured using DOCT and fit to laminar flow curves. ...

A threshold frequency was determined which filtered out all spatial frequencies of non-moving scatterers. A Gaussian window was then moved across the remaining spatial frequencies to velocity-resolve the B-scan for flow rate measurements. The shifted window was set such that the velocity-range resolved had a FWHM and shifted at increments of 24.3μm/s. Given the oversampling parameters and threshold frequency window used, the magnitude of the total detectable velocity range for the axial components of positive and negative velocities was 0.61–11.53mm/s. Velocity-resolved scatterer information was then overlaid onto structural B-scans for visualization (Fig. 3(a)). DOCT volumes were processed using standard phase-difference methods [6]. Velocity profiles for both SPFI and DOCT were fit to laminar flow curves, and measured capillary cross-section and flow rate were calculated for both imaging methods and compared.

In vivo microvessel imaging was demonstrated in normal human retina. First, a 10×10mm OCT volume of the subject was acquired, allowing for the reconstruction of a standard OCT summed voxel projection (SVP) to use as an atlas to locate smaller volumes imaged using SPFI (Fig. (4)). Several 2×2mm volumes were then densely-sampled using SPFI parameters (2500 A-scans/frame, 100 frames/volume, 100μs integration time − 25s total imaging time) at several locations across the macula (Figs. 4(a) (View 1), 4(b) (View 2), and 4(d) (View 4)), including landmarks such as the fovea (Fig. 4(c)(View 3)) and optic nerve (Fig. 4(e)(View 5)). Given the sampling parameters, threshold frequency, and assuming a 20μm scanning beam spot-size on the retina, the detectable velocity range for axially moving scatterers was 0.45–8.64mm/s. Parameters were set for a lower velocity range, compared to flow phantoms, to adequately detect slow flow in small foveal vessels (<30mm/s) [27]. The frequency window was set such that the velocity-range had a FWHM of 18.4μm/s and was shifted at increments of 82.8μm/s. Bulk motion correction using the algorithm described in [10] was implemented prior to SPFI windowing. Velocity-resolved B-scans were used to determine vascular size by fitting velocity profiles to parabolic flow curves and then calculating the zero-velocity crossing positions. Finally, blood flow rates and total retinal flow were calculated for a 2×2mm velocity-resolved volume of human fovea.

Fig. 4
10×10mm volume of in vivo human retina. (a)–(e) 2×2mm volumes sampled with 1024×2500×100pix at 100μs integration time (25s total imaging time) were acquired at several locations across the macula, (a) ( ...

SPFI acquisition of the entire macula in a single volume dataset would be impractical. Since velocity resolution and lateral sampling density are coupled in SPFI by Eq. (3), an order of magnitude increase in the lateral sampling density, A-scans/frame, would be required. This requirement limits the SPFI sampling volume size by both the available memory in the acquisition software and total imaging time.

4. Results and Discussion

Velocity-resolved bidirectional flow imaging was demonstrated in a liposyn-pumped flow phantom and compared with measurements taken with conventional DOCT (Fig. 3). Velocity profiles for both positive and negative flow were plotted for both SPFI and DOCT datasets and fit to second-order polynomials, characteristic of laminar flow, to determine the tube diameter and maximum flow velocities (Fig. 3(b)–(d)). The DOCT measured velocity profiles correlated stronger with the parabolic flow fit (mean R2=0.96) than those calculated with SPFI (mean R2=0.81). Since laminar flow velocities across the tube decrease continuously towards zero, all components below the SPFI velocity threshold contribute to velocity noise, which becomes more significant for low flow velocities. The velocities measured using SPFI shows a poor parabolic fit and a correlation of R2<0.7 for the lowest pump speed (Fig. 3(d–e) – dark blue), demonstrating the lower detection limit of the modality. This is a result of the spatial frequency overlap of the slow moving scatterers near the edge of the capillary tubes, which is filtered out in SPFI processing. As scatterer velocity increases, the fit improves and is comparable with that of DOCT (Fig. 3(f–g)). Fig. 3(d) also shows a velocity dip across all pump velocities calculated using SPFI. This is a result of a shadowing artifact as seen in Fig. 3(a), which is more pronounced in SPFI and in DOCT (Fig. 3(b)) since SPFI is an intensity sensitive modality, as opposed to phase sensitive DOCT. For each flow rate, velocity profiles were measured three times for both SPFI and DOCT, and the averaged velocity and standard deviation were calculated (Fig. 3(b)–(g)). Flow rates measured using SPFI and DOCT were plotted as functions of pump flow rates (Fig. 3(f)–(g)), showing strong correlations of both SPFI (mean R2=0.99) and DOCT (mean R2=0.99) with theoretical values.

Depth-resolved vessel maps for each volume were first reconstructed and overlaid onto structural OCT data to distinguish macular vasculature (Fig. 4(a)–(e)). Total blood flow measurements were then calculated for a single 2×2mm volume of the fovea (Fig. 4(c)(View 3)). The structural OCT data shows the foveal pit and the associated SPFI vascular map confirmed its location by resolving a circular avascular zone surrounded by a set of terminal capillaries [29]. First, the resolved vessels in the volume were identified as arteries and veins using the flow directionality information calculated by SPFI (Fig. 5(a)). Vessel orientation and Doppler angle were then measured for all 17 resolved vessels, and velocity profiles were measured at a single point for each vessel (Fig. 5(a) – dots). Doppler angles were measured by calculating the vessel cross-section displacement across sequential B-scans in the 3D datasets [4]. The velocity profiles were then fit to laminar flow curves to determine both the peak velocity and diameter of each vessel. Finally, the vessel size, orientation, and velocity information was used to calculated the total foveal blood flow [4]. A representative B-scan (Fig. 5(b)) taken across the foveal volume (dotted line) is shown with velocity-resolved flow content overlaid on top of the structural image. The velocity profiles and parabolic fits for each of the vessels are included (Fig. 5(c)–(i)) to show the strong correlation between velocities measured using SPFI and their respective laminar flow velocity profiles (mean R2=0.95). Closer inspection of vessel 17 (Fig. 5(i)) shows a blunted parabolic velocity profile, characteristic of red blood cell aggregation in microvasculature, consistent with rheological observations [3032].

Fig. 5
Velocity-resolved 2×2mm volumetric flow map of in vivo human fovea. (a) Fovea vessel map with vessels and flow direction resolved using SPFI. Arrows denote arteries and veins, and dots represent locations where velocity profiles were measured. ...

A summary of the size, peak velocity, and flow measurements for all vessels identified in the 2×2mm foveal volume are shown in Table 1. The smallest resolvable vessel was 13.64μm, which is at the resolution limit of the retinal SDOCT system. The average measured foveal vessel diameter was ~22μm, and average arterial flow velocity was greater than average venous flow velocity (Table 2), which is supported by similar measurements made using LDV [27]. The detected velocities ranged from 5.97–30.22mm/s, concurrent with human retinal vessels in the corresponding size range, also measured using LDV [27]. Finally, the total arterial and venous flow showed a net inflow of blood into the fovea. This is a result of the presence of unresolved veins in the volume and expected errors from measurements of vessel orientation angle and calculation of diameter, which would significantly impact blood flow calculations.

Table 1
Size, velocity, and flow measurements across a 2×2mm volumetric flow map of in vivo human fovea calculated using SPFI.
Table 2
Retinal flow statistics in a 2×2mm volumetric flow map of in vivo human fovea calculated using SPFI-SDOCT.

While SPFI allows for velocity-resolved volumetric blood flow imaging with velocity resolution comparable to that of DOCT, there are several advantages and drawbacks to the modality that must be noted. Since velocity resolution and spatial oversampling are coupled, this allows for the flexibility to set sampling parameters to a desired velocity range. This would allow for improved data acquisition speeds over DOCT for imaging of small spatial volumes of moderately high flow velocities, but could potentially require longer scanning times for large scan areas with low flow. As discussed previously, the velocity resolving power of SPFI falls off at the lower limits of velocity detection, as compared to DOCT, but is comparable for high flow velocities. Since velocity-resolved SPFI requires a sliding spatial frequency window, this increases the number of Fourier transforms required by the number of velocity increments desired. However, this increase in computational complexity is essentially trivial since FFT algorithms can be optimized and the transforms can be parallelized. Finally, it has been shown that spatial frequency filtering results in a sensitivity improvement [16], which gives SPFI a detection advantage over DOCT for small vessels at the resolution limit of the imaging system. The overall advantage of SPFI over DOCT is in its improved sensitivity and customizability over velocity resolution, velocity range, and acquisition time in exchange for small increases in computational complexity.

5. Conclusions

We have demonstrated velocity-resolved in vivo volumetric bidirectional flow imaging for the measurements of human retinal blood flow. Velocity-resolved SPFI-SDOCT builds on previously described SPFI imaging techniques, which provide three-dimensional flow imaging without acquiring multiple A-scans at a single lateral position. SPFI measured velocities were validated with velocities measured using conventional DOCT on flow phantoms. Velocity-resolved SPFI was then used for in vivo imaging of several volumes across the macula of healthy human retina. SPFI processing allowed for volumetric visualization of flow without the need for manual segmentation. SPFI also provided flow directionality and velocity, demonstrating the measurement of vessel size, peak velocity, and total foveal blood flow with OCT, consistent with characteristics observed in LDV and rheology.

Acknowledgments

This research was supported by NIH grant R21 EY017393.

Footnotes

[diamond with plus]Datasets for this article are available at http://hdl.handle.net/10376/1400.

References and links

1. Nassif NA, Cense B, Park BH, Pierce MC, Yun SH, Bouma BE, Tearney GJ, Chen TC, de Boer JF. In vivo high-resolution video-rate spectral-domain optical coherence tomography of the human retina and optic nerve. Opt Express. 2004. pp. 367–76. http://www.opticsinfobase.org/abstract.cfm?URI=oe-12-3-367. [PubMed]
2. Wojtkowski M, Srinivasan VJ, Ko TH, Fujimoto JG, Kowalczyk A, Duker JS. Ultrahigh-Resolution, high-speed, Fourier domain optical coherence tomography and methods for dispersion compensation. Opt Express. 2004. pp. 2404–22. http://www.opticsinfobase.org/abstract.cfm?URI=oe-12-11-2404. [PubMed]
3. Mariampillai A, Standish BA, Munce NR, Randall C, Liu G, Jiang JY, Cable AE, Vitkin IA, Yang VXD. Doppler optical cardiogram gated 2D color flow imaging at 1000 fps and 4D in vivo visualization of embryonic heart at 45 fps on a swept source OCT system. Opt Express. 2007. pp. 1627–38. http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-4-1627. [PubMed]
4. Wang Y, Bower BA, Izatt JA, Tan O, Huang D. In vivo total retinal blood flow measurement by Fourier domain Doppler optical coherence tomography. Biomed Opt. 2007;12:0412151–8. [PubMed]
5. Bower BA, Zhao M, Zawadzki RJ, Izatt JA. Real-time spectral domain Doppler optical coherence tomography and investigation of human retinal vessel autoregulation. J Biomed Opt. 2007;12:0412141–8. [PubMed]
6. Wang L, Wang Y, Guo S, Zhang J, Bachman M, Li GP, Chen Z. Frequency domain phase-resolved optical Doppler and Doppler variance tomography. Opt Commun. 2004;242:345–50.
7. Leitgeb RA, Schmetterer L, Drexler W, Fercher AF, Zawadzki RJ, Bajraszewski T. Real-time assessment of retinal blood flow with ultrafast acquisition by color Doppler Fourier domain optical coherence tomography. Opt Express. 2003. pp. 3116–21. http://www.opticsinfobase.org/abstract.cfm?URI=oe-11-23-3116. [PubMed]
8. White BR, Pierce MC, Nassif NA, Cense B, Park BH, Tearney GJ, Bouma BE, Chen TC, De Boer JF. In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical Doppler tomography. Opt Express. 2003. pp. 3490–7. http://www.opticsinfobase.org/abstract.cfm?URI=oe-11-25-3490. [PubMed]
9. Wang RK, Jacques SL, Ma Z, Hurst S, Hanson SR, Gruber A. Three dimensional optical angiography. Opt Express. 2007. pp. 4083–97. http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-7-4083. [PubMed]
10. An L, Wang RK. In vivo volumetric imaging of vascular perfusion within human retina and choroids with optical micro-angiography. Opt Express. 2008. pp. 11438–52. http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-15-11438. [PubMed]
11. Yang VXD, Gordon ML, Qi B, Pekar J, Lo S, Seng-Yue E, Mok A, Wilson BC, Vitkin IA. High speed, wide velocity dynamic range Doppler optical coherence tomography (Part I): System design, signal processing, and performance. Opt Express. 2003. pp. 794–809. http://www.opticsinfobase.org/abstract.cfm?URI=oe-11-7-794. [PubMed]
12. Tao YK, Davis AM, Izatt JA. Single-pass volumetric bidirectional blood flow imaging spectral domain optical coherence tomography using a modified Hilbert transform. Opt Express. 2008. pp. 12350–61. http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-16-12350. [PubMed]
13. Bachmann AH, Villiger ML, Blatter C, Lasser T, Leitgeb RA. Resonant Doppler flow imaging and optical vivisection of retinal blood vessels. Opt Express. 2007. pp. 408–22. http://www.opticsinfobase.org/abstract.cfm?URI=oe-15-2-408. [PubMed]
14. Szkulmowski M, Szkulmowska A, Bajraszewski T, Kowalczyk A, Wojtkowski M. Flow velocity estimation using joint Spectral and Time Domain Optical Coherence Tomography. Opt Express. 2008. pp. 6008–25. http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-9-6008. [PubMed]
15. Leitgeb RA, Michaely R, Lasser T, Sekhar SC. Complex ambiguity-free Fourier domain optical coherence tomography through transverse scanning. Opt Lett. 2007. pp. 3453–5. http://www.opticsinfobase.org/abstract.cfm?URI=ol-32-23-3453. [PubMed]
16. Yasuno Y, Makita S, Endo T, Aoki G, Itoh M, Yatagai T. Simultaneous B-M-mode scanning method for real-time full-range Fourier domain optical coherence tomography. Appl Opt. 2006. pp. 1861–5. http://www.opticsinfobase.org/abstract.cfm?URI=ao-45-8-1861. [PubMed]
17. Wang RK. In vivo full range complex Fourier domain optical coherence tomography. Appl Phys Lett. 2007;90:0541031–3.
18. Baumann B, Pircher M, Gotzinger E, Hitzenberger CK. Full range complex spectral domain optical coherence tomography without additional phase shifters. Opt Express. 2007. pp. 13375–87. http://www.opticsinfobase.org/abstract.cfm?URI=oe-15-20-13375. [PMC free article] [PubMed]
19. An L, Wang RK. Use of a scanner to modulate spatial interferograms for in vivo full-range Fourier-domain optical coherence tomography. Opt Letters. 2007. pp. 3423–5. http://www.opticsinfobase.org/abstract.cfm?URI=ol-32-23-3423. [PubMed]
20. Makita S, Fabritius T, Yasuno Y. Full-range, high-speed, high-resolution 1-mu m spectral-domain optical coherence tomography using BM-scan for volumetric imaging of the human posterior eye. Opt Express. 2008. pp. 8406–20. http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-12-8406. [PubMed]
21. Bachmann AH, Leitgeb RA, Lasser T. Heterodyne Fourier domain optical coherence tomography for full range probing with high axial resolution. Opt Express. 2006. pp. 1487–96. http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-4-1487. [PubMed]
22. Yun SH, Tearney GJ, de Boer JF, Bouma BE. Removing the depth-degeneracy in optical frequency domain imaging with frequency shifting. Opt Express. 2004. pp. 4822–8. http://www.opticsinfobase.org/abstract.cfm?URI=oe-12-20-4822. [PMC free article] [PubMed]
23. Zhang J, Nelson JS, Chen ZP. Removal of a mirror image and enhancement of the signal-to-noise ratio in Fourier-domain optical coherence tomography using an electro-optic phase modulator. Opt Lett. 2005. pp. 147–9. http://www.opticsinfobase.org/abstract.cfm?URI=ol-30-2-147. [PubMed]
24. Rubin JM, Bude RO, Carson PL, Bree RL, Adler RS. Power Doppler US: a potentially useful alternative to mean frequency-based color Doppler US. Radiology. 1994;190:853–6. [PubMed]
25. Rubin JM. Power Doppler. European Radiology. 1999;9:S318–22. [PubMed]
26. Rubin JM, Adler RS. Power Doppler expands standard color capability. Diagnostic Imaging. 1993;15:66–9. [PubMed]
27. Riva CE, Grunwald JE, Sinclair SH, Petrig BL. Blood velocity and volumetric flow rate in human retinal vessels. IOVS. 1985;26:1124–32. [PubMed]
28. Ghiglia DC, Mastin GA, Romero LA. Cellular-Automata Method for Phase Unwrapping. J Opt Soc Am. 1987. pp. 267–80. http://www.opticsinfobase.org/abstract.cfm?URI=josaa-4-1-267.
29. Engerman RL. Development of the macular circulation. Investigative Oophthalmology. 1976;15:835–40. [PubMed]
30. Bagchi P. Mesoscale simulation of blood flow in small vessels. Biophysical J. 2007;92:1858–77. [PubMed]
31. Bishop JJ, Nance PR, Popel AS, Intaglietta M, Johnson PC. Effect of erythrocyte aggregation on velocity profiles in venules. Am J Physiol. 2001;280:H222–36. [PubMed]
32. Sharan M, Popel AS. A two-phase model for flow of blood in narrow tubes with increased effective viscosity near the wall. Biorheology. 2001;38:415–28. [PubMed]