Wide velocity range Doppler optical microangiography using optimized step-scanning protocol with phase variance mask

Abstract. We propose a simple and optimized method for acquiring a wide velocity range of blood flow using Doppler optical microangiography. After characterizing the behavior of the scanner in the fast scan axis, a step-scanning protocol is developed by utilizing repeated A-scans at each step. Multiple velocity range images are obtained by the high-pass filtering and Doppler processing of complex signals between A-scans within each step with different time intervals. A phase variance mask is then employed to segment meaningful Doppler flow signals from noisy phase background. The technique is demonstrated by imaging in vivo mouse brain with skull left intact to provide bidirectional images of cerebral blood flow with high quality and wide velocity range.


Introduction
Doppler optical microangiography (DOMAG) 1 is a technological extension to Doppler optical coherence tomography (DOCT), 2 in which DOCT is combined with recently developed OMAG technique 3 to provide velocity mapping of blood flows within tissue beds in vivo. Due to its simple system and algorithm implementation, phase-resolved technique 2 is often used to deduce the axial flow velocity, V z ¼ λφ∕ð4πnTÞ, where λ is the center wavelength of the light source, n is the tissue refractive index, and ϕ and T are the phase difference and time interval between adjacent A-scans, respectively. The maximal and minimal velocities that are detectable in one B-scan [or three-dimensional (3-D) scan] are determined by the time interval T, the πambiguity and the system phase noise level. 2 The shorter T, i.e., faster imaging speed, is usually required to map the velocities of faster flow, but at an expense of slower ones and vice versa. Therefore, it is currently difficult for DOMAG/DOCT to provide the read-out of blood flows within tissue beds that have wide velocity distributions. However, capability to simultaneously capture blood flow velocities at different ranges from arteries and veins down to capillary vessels is important for comprehensive study of vascular responses in a number of diseases that have vascular components, for example, in the investigation of cerebral blood flow (CBF) responses to a neurological insult, e.g., stroke. 4,5 This is because flow velocity within cerebral cortex exhibits a wide range of values, from tens of micrometer/seconds in capillaries to tens of millimeter/seconds in middle cerebral artery (MCA).
We notice that there are a number of methods previously reported in OCT angiography 6-8 that were tuned to image the fast or slow flows. Most of these methods employ a linear scan (L-scan) in the fast axis [ Fig. 1(a)]. They are, however, not amenable to the requirement of concurrently obtaining the fast and slow velocities provided that just one 3-D scan is performed. By skipping A-scans for DOCT processing, it is possible to obtain multiple velocity ranges within one scan. 6 However, the tissue of interest is often highly scattering and heterogeneous, leading to decreased correlation degree (CD) between A-scans. 9 Note that the high degree of correlation is an utmost requirement for Doppler algorithm to accurately evaluate flow velocity. To increase the degree of correlation, normally the B-scan has to be densely sampled, which, however, translates to reducing the frame rate of the system, not suitable for in vivo imaging. The situation is even worse in the case of using skipping A-scan approach to obtain multiple velocity ranges. To mitigate this problem, Meemon and Rolland 10 used a step-scan strategy in the fast axis to obtain multiple velocity ranges, in which the scanner is driven by a stepped waveform while M-mode scan (i.e., repeated A-scan) is employed in each step. The assumption here is that the M-mode scan in each step provides 100% correlation between any two A-scans that are used to deduce the flow velocity. This is true in an ideal situation. However, in practice when the scanner switches from one step to next, it requires some time for the scanner to settle due to its mechanical inertia. Thus, the correlation between any two Ascans in each step is highly dependent upon the mechanical stability of the scanner used. In addition, due to the communication between the computer, the scanner, and the data acquisition device, this settling time is also sometimes system dependent. Consequently, the scanner behavior in the setting of entire system has to be systematically characterized so that an improved step-scan strategy can be devised to optimize the system capability to provide the imaging of wide velocity range. To meet with this requirement, this article proposes an alternative method that is able to provide multirange DOMAG (mDOMAG) with high signal quality based on a common spectral domain OCT (SD-OCT) system setup. Below, we first demonstrate the necessity for characterizing the scanner behavior during scanning for step-scan (S-scan) strategy to provide a multirange velocity map. We then propose an optimized method in which S-scan is employed in the fast axis while a number of repeated A-scans at each step are acquired only after the scanner has been stabilized.

System Setup and Methods
We employed a fiber-based SD-OCT system to demonstrate the necessity and utility of the proposed optimized step-scan strategy. The SD-OCT system was similar to that previously described. 1 The center wavelength of the light source was λ ¼ 1310 nm and both axial and lateral resolutions were ∼10 μm. The detection unit was a fast spectrometer in which a line InGaAs camera was used to capture the interferograms formed between the lights that were back reflected from the reference and sample arms. The x-y scanner was used in the sample arm to scan a focused beam spot over the sample to acquire volumetric dataset, representing tissue morphology and microcirculations. The system had a measured dynamic range of 105 dB with the light power of 5 mW at sample surface. For experiments reported below, the imaging speed of the system used was 75 kHz A-scan rate.
We first performed experiments to characterize the behavior of the scanner during M-mode scan, upon which to design the optimized scanning pattern. In experiments, we used a stationary cardboard as the sample. In this case, the complication for in vivo tissue such as tissue motion or cell swelling/shrink would be absent, facilitating the use of a correlation approach to characterize the OCT signals to achieve our purpose. The fast axis scanner was driven by a step waveform with 400 steps for one B-scan. In each step, N ¼ 50 A-scans were acquired, thus the total number of A-scans in each B-scan was 20,000. In the slow axis direction, we captured 200 B-frames to accomplish one 3-D scan covering an area of 2 × 2 mm 2 . For each B-scan, the absolute values of CD between the i'th and (i þ k)'th A-scans were calculated, where k is the number of skipped A-scans. In this study, we selected k ¼ 1, 3, and 9 to provide three ranges of velocities, corresponding to AE18.2, AE6.1, and AE2 mm∕s, respectively, assuming that the refractive index of tissue is 1.35. Finally, CDs were averaged throughout the 3-D scan according to modði þ k∕2; NÞ, which was the remainder of (i þ k∕2) divided by N. The results are given in Fig. 1(b). We found that for the system used for this study, the minimal CD corresponds to modði þ k∕2; NÞ ¼ 37.5 for any k, indicating the time point within each step when the scanner is most unstable. This is an important finding because it is contrary to the expected situation where we would expect the most unstable region is close to the beginning of the step edge. The unstable region being shifted to the second half of the step is most likely caused by the communication between the scanner, the computer, and the detection device for synchronized data acquisition. According to this finding, the system should only utilize the first half of the step to provide improved imaging of multirange velocities. Therefore, we designed the camera trigger to capture only the A-scans when modði; NÞ < 25 is satisfied. In this case, when skipping k A-scans within this set of A-scans, they would always meet the condition defined by modði þ k∕2; NÞ < 25 [see bottom of Fig. 1(b), left region of vertical line]. Note that the scanner behavior is system dependent; thus, it is recommended that each system has to be individually optimized in order to use the step-scan to obtain variable velocity ranges.
Using the above optimized S-scan, we imaged a mouse brain in vivo with skull left intact. Figure 2(a) shows a typical B-scan DOCT image with k ¼ 1 using the mature Kasai algorithm. 11 To enhance the sensitivity, a high-pass filter was used prior to Doppler algorithm, assuming a linear superposition of stationary and moving scatters. 12 We employed a differential filter reported in ultrahigh sensitive OMAG, 13 which subtracts the complex OCT signals between adjacent A-scans to eliminate the static scattering signals (1) where i is the index of the A-scan number and SðiÞ is the complex signal of the i'th A-scan. The Kasai method was then used to obtain the phase difference ϕ by skipping k A-scans: φðkÞ ¼ ∠f P i ½Fði þ kÞF Ã ðiÞg, where * denotes the complex conjugate. The result is shown in Fig. 2(b) where the nonflow area contains significant phase noise, resulted from the subtraction operation of Eq. (1). To extract meaningful Doppler signal, phase variance 14,15 is used to eliminate the noisy background, including phasewrapped signals. Specifically, we calculated the phase variance  of the filtered signal: PðkÞ ¼ σ 2 i f∠½Fði þ kÞF Ã ðiÞg to define a binary mask with a constant threshold RðkÞ. The values of RðkÞ are determined by the phase stability of the system, the flow variation status in the tissue, and the number of A-scans. With the current system setup, we determined RðkÞ ¼ 2.4, 1.4, and 0.4 for k ¼ 1, 3, and 9, respectively. Thus, the final mDOMAG velocity image is given by This operation has an additional advantage of eliminating fast flows that are beyond the capability of the system at a specific range of interest, causing phase wrapping in the phase signals. The final results are shown in Figs. 2(c)-2(e) for k ¼ 1, 3, and 9, respectively, where the Doppler signals for three velocity ranges were successfully extracted. The vessels with the fast flow are imaged in the k ¼ 1 range [arrows in Fig. 2(c)], while these vessels are disappeared in the k ¼ 3 range [Fig. 2(d)] due to their phase wrapping. The vessels with intermediate velocity (around 6 mm∕s) start to emerge in the k ¼ 3 range [arrows in Fig. 2(d)]. In the slowest range of k ¼ 9, numerous capillaries signals were seen [arrows in Fig. 2(e)], which are, however, beyond the detection capability of other ranges.

Results
Finally, we demonstrate the utility of the mDOMAG in the imaging of 3-D CBF within a mouse brain in vivo. In this experiment, we imaged 2 × 2 tiles of 2 × 2 mm 2 areas of the mouse brain to obtain CBF mapping within an area of 3.6 × 3.6 mm 2 . The mosaic bidirectional en face maximum projection images for k ¼ 1, 3, and 9 are given in Figs. 3(a), 3(b), and 3(c), showing the functional vessels in three velocity ranges of AE18.2, AE6.1, and AE2 mm∕s, respectively. At k ¼ 1, the velocity image mainly provides the vessels with the fast flows in the branches of MCA [pointed by arrows in Fig. 3(a)]. These branched arterioles are mostly disappeared in the velocity range of k ¼ 3 [ Fig. 3(b)] due to their phase wrapping effect causing large phase variance, which is subsequently masked out by Eq. (2). The diving arterioles and rising venules appear as red and green isolated spots in Fig. 3(a) (e.g., marked by circles). This is because their flow directions are mainly parallel to the beam axis, giving their axial velocity large enough to be detected at this range. However, the surface vessels connected to them flow nearly perpendicularly to the probe beam, escaped the detection at k ¼ 1 but fell within the range for k ¼ 3 (e.g., hollow arrow heads). The venules that run almost parallel to the branched arterioles of the MCA over the surface of cortex are detected in Fig. 3(b) (pointed by the arrows), which have a velocity slower than the arterioles. At the smallest range [ Fig. 3(c)], phase variance threshold is sufficient to keep all the slow flow signals, including capillaries (e.g., pointed by arrow heads). However, some out-of-range flow could not be completely removed, which is observed as phase-wrapped signals, seen as yellow color (combination of red and green). From the images with three different velocity ranges, we would be able to restore the complete vasculature. To restore, we combined the three 3-D mDOMAG maps into a single map with wide velocity range using a simple criterion: mDOMAG ¼ Vðk C Þ, where k C is the argument of the maximum of jVðkÞj. Figures 3(d)-3(f) show the combined data displayed with the same velocity ranges as in Figs. 3(a)-3(c), respectively. The vessel continuity is much improved when compared with the uncombined data. The slowest velocity range gives the best connectivity of the vessels, but the phase noise from the larger range data is amplified, resulting in some vessels appearing yellowish [arrow head in Fig. 3(f)]. This suggests that the monitoring of CBF within the cortex requires the investigation of both original and combined data. Overall, mDOMAG presents results with low noise, giving the maximal projection images with bidirectional multirange flow velocities and with high quality.
Note that flow velocity obtained after high-pass filtering deviates from the expected parabolic profile (Figs. 2 and 3). This phenomenon results from the overestimation of small Doppler signals affected by the differential filter. 12 Decreasing the cut-off frequency of the filter could lower overestimation, but at an expense of detecting more low-frequency components, which might instead cause underestimation. Also the assumption of linear contributions from different frequency components might not be optimal, leading to a biased Doppler estimation under high-pass filtering. Our future work will focus on compensating these errors or designing a new high-pass filter with less influence on the Doppler estimation. Also note that in the current study, we did not consider the relation between the degree of correlation and the fidelity of the Doppler signal 16 when imaging the CBF. Such relation can be important for in vivo applications, which might be useful to guide the optimization of scanning patterns to solve the limitations of the current approach.

Conclusion
In conclusion, we have demonstrated an optimized multirange DOMAG in which an S-scan protocol was used to achieve ∼100% correlation between repeated A-scans. High sensitivity and high quality bidirectional mDOMAG images with a wide velocity range were shown in mouse brain in vivo with skull left intact. It is recommended that to perform multiple velocity imaging using skipping, A-scan approach requires the characterization of the scanner in each individual system. Fig. 3 Bidirectional projection mDOMAG of a mouse brain in vivo. (a)-(c) The separately processed images at k ¼ 1, 3, and 9, respectively. (d)-(f) The combined velocity map but displayed with the same color bars as in (a)-(c), respectively. The imaging area is located between bregma and lambda at the right brain. A, anterior; L, lateral. White box ¼ 0.5 × 0.5 mm 2 .