13 February 2015 One specific velocity color mapping using optical coherence tomography
Author Affiliations +
J. of Biomedical Optics, 20(5), 051034 (2015). doi:10.1117/1.JBO.20.5.051034
Depth resolved coherence gating along with Doppler shift detection of the carrier frequency is used for one predetermined velocity mapping in different flows. Bidirectional rapid scanning optical delay of optical coherence tomography system is applied in the reference arm. Tilted capillary entry is used as a hydrodynamic phantom to model a sign-variable flow with complex geometry. Structural and one specific velocity images are obtained from the scanning interferometer signal processing in the frequency domain using analog and digital filtering. A standard structural image is decomposed into three parts: stationary object, and positive and negative velocity distributions. The latter two show equivelocity maps of the flow. The final image is represented as the complexation of the three.
Proskurin, Potlov, and Frolov: One specific velocity color mapping using optical coherence tomography



Direction sensitive investigation of reversible flows using phase sensitive technique was applied in acoustic reflectometry1 and optical coherence tomography (OCT).2. Velocity image reconstruction based on estimation of phase change between sequential A-scans is represented as the color or Doppler image in both. For that purpose, cross-3 and auto-correlation45. functions of the reflected signal are used. For example, the OCT system presented in Ref. 3 is based on the calculation of cross-correlation function for a certain number, N, of adjacent A-scans. The mean phase shift variation is used as


where «Arg» and «*»—denote the operation of taking argument of a complex signal and complex conjugation, respectively. Si,T and Si+1,T* are the complex signals of adjacent A-scans. T is the temporal interval between the adjacent A-scans. The Doppler frequency is then represented as


where fA=(1/T) is the reference arm scanning frequency. Flow velocity, V, is then estimated as follows:


where λo is the central wavelength of the low coherence source, n is the refractive index of the media, and α is the Doppler angle.

The OCT system presented in Ref. 5 operates in the same way. The difference is that the cross-correlation function is used5,6 in place of the autocorrelation Kasai function.1 Flow velocity is then calculated from the phase shift of the consecutive A-scans. This allows obtaining high sensitivity to the motion of sharp surfaces or sparse particles in a flow, but leads to the 2π-ambiguity.

To reduce sensitivity to the changes of the incidence angle and to the velocity pulsations of the moving blood, Doppler spectra variance, σ2, can be used as8,9


In Refs. 10 and 11, improved Doppler optical microangiography is described. In this method, the back reflected signal only from the moving particles is registered. This leads to a considerable reduction of the noise to compare with the standard backscattered signal processing. Blood flow velocity is determined from the phase shift of the signal carrier frequency.

A double correlation technique is notable for motion artifacts’ reduction and background noise removal,12 for which adaptive Wiener filtering is applied. It is also used to increase the structural image quality.

High sensitivity biological fluids’ velocity measurement13,14 is based on the phase synchronization of the reference signal, ΔΦR, and in the sample, ΔΦS. An electro-optic phase modulator with periodic phase tuning is used in the reference arm.14 In practice, this time domain system allows the reduction of motion artifacts which causes interference pattern blurring, that, in turn, causes the sensitivity reduction of the flow registration. This approach is supposed to quantitatively register velocity change, but the Doppler frequency is, again, represented using the phase shift

where k0 is the central wave number and VR is the reference mirror velocity, which is defined as
where ΔZR is the reference mirror shift for the time, τ. To quantify ΔZR, it is necessary to determine electrical values
where β is the constant part of the phase, ΔV is the change of voltage for τ, and Vπ is the voltage at which the central wave number14 is shifted by π.

Another approach to obtain velocity information is based on the speckle flickering frequency in the structural OCT images.15 Speckle fluctuations depend on the average flow velocity. Therefore, the ratio of the high and low components of the intensity fluctuation power spectrum could show the average velocity of a flow. The Doppler shift is then calculated using the following equation:16,17

where iAC2 is the average value of the alternating photocurrent, iDC2 is the average value of the direct photocurrent, and NS is the number of the detected speckles. And, it is assumed, the detector’s surface should be uniformly illuminated.

The description method for one specific velocity (OSV) mapping is based on the registration of the Doppler shift of the carrier, fc, itself. This gives the average integrated information over the probe volume and does not depend on Doppler spectrum broadening. The latter does not give the absolute value of the flow velocity, but primarily has a statistical relationship with it rather than a functional one. In addition, OSV mapping does not have the 2π ambiguity disadvantage. The signal processing of Doppler spectra themselves is performed in the frequency domain after a short time Fourier transform (STFT).

Since it is difficult to register the Doppler shift of the carrier itself when applying ultrafast scanning techniques, fD=±2Vcos(α)/λ0, we use the modified double pass rapid scanning optical delay (RSOD)18,19 at lower scanning frequencies. RSOD allows the decoupling phase and group delay, as well as a shift of fc to the values of 20 to 30 kHz. The maximum signal-to-noise ratio of the detected photocurrent appears to be in the frequency range between the decreasing 1/f noise and the increasing white noise. Using a time-domain optical delay and processing the signal in the frequency domain, it is possible to decompose a standard motionless structural image, and OSV images of a chosen predetermined velocity, V, with a certain predetermined accuracy, ΔV.

In the current work, the method of one predetermined velocity mapping is applied to study a flow with complex geometry in a tilted capillary entry (tilted die entry)20 at a steeper angle, α=77deg. The initial OCT image is decomposed into a structural image and two velocity images corresponding to the positive and negative directions of the flow. The final image is obtained after color coding and complexation of the three.


Materials and Methods

The experimental setup (Fig. 1) is a scanning fiber-based Michelson interferometer. Superluminescent diode (SLD) irradiation (λo=1300nm, Δλ=70nm) after the fiber coupler goes to the reference arm—an optical delay line (modified RSOD19). Another part of the irradiation goes to the sample arm for the investigated flow. The axial resolution of the system is about 11μm. Tilted capillary entry is used to model a flow with complex geometry (Fig. 2). On the way back, irradiation is focused on the balanced detectors. After an amplifier and a tunable band-pass filter (2 to 200 kHz), the electric signal goes to a 10-bit analog-to-digital converter. Final processing is performed in the frequency domain using digital STFT.

Fig. 1

Schematic representation of the experimental setup—scanning Michelson interferometer. SLD—broadband source, superluminescent diode; D—detectors in the balance scheme; ODL—modified rapid scanning optical delay line.


Fig. 2

Tilted capillary entry, α=77deg.


For α=77deg, the speed range in the hydrodynamic phantom is equal to 0 to 28mm/s. A 7.5-kHz Doppler shift corresponds to a velocity of V=17mm/sec. The dynamic range of the system is about 100 dB. The system has been tested on the parabolic velocity profile in a cylindrical capillary with Doppler spectrum broadening Δf0, 3–1, 6-kHz depending on the frequency of the centroid, fo


The standard deviation of the parabolic profile ΔV/V is about 5% to 7%. A hydrodynamic phantom has been made of a transparent plastic material with two different inner diameters d1=2.2mm, d2=0.55mm. This gives a complex flow with a converging ratio of 41. The outer diameter of the phantom is 5 mm. This allows the study of the flow upstream (before the entry) and downstream (after the entry), with a resolution which cannot be detected by Doppler ultrasound systems. A 1% intralipid solution in distilled water was used for obtaining the backscattered signal and velocity images of the flow. The flow rate was equal to 0.11ml/min, the average velocity of the fluid in the large lumen (d1) was about 0.48mm/s and, assuming a parabolic velocity profile, the axial flow velocity was V1max=0.97mm/s. In the lumen with the smaller diameter (d2), the average flow rate was about 7.7ml/s, and V2max=15.4mm/s. Since the percentage of the added intralipid was low, it was assumed that the density ρ and the dynamic viscosity of the solution were similar to those in water at 20°C.

The applied optical delay allows shifting the carrier to the region with minimal noise with a sufficiently fast repetition rate of A-scans (several kilohertz). This allows obtaining good enough structural images and registers the Doppler shift of the carrier corresponding to the indicated velocities. Odd A-scans always have a positive Doppler shift, while in the even A-scans, it is always negative. This approach gives the possibility of splitting the initial image into the three images corresponding to the steady object and images showing the positive and negative directions of the flow (positive and negative OSV images).


Velocity Mapping Algorithm

Using the aforementioned approach, one specific velocity mapping algorithm was suggested. It includes: (1) separation of the input data into two parts which correspond to the positive and negative shifts of the carrier frequency; (2) the independent construction of a two-dimensional (2-D) structural OCT image and (3) 2-D OSV images with (4) consecutive color coding.

x[s] is the input set of data, where s=1,2,,S and S is the number of samples in the signal. It is received from the raw data file and represents the amplified and digitized mean intensity of the radiation detected by the system.

Then it is separated to the functions xodd[p,k] and xeven[p,k] corresponding to the positive and the negative shifts of the carrier frequency as follows:

where p=1,2,,P, p is the count number within the A-scan, k=S/p, and P is the number of samples per one A-scan.

The resulting signal xsum[p,k] is the algebraic sum of xodd[p,k] and xeven[p,k]:


Note that from the xsum[p,k] set the analog band-pass filter of the system separates the high- and low-parts of the spectrum.

To control the mapping process, two parameters have been chosen: the upper Vupper and lower Vlower limits of the mapped velocity


On that basis, the upper ωupper and the lower ωlower cut-off frequencies of the band-pass filters were determined:


They are necessary for determination of the frequency band from the functions xodd[p,k] and xeven[p,k].

To compare with the previously described algorithm,20 this algorithm has better tuning of the chosen velocity diapason. It is determined not using the central velocity and its accuracy, but directly by independently choosing the upper and the lower velocity values. In addition, the ωupper and ωlower determinations have been also digitally improved. Currently, they are directly determined from the Doppler angle and the source central wavelength.

Functions xodd[p,k], xeven[p,k], and xsum[p,k] are processed in three stages: (1) short-time Fourier transform; (2) Hilbert transform, and (3) taking the logarithm of the signal envelope. After that functions xodd[m,k], xeven[m,k], and xsum[m,k] are obtained and considered. They represent three B-scans and correspond to the three parts of the image, Figs. 3(a) to 3(c) and 4(a) to 4(c), respectively. Note that m=1,2,,M, where M is the number of samples per one A-scan after the above processing.

Fig. 3

Structural and one specific velocity (OSV) images of chosen velocities in a flow with a parabolic profile. Initial image (a); structural image of the tube (b); set aside velocity image (c); image of the parabolic profile with band-pass filters 6 to 9 kHz (d); 9 to 12 kHz (e); 12 to 15 kHz (f); 14 to 17 kHz (g); complexed image for 9 to 12 kHz filter (h). Dimensions of the pictures are 0.7×0.7mm.


Fig. 4

Standard structural image (a); image of odd A-scans (b); image of even A-scans (c); OSV image of positively directed flow (d); OSV image of negatively directed flow (e); OSV image of bidirectional flow (f); the resulting image—(g). Dimensions of the pictures are 0.7mm×0.7mm.


Automatic identification of the directions of the flow xtrends[m,k] is performed as the element-wise subtraction of xeven[m,k] from xodd[m,k] with a shift by one column:

and subsequent logical analysis of the obtained results. Functions characterizing positively xpos[m,k] and negatively xneg[m,k] directed flows are determined as follows:

Then B-scans corresponding to the flows xpos[m,k] and xneg[m,k] are color coded considering the experimentally determined minimal level of the constant signal, L, which is separated from the noise (Figs. 3 and 4):


After that functions xmpos[m,k] and xmneg[m,k] are complexed into a common array (OSV map) representing the flow direction and location xcart[m,k] [Fig. 4(f)]:


Function xsum[m,k] is encoded to be at least 20 times higher than the functions xcart[m,k] and then is complexed with them:


The resulting function xlast[m,k] is normalized to the color scale, which includes shades of green to visualize the image of the phantom structure; red and blue colors are used to visualize the direction of the flow; the white color is used as a reference point (no movement, no reflection). To minimize complexation distortions, a proportionality factor between the xsum[m,k] and xcart[m,k] values equal to 20 was estimated as the best one. This allows simultaneously increasing the OSV accuracy values by 18% and the computation speed by 13%. The resulting image is presented in Figs. 3(h) and 4(g).


Results and Discussion

The described algorithm is realized in a form of a software package in graphical programming language «G» of LabVIEW. A simple model of a blood vessel was used for software testing and debugging. It is a transparent tube having an inner diameter of 0.55 mm. Figure 3 shows a structural image obtained by the system—(a); structural image of the tube—(b); structural image of a selected speed in the flow—(c) and final complexed image—(h).

The OSV image of the parabolic profile with other band-pass filters: 6 to 9 kHz (3.7 to 5.5mm/s)—Fig. 3(d), 9 to 12 kHz (5.5 to 7.3mm/s)—Fig. 3(e), 12 to 15 kHz (7.3 to 9.1mm/s)—Fig. 3(f), 14 to 17 kHz (8.5 to 10.4mm/s)—Fig. 3(g), demonstrate the capabilities of the proposed algorithm.

Figure 4(a) shows the results obtained when scanning was performed at a distance of 0.5 mm from the entry. The original structural image contains information about the steady parts as well as the positive and negative directions of the flow (bright stripes). But it is not possible to visually distinguish one from the others, even by changing the analog and digital filters. In the automated mode, the proposed algorithm allows the separation of initial image into three parts for about 1 to 2 s. The proposed method of OSV Doppler mapping splits the original B-scan into three parts and independently reconstructs the structural image and 2-D color-coded OSV images [Figs. 4(b) to 4(e)] with their subsequent complexation [Figs. 4(f) and 4(g)]. By choosing the upper, Vupper, and the lower, Vlower, limits of the mapped velocity, it is possible to change the predetermined specific velocity and the accuracy of the equivelocity map.



The method of Doppler mapping of a chosen velocity is applied on the basis of low coherence depth discrimination of the reflected signal and detection of the introduced carrier frequency changes. It gives 2-D maps with a quantitative functional (rather than statistical) relation between the Doppler shifts of the carrier and flow velocity without a 2π ambiguity. Additional color coding enables automated separation of the standard OCT structural image in to three parts: steady image and two velocity images, even in the more complex flow used to compare to the previously described image.20 Note that velocity images are predetermined and quantitatively mapped and are proportional to the Doppler shift of the carrier itself. Colors are used as an additional option of the algorithm to separate a stationary object, as well as the positive and negative directions of the flow. It is possible to apply this approach in the study of small animals’ hearts where similar flows with complex geometry appear in nature. Since color (Doppler) ultrasound systems have similar disadvantages as the low coherence optical ones, we suggest applying this algorithm in Doppler ultrasound diagnostics to study bigger biomedical objects.


This work was supported by the Government Contract 12.849.2014/K in Biomedical Diagnostics.



C. Kasai, K. Namekawa and A. Koyano, “Omoto real-time two-dimensional blood flow imaging using an autocorrelation technique,” IEEE Trans. Sonics. Ultrason. 32(3), 458–464 (1985).IESUAU0018-9537http://dx.doi.org/10.1109/T-SU.1985.31615Google Scholar


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).OPEXFF1094-4087http://dx.doi.org/10.1364/OL.25.000114Google Scholar


L. Wang et al., “Frequency domain phase-resolved optical Doppler and Doppler variance tomography,” Opt. Commun. 242(4–6), 345–350 (2004).OPCOB80030-4018http://dx.doi.org/10.1016/j.optcom.2004.08.035Google Scholar


Z. Ding et al., “Real-time phase-resolved optical coherence tomography and optical Doppler tomography,” Opt. Express 10(5), 236–245 (2002).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.10.000236Google Scholar


V. X. D. Yang et al., “Improved phase-resolved optical Doppler tomography using the Kasai velocity estimator and histogram segmentation,” Opt. Commun. 208(4–6), 209–214 (2002).OPCOB80030-4018http://dx.doi.org/10.1016/S0030-4018(02)01501-8Google Scholar


D. Morofke, M. C. Kolios and V. X. D. Yang “Two-dimensional velocity estimation for doppler optical coherence tomography,” Proc. SPIE 6429, 64292F (2007).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.701372Google Scholar


D. Morofke et al., “Wide dynamic range detection of bidirectional flow in Doppler optical coherence tomography using a two-dimensional Kasai estimator,” Opt. Lett. 32(3), 253–5 (2007).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.32.000253Google Scholar


Y. Zhao et al., “Doppler variance imaging for clinical monitoring of in vivo human skin blood flow,” Opt. Lett. 25(18), 1358–1360 (2000).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.25.001358Google Scholar


B. R. White et al., “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical Doppler tomography,” Opt. Express 11(25), 3490–3497 (2003).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.11.003490Google Scholar


R. K. Wang and L. An “Doppler optical micro-angiography for volumetric imaging of vascular perfusion in vivo,” Opt. Express 17(11), 8926–8940 (2009).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.17.008926Google Scholar


L. An, J. Qin and R. K. Wang “Ultrahigh sensitive optical microangiography for in vivo imaging of microcirculations within human skin tissue beds,” Opt. Express 18(8), 8220–8228 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.008220Google Scholar


A. Doronin and I. Meglinski “Imaging of subcutaneous microcirculation vascular network by double correlation optical coherence tomography,” Laser Photon. Rev. 7(5), 797–800 (2013).1863-8880http://dx.doi.org/10.1002/lpor.2013.7.issue-5Google Scholar


A. H. Bachmann et al., “Resonant Doppler flow imaging and optical vivisection of retinal blood vessels,” Opt. Express 15(2), 408–422 (2007).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.15.000408Google Scholar


E. Koch et al., “Resonant Doppler imaging with common path OCT,” Proc. SPIE 7372, 737220 (2009).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.831852Google Scholar


J. K. Barton and S. Stromski “Flow measurement without phase information in optical coherence tomography images,” Opt. Express 13(14), 5234–5239 (2005).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.13.005234Google Scholar


R. Vinayakrishnan “Speckles in laser Doppler perfusion imaging,” PhD Thesis, Univ. of Twente, Enschede, The Netherlands (2007).Google Scholar


J. M. Schmitt, S. H. Xiang and K. M. Yung “Speckle in optical coherence tomography,” J. Biomed. Opt. 4(1), 95–105 (1999).JBOPFO1083-3668http://dx.doi.org/10.1117/1.429925Google Scholar


S. G. Proskurin and I. V. Meglinski “Optical coherence tomography imaging depth enhancement by superficial skin optical clearing,” Laser Phys. Lett. 4(11), 824–826 (2007).1612-2011http://dx.doi.org/10.1002/lapl.200710056Google Scholar


S. G. Proskurin and R. K. Wang “One specific velocity visualization in flows with complex geometry,” Proc. SPIE 5696, 129–135 (2005).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.584617Google Scholar


S. G. Proskurin, A. Yu. Potlov and S. V. Frolov, “Doppler mapping of an alternating-sign flow with a complex structure using optical coherence tomography,” Quantum Electron. 44(1), 54–58 (2014).QUELEZ1063-7818http://dx.doi.org/10.1070/QE2014v044n01ABEH015232Google Scholar


Sergey G. Proskurin is an associate professor of biomedical engineering at Tambov State Technical University, Russia. He graduated and postgraduated from the Nonlinear Optics of Physics Department of M. V. Lomonosov Moscow State University, Russia. Then he received several international awards and grants working at the Mechanical Engineering Laboratory and Hamamatsu Photonics KK, Tsukuba, Japan, Uppsala University (Sweden), Cranfield University (United Kingdom). Key research interests are biophotonics and biomedical optics.

Anton Yu. Potlov is a postgraduate student and assistant lecturer of biomedical engineering at Tambov State Technical University, Russia. He received his BSc degree in biomedical engineering from Tambov State Technical University in 2010 and his MSc degree in biomedical engineering from Tambov State Technical University in 2012. Since 2012, he has been a PhD student in the specialty of medical tools, devices, and systems. He has been an SPIE member since 2014.

Sergej V. Frolov has been a professor and chairman of biomedical engineering at Tambov State Technical University, Russia, since 2006. He graduated and postgraduated from the Tambov Institute of Chemical Engineering with Doctor of Engineering degrees in two specialties: automation of technological processes and production, and application of computer technology, mathematical modeling and mathematical methods in scientific research in 1999. His thesis is titled “Methods, algorithms and interval-system guarantee optimum process control.”

Sergey G. Proskurin, Anton Y. Potlov, Sergey V. Frolov, "One specific velocity color mapping using optical coherence tomography," Journal of Biomedical Optics 20(5), 051034 (13 February 2015). http://dx.doi.org/10.1117/1.JBO.20.5.051034
Submission: Received ; Accepted

Doppler tomography

Optical coherence tomography

Signal processing

Bandpass filters

Doppler effect


Phase shifts

Back to Top