17 October 2017 Spectroscopic Doppler analysis for visible-light optical coherence tomography
Author Affiliations +
J. of Biomedical Optics, 22(12), 121702 (2017). doi:10.1117/1.JBO.22.12.121702
Abstract
Retinal oxygen metabolic rate can be effectively measured by visible-light optical coherence tomography (vis-OCT), which simultaneously quantifies oxygen saturation and blood flow rate in retinal vessels through spectroscopic analysis and Doppler measurement, respectively. Doppler OCT relates phase variation between sequential A-lines to the axial flow velocity of the scattering medium. The detectable phase shift is between π and π due to its periodicity, which limits the maximum measurable unambiguous velocity without phase unwrapping. Using shorter wavelengths, vis-OCT is more vulnerable to phase ambiguity since flow induced phase variation is linearly related to the center wavenumber of the probing light. We eliminated the need for phase unwrapping using spectroscopic Doppler analysis. We split the whole vis-OCT spectrum into a series of narrow subbands and reconstructed vis-OCT images to extract corresponding Doppler phase shifts in all the subbands. Then, we quantified flow velocity by analyzing subband-dependent phase shift using linear regression. In the phantom experiment, we showed that spectroscopic Doppler analysis extended the measurable absolute phase shift range without conducting phase unwrapping. We also tested this method to quantify retinal blood flow in rodents in vivo.
Shu, Liu, Duan, and Zhang: Spectroscopic Doppler analysis for visible-light optical coherence tomography

1.

Introduction

Retinal metabolic rate of oxygen (rMRO2) is a parameter critical for the fundamental investigation and clinical diagnosis of several blinding diseases such as age-related macular degeneration, glaucoma and, most importantly, diabetic retinopathy (DR).12.3 Previous studies suggested that retinal metabolic abnormalities play an important role in the development of DR2,4,5 and that rMRO2 is a potential biomarker for early diagnosis of DR before irreversible damage occurs.67.8 In order to calculate rMRO2, we need to measure oxygen saturation (sO2) and the flow rate in major retinal blood vessels, both of which can be effectively quantified simultaneously by visible-light optical coherence tomography (vis-OCT).910.11 As a single imaging modality for rMRO2 measurement, vis-OCT detects sO2 by retrieving the absorption spectrum of blood from spectroscopic analysis and measures the flow rate using the established Doppler OCT method.1213.14 Compared with conventional OCT using near-infrared (NIR) light, vis-OCT provides reliable sO2 measurement and higher resolution, but has reduced measurement range of flow velocity.13

This drawback of vis-OCT in flow measurement originates from the underlying mechanism of Doppler OCT, which relates phase shift between sequentially acquired OCT A-lines to the axial displacement of scatterers. The axial displacement is converted to flow velocity based on the A-line rate of system and Doppler angle, which is the angle between the flow direction and light probing axis.1516.17.18 Since phase variation is periodic, the actual detectable Doppler phase shift always lies between π and π. Phase wrapping occurs when flow induced axial displacement exceeds the center wavelength of the probing beam.15,17 Therefore, using shorter wavelengths, vis-OCT is more sensitive to low flow velocity, but more susceptible to phase wrapping in measuring high flow velocity.

There are a few potential ways to extend the measurable velocity range of vis-OCT. First, phase unwrapping can be applied to increase the range beyond π and π. However, although various unwrapping techniques have been investigated, it remains challenging to remove artifacts induced by phase wrapping in clinical settings.19,20 Second, the phase shift induced by a certain flow velocity can be reduced by decreasing the time interval between two sequential observations, which is the inverse of A-line acquisition rate. This can be done by decreasing the spectrometer integration time of spectral-domain OCT (SD-OCT) or by switching to a fast swept-source system. However, a reduced spectrometer integration time leads to a lower signal-to-noise ratio (SNR) and a swept-source for visible light is not yet available. Third, the imaging operator can intentionally position the OCT probing beam and retinal blood vessel for a larger Doppler angle to decrease the projected axial velocity. However, this procedure can significantly increase the operation complexity and may be impractical for the region around the optic nerve head.

We seek to extend the flow velocity measurement range for vis-OCT without using sophisticated phase unwrapping. In most existing Doppler OCT algorithms, the center wavelength or wavenumber of the broadband light source is used to translate the detected phase shift to the axial displacement of scatters.15,21 However, if we split the broad spectrum into a number of narrower subbands for OCT reconstruction, a series of Doppler phase shifts can be obtained and these phase shifts follow a linear relationship with the center wavenumbers of the subbands.12,22,23 The slope of this linear function, instead of the actual phase shift values, can be used to retrieve the axial flow velocity of the scatterers, which is independent of phase wrapping.15,21 Hence, spectroscopic Doppler analysis may extend the maximum measurement limit beyond π and π, which is especially beneficial in vis-OCT.

In this paper, we report a spectroscopic Doppler analysis method to quantify flow velocity in spectral-domain vis-OCT. We performed a series of short-time Fourier transforms (STFT) of the original broadband OCT interferogram and obtained Doppler phase shifts from different subbands. We applied linear regression on the wavenumber-dependent phase shift values and retrieved the axial flow velocity from the fitted slope. We validated our algorithm using a vis-OCT system with a 100-nm bandwidth at a center wavelength of 560 nm. Phantom experiments showed that both conventional Doppler and the spectroscopic Doppler analysis methods measured flow velocities accurately when the velocity was low; but spectroscopic Doppler analysis increased the measurable range of the flow velocity. We also demonstrated the spectroscopic Doppler analysis method in measuring in vivo retinal blood flow in rodents.

2.

Method of Spectroscopic Doppler Analysis

2.1.

Principle of Doppler Optical Coherence Tomography

SD-OCT measures flow velocity by detecting the phase difference between two sequential A-lines.21 If a single moving scatterer appears at depths of z0 and z0+δz in the first and second A-lines, we can describe the corresponding interferograms i(k) of the two measurements as

(1)

i0(k)=ξs(k)RRRScos(2nkz0),
and

(2)

i1(k)=ξs(k)RRRScos[2nk(z0+δz)],
where ξ is a constant determined by the detector quantum efficiency; s(k) is the spectral power density function of the source; RR and RS are the reflectivities of the reference mirror and the sample scatterer, respectively; n is the refractive index; and k is the angular wavenumber, defined as 2π/λ (λ is the wavelength). Applying Fourier transform to i0(k) and i1(k), we can obtain the complex OCT A-lines I0(z) and I1(z). If δz is much smaller than OCT axial resolution, the modulus of complex A-line, corresponding to OCT structure image, is insensitive to the scatter displacement. However, a phase shift term is added to I1(z), which is proportional to δz

(3)

I0(z)=ξRRRS[S(z2nz0)+S(z+2nz0)],

(4)

I1(z)=ξRRRS[S(z2nz0)ej2nkequδz+S(z+2nz0)ej2nkequδz],

(5)

ΔφD=2nkequδz,
where S(z) is the autocorrelation of s(k); and kequ is the equivalent wavenumber of the light source, which can be calculated based on Eq. (20) in Ref. 24 as

(6)

kequ=k·s(k)s(k).

For a light source with a symmetric spectrum, center wavenumbers kc and kequ are interchangeable in Eq. (6). However, in real cases when the source spectrum is asymmetric, using kc directly leads to bias in flow velocity measurement.

We use T to represent the acquisition time interval between two sequential A-lines, which is usually determined by the spectrometer integration time in SD-OCT or the wavelength sweeping cycle in a swept-source OCT (SS-OCT). We can rearrange Eq. (5) and obtain the axially projected flow velocity v by taking the ratio of the axial displacement δz and T

(7)

v=ΔφD2nTkequ.

Phase wrapping happens when the axial displacement induced absolute phase shift exceeds π, which results in ambiguous axial flow velocity measurement. The phase wrapping limited maximum absolute axial velocity (PLV) is given by15,18,21

(8)

v,wrap=π2nTkequ.

Another factor that limits the flow detection is the washout effect, which describes the signal loss due to averaging of motion-varied interferometric fringes over the integration time of the optical detector. For SD-OCT, the washout effect starts at around PLV and will cause a Doppler phase signal unrecoverable beyond twice PLV.19

2.2.

Principle of Spectroscopic Doppler Optical Coherence Tomography

STFT has been used in spectroscopic analysis of OCT signals to extract wavenumber-dependent absorption and scattering properties in tissue,2425.26 which can be described as

(9)

STFT{i(k)}(τ,ω)=DFT[i(k)w(k,τ,ω)],
where τ and ω are the center wavenumber and bandwidth of each OCT subband, respectively; DFT is discrete Fourier transform; and w(·) represents the Gaussian window function. Since the Gaussian window gives more weight to the central wavelength in the subband, the equivalent wavenumber defined in Eq. (6) can be estimated by the center wavenumber τ. If the axial flow velocity is v, the detected Doppler phase shift in each subband can be calculated from Eq. (7) as

(10)

ΔφD(τ)=2nTvτ,
which shows that a series of ΔφD(τ) can be obtained by sliding the window function across the broadband OCT spectrum, and that the Doppler phase shift is a linear function of the center wavenumber with its slope proportional to the axial flow velocity. If we use a to represent the slope of this linear function, the axial flow velocity can be calculated as

(11)

v=a2nT.

Therefore, we can use spectroscopic analysis and linear fitting to obtain a and, therefore, v.

We performed numerical simulation to better demonstrate how Doppler phase shift varies with respect to center wavenumber of the OCT subband. Though OCT reconstruction and Doppler analysis are achieved in wavenumber domain, researchers often refer to OCT bandwidth in wavelength. Therefore, we provided both wavenumber and wavelength values in the following description. The simulated OCT source had a square-shaped spectrum with (0.19×105)rad/cm (100 nm) bandwidth and (1.10×105)rad/cm (570 nm) center wavelength. The Gaussian windows used in STFT had (0.02×105)rad/cm (9 to 12 nm, depending on the wavelength) standard deviation, and the center wavenumbers ranged from (1.04×105)rad/cm (604 nm) to (1.16×105)rad/cm (542 nm) at 0.03  rad/cm interval, as shown in Fig. 1(a). In Sec. 2.1, we use z0 and z0+δz to represent the depth locations of a moving scatter in the two sequential OCT A-lines and δz is the axial displacement. In the simulation, we simulated seven flow velocities whose axial displacements between two sequential A-lines were 0,δz0,2δz0,,6δz0, where δz0 was 35.7 nm. According to Eq. (5), 6δz0 would lead to a 2π Doppler phase shift assuming that the center wavelength was 570 nm and that the refractive index was 1.33. Figure 1(b) shows the subband-dependent Doppler phase shifts under different axial displacements, where the colored dots indicate the measured phase shifts and the dashed black lines show the fitted linear relationship. Each color represents one of the OCT subband in Fig. 1(a). One would note that 0 and 6δz0 displacements would both yield zero flow velocity if measuring the Doppler phase shift with the whole spectrum due to phase wrapping. However, the slope of subband-dependent phase shift array is uniquely related to displacement values, unaffected by phase wrapping. Another interesting phenomenon occurs in the case of 3δz0, where a 2π shift exists in the B3 subband. This is because the displacement induced phase shifts are so close to π that phase wrapping occurs among different subbands. We call this phenomenon spectroscopic phase shift discontinuity, which needs to be corrected before linearly fitting different phase shift values with subband center wavenumbers. Although this phase shift discontinuity can also happen when flow induced phase shifts are around 3π, 5π, and etc., washout effect limits the measurable axial flow velocity within twice PLV, corresponding to 2π Doppler phase shift. Therefore, we only need to correct the spectroscopic phase shift discontinuity that happens around π.

Fig. 1

Simulation of spectroscopic Doppler OCT analysis. (a) Selection of Gaussian windows for subbands. (b) Subband-dependent Doppler phase shifts under various axial displacements, from 0 to 6δz0. The dotted and the dashed lines are the phase shifts and the fitted lines for each single displacement. The colors represent the different Gaussian windows corresponding to (a).

JBO_22_12_121702_f001.png

2.3.

Flow Velocity Calculation in Spectroscopic Doppler Optical Coherence Tomography

The major application of Doppler OCT is to measure blood flow velocity by averaging Doppler phase shift values of all pixels within a blood vessel cross section. In cases where the flow velocity exceeds PLV, phase unwrapping needs to be implemented in conventional Doppler OCT to obtain correct the measurement result. If the flow velocity is constant throughout the blood vessel, phase unwrapping can be done simply by adding or subtracting 2π to all the pixels. However, the actual flow profile within the vessel follows an approximate three-dimensional paraboloid. When only the flow velocity in the center of a blood vessel exceeds PLV, it is challenging to find the boundary between the wrapped and nonwrapped areas. The state-of-the-art two-dimensional algorithms to correct this type of phase wrapping are usually iterative processes and often generate biased offset error in noisy datasets.2728.29

In spectroscopic Doppler OCT, we can retrieve flow velocity based on the linear relationship between phase shifts and center wavenumbers of OCT subbands. Phase wrapping only affects the intercept of the fitted line without changing its slope. Since we extract the flow velocity from the slope of the fitted line, phase unwrapping is not required. But we do have to correct the spectroscopic phase shift discontinuity, which is demonstrated in Fig. 1(b). We summarize the major steps of spectroscopic Doppler analysis in Fig. 2.

  • Step 1: STFT

    We perform STFT on the broadband OCT B-scan and obtain a set of subband B-scans. We calculate phase shift values in each subband image in the same way as conventional Doppler OCT processing. We segment the region within blood vessel from one subband phase shift image and apply the obtained boundary to all subbands.

  • Step 2: Correction for phase shift discontinuity

    In this step, we aim to calculate the average phase shifts within the blood vessel for each subband. If there is no spectroscopic phase shift discontinuity, the average phase shift Δφ¯D(τ) can be calculated as follows:

    (12)

    Δφ¯D(τ)=ΔφD[τ,(x,z)]N,
    where τ is the center wavenumber of subband used in STFT, N is the number of pixels inside a blood vessel used for spatial averaging, and (x,z) is the coordinate of the pixel in the B-scan image.

    However, if the phase shift values of some pixels are close to π, spectroscopic phase shift discontinuity may occur. To preserve the linear relationship, the average phase shift should be calculated as follows instead

    (13)

    Δφ¯D(τ)=ΔφD[τ,(x,z)]N+Δφcorr[τ,(x,z)]N,
    where a correction term is added to compensate for the discontinuity. To calculate the correction term, we select a phase shift image reconstructed on an arbitrary subband τ0 as the reference and subtract the reference image from phase shift images of all the subbands as

    (14)

    ψ[τ,(x,z)]=|ΔφD[τ,(x,z)]ΔφD[τ0,(x,z)]|.

    ψ[τ,(x,z)] is expected to be small if no discontinuity occurs. Therefore, we compensate the discontinuity by adding or subtracting 2π when ψ[τ,(x,z)] exceeds a threshold value

    (15)

    Δφcorr[τ,(x,z)]={2π·sgn{ΔφD[τ0,(x,z)]},ψ[τ,(x,z)]>th,0,ψ[τ,(x,z)]th,
    where sgn(·) is the sign function. We mentioned in Sec. 2.2 that the spectroscopic phase shift discontinuity only happens when the flow induced phase shift is around π and that the maximum measurable phase shift should be smaller than 2π due to the washout effect. Therefore, we define the threshold value th as

    (16)

    th=2π(ττ0)τ0.

  • Step 3: Linear fitting

    We do linear regression to the subband-dependent average phase shift values using least-squares fit and extract the slope of the fitted linear relationship.

    (17)

    a=dΔφ¯D(τ)dτ.

Fig. 2

Flowchart of spectroscopic Doppler analysis.

JBO_22_12_121702_f002.png

Then, the axial flow velocity is calculated by Eq. (11).

3.

Experimental System

We demonstrated flow velocity measurement in vis-OCT using spectroscopic Doppler analysis. The system schematic is shown in Fig. 3(a). We used a broadband supercontinuum laser (Superk EXTREME, NKT Photonics) to generate a 90-nm bandwidth visible light source centered at 565 nm, which is also used in clinical settings.11,30 The source beam is coupled into a 2×2 fiber coupler (Nufern 460-HP, 50:50 splitting ratio, GouldFiber Optics) which delivers the light to a sample arm and a reference arm. In the sample arm, the illumination light is collimated and scanned by a set of galvanometer mirrors (6210H, Cambridge Technology). In phantom experiments, we focused the illumination beam using an achromatic doublet (Edmund Optics) with a focal length of 35 mm [Fig. 3(b)]. In in vivo experiments, we used a pair of achromatic doublets with focal lengths of 75 and 15 mm, respectively, to relay the illumination beam onto the cornea [Fig. 3(c)]. In the reference arm, after collimation, the beam propagates through a BK7 glass block for dispersion compensation and is reflected by a mirror. The backscattered light from the sample arm interferes with the back-reflected light from the reference arm in the fiber coupler and the interferogram is collected by a home-built spectrometer, consisting of a transmission grating (1800  line/mm, Wasatch Photonics), an SLR lens (85 mm, f/1.4, Samyang), and a line-scan camera (spL4096-70 km, Basler). Though there are 4096 linear pixels in the camera, we only use the center 2048 for spectrum acquisition.

Fig. 3

(a) Schematic of vis-OCT system. SC, supercontinuum source; SM, spectrometer; G, grating; CL, camera lens; LC, line-scan camera; FC, fiber coupler; C1 to C3, collimator; M, mirror; GB, glass block; GS, galvanometer scanner; SH, scan head. (b) Microscopic setup for phantom experiment. (c) Ophthalmoscopic setup for in vivo rodent eye imaging. L1 to L3, lenses.

JBO_22_12_121702_f003.png

4.

Experimental Results

4.1.

Phantom Validation

We imaged a flow phantom made from 0.5% intralipid (Sigma, I141-100ML) in a plastic tube with a 125-μm inner diameter. The refractive index of intralipid is 1.34.31 Flow rate was controlled by a motorized syringe pump (Model A99-EM, Razel Scientific Instruments) within a speed ranging from 0 to 0.25  μL/s. The Doppler angle of the tube was 83.5 deg. The OCT A-line rate was 20 kHz. We acquired 2-mm wide B-scans crossing the plastic tube, which consisted of 2048 A-lines. In each measurement, we performed 16 repeated B-scans for averaging. Note that in the phantom experiment, we used flow rate instead of flow velocity because flow rate, containing both velocity and cross-sectional area of flow, is of great interest in hemodynamic and metabolic studies. Using flow rate is also more convenient since the true flow rate values can be obtained from the syringe pump without the need to accurately measure the diameter of plastic tube. Therefore, we converted the flow velocity provided by Doppler OCT to flow rate in the following phantom study.

We measured flow induced phase shift using both conventional and spectroscopic Doppler analyses. In conventional Doppler analysis, we obtained the averaged phase shift inside the tube both with and without phase unwrapping. Before phase unwrapping, we further filtered the averaged phase shift images with a 3×15 (axial×lateral) pixels Doppler filter following the method reported in Ref. 22, since phase unwrapping is sensitive to noise. We selected an area close to the outer boundary of the flow area on the phase shift B-scan as a reference, where no phase wrapping occurs due to its low flow velocity. Then, we searched throughout the flow area and labeled the phase wrapped pixels if their values changed more than π compared with adjacent pixels. Finally, we added or subtracted 2π to the labeled pixels and recovered the unwrapped phase shift values.27 In spectroscopic Doppler analysis, we generated 16 Gaussian windows for STFT, whose center wavenumbers were evenly distributed from (1.06×105)rad/cm (593 nm) to (1.16×105)rad/cm (542 nm). The windows had a standard deviation of (0.02×105)rad/cm (10 nm at 565-nm wavelength). We applied STFT, corrected for possible spectroscopic phase discontinuity, and calculated subband-dependent phase shifts.

Figure 4(a) shows phase shift B-scans under two flow rates before and after phase unwrapping. There is no phase wrapping in flow 1 (0.021  μL/s) and obvious phase wrapping in flow 2 (0.218  μL/s). Therefore, before phase unwrapping, the averaged phase shift values throughout the flow area are similar between flow 1 and flow 2, which can be shown in Fig. 4(b) by the vertical positions of linear functions (data 1 and data 2). However, there is a clear difference between the slopes of the two linear functions. After phase unwrapping, the difference in flow rates between flow 1 and flow 2 can be better appreciated in Fig. 4(b) from both the vertical positions and the slopes of the linear functions (data 3 and data 4). The slopes of linear functions for data 2 and data 4 are the same, indicating that spectroscopic Doppler analysis can retrieve flow velocity without phase unwrapping.

Fig. 4

Comparing traditional and spectroscopic Doppler analyses in phantoms. (a) Averaged Doppler phase shift images of two selected flow velocities before (data 1 and data 2) and after phase unwrapping (data 3 and data 4) using traditional Doppler analysis. (b) Subband-dependent phase shifts using spectroscopic Doppler analysis.

JBO_22_12_121702_f004.png

To further compare conventional and spectroscopic Doppler analyses, we used both methods to measure seven different flow rates in phantom without phase unwrapping, as shown in Fig. 5. To increase flow measurement accuracy, we applied spatial averaging across the entire flow cross section in both methods, which provided a fair comparison. The flow rates ranged from 0 to 0.24  μL/s with a step size of 0.04  μL/s. Phase wrapping occurred in flow rates greater than 0.12  μL/s. Figure 5(a) shows subband-dependent phase shifts at different flow rates. In order to avoid crossing among different phase shift lines and allow better illustration of their slopes, we manually offset these lines according to the preset flow rates in Fig. 5(a). Therefore, the intercepts of different lines do not represent real calculated values. Figure 5(b) compares measured flow rates using both methods. Since no phase wrapping exists in the four lower flow rates, both methods can accurately retrieve the preset values. For the three higher rates, conventional Doppler OCT fails to provide accurate measurement results due to phase wrapping. However, spectroscopic Doppler OCT is still able to quantify 0.18 and 0.21  μL/s. In the case of 0.24  μL/s, a severe washout effect occurs in the center of the flow area, which decreases the SNR. Therefore, the subband-dependent phase shift array is unstable, and its fitted slope value fails to provide an accurate estimate for the flow rate.

Fig. 5

(a) Subband-dependent phase shifts for seven flow rates. Note that we manually shift these lines according to the preset flow rates to avoid crossing and allow better visualization. (b) Measured flow rates using conventional and spectroscopic Doppler analysis. The black dashed line represents the true flow rate values.

JBO_22_12_121702_f005.png

4.2.

In Vivo Retinal Blood Flow Measurements

We used both conventional and spectroscopic Doppler analyses to compare measured retinal blood flow velocity without phase unwrapping in rats (Long Evans, 300 g, Charles River) in vivo. The details of animal preparation and imaging protocol have been described elsewhere.32 All experimental procedures complied with the Association for Research in Vision and Ophthalmology statement for the use of animals in ophthalmic and vision research. The laboratory animal study protocol was approved by the Institutional Animal Care and Use Committee at Northwestern University.

In brief, vis-OCT measurements were conducted in rat eyes using a circular scanning protocol.33,34 The diameter of the scanning circle was 0.5  mm on the rat retina. Each circular scan contained 4096 A-lines. The vis-OCT probing power was 0.8 mW on the cornea. We first used conventional Doppler analysis, where we calculated phase shift images using the whole broadband spectrum. Then, we randomly selected 24 vessels whose phase shifts were roughly distributed between 0.2 and 1 rad. We applied spectroscopic analysis to these selected vessels, where we used the same set of Gaussian windows as in the phantom experiment, but adjusted the standard deviation to 0.03×105  rad/cm (around 15 nm at 565-nm wavelength). We applied STFT on the raw OCT interferograms and obtained a series of phase shift images. We manually segmented the flow areas within retinal blood vessels from subband phase shift images. We obtained axial flow velocities, converted them to phase shift values, and compared them with the results from the conventional Doppler method. There are two reasons why we did not provide flow rates as in the phantom study. First, in the phantom study, since the true flow rate is continuously controllable, it is convenient to compare between the preset flow rate values and the measured values. Second, in the animal experiment, the flow rates across different vessels may be similar. However, variations in the Doppler angles among different blood vessels give rise to distinct axial flow velocities. Therefore, by using axial flow velocity without calculating the flow rate, we can compare the performance of conventional and spectroscopic Doppler methods at a wide range of phase shift values.

Figure 6 compares the rat retinal axial flow velocities represented by Doppler phase shifts using conventional and spectroscopic Doppler analyses without phase unwrapping. The scatterplot shows a good correlation between phase shifts measured by two methods when the mean Doppler phase shift is smaller than 1 rad. However, there is a larger variation as compared with the phantom results shown in Fig. 5 due to a reduced SNR. We can also observe in Fig. 6 that the five highlighted data points all lie below the diagonal line, suggesting lower flow velocity estimation by the conventional Doppler analysis. This is consistent with the phantom experiment that indicates that the conventional Doppler OCT method tends to underestimate flow velocity when phase wrapping occurs.

Fig. 6

Axial retinal blood flow velocities measured by conventional and spectroscopic Doppler analyses, both converted to radian values for comparison. The dashed line suggests identical values measured by both methods. The dashed ellipse highlights five cases where phase wrapping possibly occur. The two insets are cross-sectional phase images of two selected blood vessels.

JBO_22_12_121702_f006.png

5.

Discussion and Conclusion

We developed a spectroscopic Doppler analysis for flow velocity measurement. Instead of calculating Doppler phase shift using broadband OCT interferogram as in conventional Doppler OCT, we used STFT to reconstruct a series of complex OCT B-scans in different subbands, whose Doppler phase shifts were then calculated separately. Instead of averaging the phase shift values for improved SNR as proposed previously by Tan et al.,35 we fit these phase shift values against the center wavenumbers of the subbands using a linear function whose slope was related to the axial flow velocity. We demonstrated that spectroscopic Doppler analysis was unaffected by phase wrapping since the absolute shift values were not directly used in velocity calculation. Therefore, spectroscopic Doppler analysis does not rely on phase unwrapping to measure a vascular flow whose center velocity exceeds PLV. Using subbands with reduced bandwidth will decrease the axial resolution in the Doppler phase shift images, which affects measurement accuracy if the flow area cannot be clearly resolved at the reduced axial resolution.

Compared with conventional Doppler OCT and its phase unwrapping techniques, spectroscopic Doppler analysis may be a good alternative for flow velocity measurement in vis-OCT. Working on shorter wavelengths, vis-OCT has much improved axial resolution, but is more susceptible to phase wrapping. Due to higher resolution, the phase shift images of vis-OCT retain sufficient information even when only reconstructed on a spectral subband. Due to a much lower PLV, vis-OCT requires more careful phase unwrapping procedures when measuring flow velocity than using a conventional Doppler OCT method. Though phase unwrapping can provide satisfactory results when the image SNR is guaranteed, as shown in our phantom experiment, it is not trivial, and is even biased if the image quality is moderate, especially in clinical settings. In addition, the unwrapping errors may accumulate during the iterative process adopted in most unwrapping algorithms.2728.29 Therefore, spectroscopic analysis is especially helpful in vis-OCT, though it can also be applied to conventional NIR OCT as well. However, we expect that spectroscopic analysis may be less valuable for OCT systems working beyond 1  μm because (1) these OCT systems have a much higher PLV and can depend on a conventional Doppler method and (2) the phase shift values do not vary much with the center wavenumber of different subbands and the fitted slope value may be more vulnerable to noise.

The total number and bandwidth of the spectral subbands affect the flow measurement accuracy. A limited number of subbands leads to fewer points for linear regression, which makes the fitted slope susceptible to noise. Too many subbands, on the other hand, requires more intensive computation with a marginal benefit. Therefore, we selected 16 subbands to perform spectroscopic analysis, which was a tradeoff between accuracy and computation efforts. We further confirmed our selection using numerical simulation. Selecting proper bandwidth for subbands is also important. A large bandwidth decreases the spectral fitting range and a small bandwidth leads not only to reduced resolution, but also lower SNR.36 In our experiment, we empirically used (0.04×105)rad/cm for phantom experiments and (0.06×105)rad/cm for animal experiments. (Note that the subband width is twice the standard deviation of the Gaussian window used in STFT.) We used a larger bandwidth in animal experiments to limit the phase noise,36 since the SNR of in vivo data is naturally lower than that of the phantom data. A more thorough investigation will be necessary to provide comprehensive evaluation of the effects of subband selection for spectroscopic Doppler analysis.

Though spectroscopic Doppler analysis is robust to phase wrapping, it is still limited by the fringe washout effect in spectrometer-based SD-OCT. The subband-dependent phase shift becomes noisy when significant washout occurs and the slope extracted from linear fitting fails to reveal the flow velocity. SS-OCT is almost unaffected by the fringe washout effect since it increases the washout threshold by a factor that equals the total number of spectral samples in each A-line.19 Therefore, although not tested here, our spectroscopic Doppler analysis can potentially extend the velocity measurement range in SS-OCTs.

We notice that in the current study, sufficient SNR is critical to obtain reliable flow velocity using spectroscopic Doppler analysis. To guarantee sufficient SNR, we applied spatial averaging across the entire flow cross section in both phantom and animal studies. Since the flow velocity is translated from the slope of the fitted straight line, accurate measurement requires consistency of phase shift values calculated from spectral subbands, which themselves could be more susceptible to influence from noise due to reduced spectral bandwidth. Therefore, our future work will include comprehensively investigating the validity of spectroscopic Doppler analysis under different SNR conditions.

In summary, we present a spectroscopic Doppler analysis for flow velocity measurement in vis-OCT and tested it in both phantom and in vivo experiments. This method uses spectroscopic information of each phase-shifted image pixel to unambiguously retrieve flow velocity instead of depending solely upon the pixel spatial relationship to correct the phase wrapping as in the conventional Doppler OCT method. Though we demonstrated our method using spectrometer-based OCT only, we expect that SS-OCT will benefit from this method as well.

Disclosures

Hao F. Zhang, and Wenzhong Liu have financial interests in Optician Health, which did not support this work.

Acknowledgments

We acknowledge the generous financial support from the National Institutes of Health (NIH) Grant Nos. DP3DK108248, R24EY022883, R01EY026078, and R21EY027502.

References

1. W. Sickle, “Retinal metabolism in dark and light,” in Physiology of Photoreceptor Organs, pp. 667–727, Springer, Berlin, Heidelberg (1972). Google Scholar

2. R. A. Kowluru, J. Tang and T. S. Kern, “Abnormalities of retinal metabolism in diabetes and experimental galactosemia VII. Effect of long-term administration of antioxidants on the development of retinopathy,” Diabetes 50(8), 1938–1942 (2001).DIAEAZ0012-1797 http://dx.doi.org/10.2337/diabetes.50.8.1938 Google Scholar

3. R. A. Linsenmeier and H. F. Zhang, “Retinal oxygen: from animals to humans,” Prog. Retinal Eye Res. 58, 115–151 (2017). http://dx.doi.org/10.1016/j.preteyeres.2017.01.003 Google Scholar

4. T. S. Kern, R. A. Kowluru and R. L. Engerman, “Abnormalities of retinal metabolism in diabetes or galactosemia: ATPases and glutathione,” Invest. Ophthalmol. Visual Sci. 35(7), 2962–2967 (1994).IOVSDA0146-0404 Google Scholar

5. R. L. Engerman and T. S. Kern, “Hyperglycemia as a cause of diabetic retinopathy,” Metabolism 35(4), 20–23 (1986).METAAJ0026-0495 http://dx.doi.org/10.1016/0026-0495(86)90182-4 Google Scholar

6. R. A. Kowluru and P.-S. Chan, “Oxidative stress and diabetic retinopathy,” Exp. Diabetes Res. 2007, 43603 (2007). http://dx.doi.org/10.1155/2007/43603 Google Scholar

7. S. H. Hardarson and E. Stefánsson, “Retinal oxygen saturation is altered in diabetic retinopathy,” Br. J. Ophthalmol. 96(4), 560–563 (2012).BJOPAL0007-1161 http://dx.doi.org/10.1136/bjophthalmol-2011-300640 Google Scholar

8. W. Liu et al., “Increased retinal oxygen metabolism precedes microvascular alterations in type 1 diabetic miceretinal oxygen consumption increases in diabetic mice,” Invest. Ophthalmol. Visual Sci. 58(2), 981–989 (2017).IOVSDA0146-0404 http://dx.doi.org/10.1167/iovs.16-20600 Google Scholar

9. J. Yi et al., “Visible light optical coherence tomography measures retinal oxygen metabolic response to systemic oxygenation,” Light Sci. Appl. 4(9), e334 (2015). http://dx.doi.org/10.1038/lsa.2015.107 Google Scholar

10. S. Chen, J. Yi and H. F. Zhang, “Measuring oxygen saturation in retinal and choroidal circulations in rats using visible light optical coherence tomography angiography,” Biomed. Opt. Express 6(8), 2840–2853 (2015).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.6.002840 Google Scholar

11. S. Chen et al., “Retinal oximetry in humans using visible-light optical coherence tomography [Invited],” Biomed. Opt. Express 8(3), 1415–1429 (2017).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.8.001415 Google Scholar

12. J. Yi and X. Li, “Estimation of oxygen saturation from erythrocytes by high-resolution spectroscopic optical coherence tomography,” Opt. Lett. 35(12), 2094–2096 (2010).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.35.002094 Google Scholar

13. S. Chen et al., “Dual-band optical coherence tomography using a single supercontinuum laser source,” J. Biomed. Opt. 21(6), 066013 (2016).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.21.6.066013 Google Scholar

14. S. Chen et al., “Imaging hemodynamic response after ischemic stroke in mouse cortex using visible-light optical coherence tomography,” Biomed. Opt. Express 7(9), 3377–3389 (2016).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.7.003377 Google Scholar

15. J. A. Izatt et al., “In vivo bidirectional color Doppler flow imaging of picoliter blood volumes using optical coherence tomography,” Opt. Lett. 22(18), 1439–1441 (1997).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.22.001439 Google Scholar

16. S. Yazdanfar, M. Kulkarni and J. Izatt, “High resolution imaging of in vivo cardiac dynamics using color Doppler optical coherence tomography,” Opt. Express 1(13), 424–431 (1997).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.1.000424 Google Scholar

17. Y. Zhao et al., “Phase-resolved optical coherence tomography and optical Doppler tomography for imaging blood flow in human skin with fast scanning speed and high velocity sensitivity,” Opt. Lett. 25(2), 114–116 (2000).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.25.000114 Google Scholar

18. Z. Chen et al., “Optical Doppler tomographic imaging of fluid flow velocity in highly scattering media,” Opt. Lett. 22(1), 64–66 (1997).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.22.000064 Google Scholar

19. H. C. Hendargo et al., “Doppler velocity detection limitations in spectrometer-based versus swept-source optical coherence tomography,” Biomed. Opt. Express 2(8), 2175–2188 (2011).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.2.002175 Google Scholar

20. Y. Wang et al., “In vivo total retinal blood flow measurement by Fourier domain Doppler optical coherence tomography,” J. Biomed. Opt. 12(4), 041215 (2007).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2772871 Google Scholar

21. R. A. Leitgeb et al., “Doppler optical coherence tomography,” Prog. Retinal Eye Res. 41(100), 26–43 (2014). http://dx.doi.org/10.1016/j.preteyeres.2014.03.004 Google Scholar

22. D. J. Faber et al., “Light absorption of (oxy-) hemoglobin assessed by spectroscopic optical coherence tomography,” Opt. Lett. 28(16), 1436–1438 (2003).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.28.001436 Google Scholar

23. C. Xu et al., “Near-infrared dyes as contrast-enhancing agents for spectroscopic optical coherence tomography,” Opt. Lett. 29(14), 1647–1649 (2004).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.29.001647 Google Scholar

24. J. Yi and V. Backman, “Imaging a full set of optical scattering properties of biological tissue by inverse spectroscopic optical coherence tomography,” Opt. Lett. 37(21), 4443–4445 (2012).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.37.004443 Google Scholar

25. F. E. Robles et al., “Molecular imaging true-colour spectroscopic optical coherence tomography,” Nat. Photonics 5(12), 744–747 (2011).NPAHBY1749-4885 http://dx.doi.org/10.1038/nphoton.2011.257 Google Scholar

26. A. L. Oldenburg, C. Xu and S. A. Boppart, “Spectroscopic optical coherence tomography and microscopy,” IEEE J. Sel. Top. Quantum Electron. 13(6), 1629–1640 (2007).IJSQEN1077-260X http://dx.doi.org/10.1109/JSTQE.2007.910292 Google Scholar

27. D. C. Ghiglia and L. A. Romero, “Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods,” J. Opt. Soc. Am. A 11(1), 107–117 (1994). http://dx.doi.org/10.1364/JOSAA.11.000107 Google Scholar

28. J. M. Huntley and H. O. Saldner, “Error-reduction methods for shape measurement by temporal phase unwrapping,” J. Opt. Soc. Am. A 14(12), 3188–3196 (1997).JOAOD60740-3232 http://dx.doi.org/10.1364/JOSAA.14.003188 Google Scholar

29. R. M. Goldstein, H. A. Zebker and C. L. Werner, “Satellite radar interferometry: two-dimensional phase unwrapping,” Radio Sci. 23(4), 713–720 (1988). http://dx.doi.org/10.1029/RS023i004p00713 Google Scholar

30. J. Yi et al., “Human retinal imaging using visible-light optical coherence tomography guided by scanning laser ophthalmoscopy,” Biomed. Opt. Express 6(10), 3701–3713 (2015).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.6.003701 Google Scholar

31. H. Ding et al., “Determination of refractive indices of porcine skin tissues and intralipid at eight wavelengths between 325 and 1557 nm,” J. Opt. Soc. Am. A 22(6), 1151–1157 (2005).JOAOD60740-3232 http://dx.doi.org/10.1364/JOSAA.22.001151 Google Scholar

32. W. Liu et al., “Measuring retinal blood flow in rats using Doppler optical coherence tomography without knowing eyeball axial length,” Med. Phys. 42(9), 5356–5362 (2015).MPHYA60094-2405 http://dx.doi.org/10.1118/1.4928597 Google Scholar

33. H. Wehbe et al., “Automatic retinal blood flow calculation using spectral domain optical coherence tomography,” Opt. Express 15(23), 15193–15206 (2007).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.15.015193 Google Scholar

34. Y. Wang et al., “Retinal blood flow measurement by circumpapillary Fourier domain Doppler optical coherence tomography,” J. Biomed. Opt. 13(6), 064003 (2008).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2998480 Google Scholar

35. O. Tan et al., “En face Doppler total retinal blood flow measurement with 70 kHz spectral optical coherence tomography,” J. Biomed. Opt. 20(6), 066004 (2015).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.20.6.066004 Google Scholar

36. H. C. Hendargo et al., “Synthetic wavelength based phase unwrapping in spectral domain optical coherence tomography,” Opt. Express 17(7), 5039–5051 (2009).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.17.005039 Google Scholar

Biographies for the authors are not available.

© 2017 Society of Photo-Optical Instrumentation Engineers (SPIE)
Xiao Shu, Wenzhong Liu, Lian Duan, Hao F. Zhang, "Spectroscopic Doppler analysis for visible-light optical coherence tomography," Journal of Biomedical Optics 22(12), 121702 (17 October 2017). http://dx.doi.org/10.1117/1.JBO.22.12.121702 Submission: Received 3 May 2017; Accepted 18 August 2017
Submission: Received 3 May 2017; Accepted 18 August 2017
JOURNAL ARTICLE
8 PAGES


SHARE
Back to Top