Doppler optical coherence tomography in cardiovascular applications

The study of flow dynamics in complex geometry vessels is highly important in various biomedical applications where the knowledge of the mechanic interactions between the moving fluid and the housing media plays a key role for the determination of the parameters of interest, including the effect of blood flow on the possible rupture of atherosclerotic plaques. Doppler Optical Coherence Tomography (DOCT), as a functional extension of Optical Coherence Tomography (OCT), is an optic, non-contact, noninvasive technique able to achieve detailed analysis of the flow/vessel interactions. It allows simultaneous high resolution imaging (∼10 μm typical) of the morphology and composition of the vessel and determination of the flow velocity distribution along the measured cross-section. We applied DOCT system to image high-resolution one-dimensional and multi-dimensional velocity distribution profiles of Newtonian and non-Newtonian fluids flowing in vessels with complex geometry, including Y-shaped and T-shaped vessels, vessels with aneurism, bifurcated vessels with deployed stent and scaffolds. The phantoms were built to mimic typical shapes of human blood vessels, enabling preliminary analysis of the interaction between flow dynamics and the (complex) geometry of the vessels and also to map the related velocity profiles at several inlet volume flow rates. Feasibility studies for quantitative observation of the turbulence of flows arising within the complex geometry vessels are discussed. In addition, DOCT technique was also applied for monitoring cerebral mouse blood flow in vivo. Two-dimensional DOCT images of complex flow velocity profiles in blood vessel phantoms and in vivo sub-cranial mouse blood flow velocities distributions are presented.     

Self-referenced Doppler optical coherence tomography

Doppler optical coherence tomography (DOCT) allows simultaneous micrometer-scale resolution cross-sectional imaging of tissue structure and blood flow. We demonstrate a fiber-optic polarization-diversity-based differential phase contrast DOCT system as a method to perform self-referenced velocimetry in highly scattering media. Using this strategy, we reduced common-mode interferometer noise to <1 Hz and improved Doppler estimates in a scattering flow phantom by a factor of 5.     

Vessel boundary extraction based on a global and local deformable physical model with variable stiffness

Reliable and efficient vessel cross-sectional boundary extraction is very important for many medical magnetic resonance (MR) image studies. General purpose edge detection algorithms often fail for medical MR images processing due to fuzzy boundaries, inconsistent image contrast, missing edge features, and the complicated background of MR images. In this regard, we present a vessel cross-sectional boundary extraction algorithm based on a global and local deformable model with variable stiffness. With the global model, the algorithm can handle relatively large vessel position shifts and size changes. The local deformation with variable stiffness parameters enable the model to stay right on edge points at the location where edge features are strong and at the same time, fit a smooth contour at the location where edge features are missing. Directional gradient information is used to help the model to pick correct edge segments. The algorithm was used to process MR cine phase-contrast images of the aorta from 20 volunteers (over 500 images) with excellent results.     

Doppler optical coherence tomography with a micro-electro-mechanical membrane mirror for high-speed dynamic focus tracking

An elliptical microelectromechanical system (MEMS) membrane mirror is electrostatically actuated to dynamically adjust the optical beam focus and track the axial scanning of the coherence gate in a Doppler optical coherence tomography (DOCT) system at 8 kHz. The MEMS mirror is designed to maintain a constant numerical aperture of approximately 0.13 and a spot size of approximately 6.7 microm over an imaging depth of 1mm in water, which improves imaging performance in resolving microspheres in gel samples and Doppler shift estimation precision in a flow phantom. The mirror's small size (1.4 mm x 1 mm) will allow integration with endoscopic MEMS-DOCT for in vivo applications.     

Wide dynamic range detection of bidirectional flow in Doppler optical coherence tomography using a two-dimensional Kasai estimator

We demonstrate extended axial flow velocity detection range in a time-domain Doppler optical coherence tomography (DOCT) system using a modified Kasai velocity estimator with computations in both the axial and transverse directions. For a DOCT system with an 8 kHz rapid-scanning optical delay line, bidirectional flow experiments showed a maximum detectable speed of >56 cm/s using the axial Kasai estimator without the occurrence of aliasing, while the transverse Kasai estimator preserved the approximately 7 microm/s minimum detectable velocity to slow flow. By using a combination of transverse Kasai and axial Kasai estimators, the velocity detection dynamic range was over 100 dB. Through a fiber-optic endoscopic catheter, in vivoM-mode transesophageal imaging of the pulsatile blood flow in rat aorta was demonstrated, for what is for the first time to our knowledge, with measured peak systolic blood flow velocity of >1 m/s, while maintaining good sensitivity to detect aortic wall motion at <2 mm/s, using this 2D Kasai technique.     

Experimental study of the multiple scattering effect on the flow velocity profiles measured in Intralipid phantoms by DOCT

Time domain Doppler Optical Coherence Tomography (DOCT) technique was applied to measure flow velocity profiles in highly scattering media. We analyzed the distortions of the measured velocity profiles of the 1% Intralipid solution flow embedded into the scattering medium at different embedding depths. For this purpose a tissue phantom consisting of a plain glass capillary (inner diameter 0.3 mm) embedded into a slab of Intralipid solution mimicking human skin was designed. The measured flow velocity profiles and behavior of distortions caused by multiple scattering are shown.     

光学相干层析多普勒成像功能拓展研究

光学多普勒成像(Optical Doppler tomography,ODT)是一种结合了光学相干层析成像技术(Opticalcoherence tomography,OCT)和多普勒流速仪的非侵入、非接触的成像技术,能够实现对高散介质组织内部的血管分布和血液流速的探测。阐述了基于数字希尔伯特变换的相位分离多普勒光学相干层析成像技术的工作原理,并且通过对玻璃毛细管和生物芯片微通道管中聚苯乙烯溶液流速的实验测量,准确测量管内微粒缓慢移动时的多普勒频移量,获得了玻璃管内和生物芯片微通道管中流速分布曲线,证实了所提方法的可行性。获取的多普勒图像具有较高的空间分辨力和速度分辨力,在未来的临床应用中有潜在的应用价值。     

Non-stationary multifractality in stock returns

We perform an extensive empirical analysis of scaling properties of equity returns, suggesting that financial data show time varying multifractal properties. This is obtained by comparing empirical observations of the weighted generalised Hurst exponent (wGHE) with time series simulated via Multifractal Random Walk (MRW) by Bacry et al. [E. Bacry, J. Delour, J.-F. Muzy, Multifractal random walk, Physical Review E 64 (2) (2001) 026103]. While dynamical wGHE computed on synthetic MRW series is consistent with a scenario where multifractality is constant over time, fluctuations in the dynamical wGHE observed in empirical data are not in agreement with a MRW with constant intermittency parameter. We test these hypotheses of constant multifractality considering different specifications of MRW model with fatter tails: in all cases considered, although the thickness of the tails accounts for most of the anomalous fluctuations of multifractality, it still cannot fully explain the observed fluctuations.     

Quantitative retinal-blood flow measurement with three-dimensional vessel geometry determination using ultrahigh-resolution Doppler optical coherence angiography

Retinal blood flow quantification by retinal vessel segmentation with Doppler optical coherence angiography is presented. Vessel diameter, orientation, and position are determined in an en face vessel image and two representative cross-sectional flow images of the vessel. Absolute blood flow velocity is calculated with the help of the measured Doppler frequency shift and determined vessel angle. The volumetric flow rate is obtained with the position and the region of the vessel lumen. The volumetric blood flow rate of retinal arteries before and after a bifurcation is verified in a healthy human eye.     

Modeling electricity spot prices using mean-reverting multifractal processes

We discuss stochastic modeling of volatility persistence and anti-correlations in electricity spot prices, and for this purpose we present two mean-reverting versions of the multifractal random walk (MRW). In the first model the anti-correlations are modeled in the same way as in an Ornstein–Uhlenbeck process, i.e. via a drift (damping) term, and in the second model the anti-correlations are included by letting the innovations in the MRW model be fractional Gaussian noise with H<1/2$H<1/2$. For both models we present approximate maximum likelihood methods, and we apply these methods to estimate the parameters for the spot prices in the Nordic electricity market. The maximum likelihood estimates show that electricity spot prices are characterized by scaling exponents that are significantly different from the corresponding exponents in stock markets, confirming the exceptional nature of the electricity market. In order to compare the damped MRW model with the fractional MRW model we use ensemble simulations and wavelet-based variograms, and we observe that certain features of the spot prices are better described by the damped MRW model. The characteristic correlation time is estimated to approximately half a year.     

改进的经验模态分解法分离超声多普勒血流与管壁信号

超声多普勒血流信号常包含管壁信号的干扰,准确分离二者对提高血流检测的精度具有重要作用。本文提出两种改进的经验模态分解(EMD)方法,先将含管壁信号的超声多普勒信号分解成多层本征模态函数(IMF),然后根据血流信号与管壁信号的不同特性,对既含管壁信号又含血流信号的IMF分量进行分离处理,最后将各层IMF分量中的管壁成分叠加得到管壁信号的估计,而血流信号可通过原信号减去估计的管壁信号而得到。将本方法用于计算机仿真信号和人体实测的超声多普勒信号,并与高通滤波器法、空间选择性降噪法和原EMD法进行比较,结果表明:本文提出的两种方法能在较大的管壁搏动速度范围内准确地分离血流信号和管壁信号,其平均相对误差比高通滤波器的结果降低了约52%和57%。可见,本文提出的两种方法有望用于血流信号与管壁信号的准确分离。      

Imaging of subcutaneous blood vessels and flow velocity profiles by optical coherence tomography

We have applied a compact low power rapid scanning Doppler Optical Coherence Tomography system to monitor multi-dimensional velocity profiles within the complex vessels and simultaneous real-time non-invasive imaging of skin tissues morphology in vivo, in the wavelength range of 1.3–1.5 nm. Optical clearing of skin tissues has been utilized to achieve depth of OCT images up to 1.7 mm. Current approach enables applying low-power (0.4–0.5 mW) and low-noise broadband near-infrared light sources and obtaining OCT images with down to 12 μm spatial resolution. Two-dimensional time-domain OCT images of complex flow velocity profiles in blood vessel phantom and in vivo subcutaneous human skin tissues are presented. The effect of optical clearing on in vivo images is demonstrated and discussed.     

Doppler ultrasound signals simulation from vessels with various stenosis degrees

As a non-invasive method, the Doppler ultrasound technique is used to detect the vessel stenosis. To search for characteristics of Doppler ultrasound signals sensitive to the stenosis, a computer simulation approach is proposed in this paper to generate Doppler ultrasound signals from vessels with various stenosis degrees. The blood flow velocity distribution in a stenosed vessel is firstly calculated using the transient finite element method (FEM). Then the power spectral density of Doppler signals is estimated using the overall-distribution nonparametric estimation method. Finally Doppler signals are generated using the cosine-superposed method. The proposed approach is proved to be useful for simulating Doppler ultrasound signals from vessels with various stenosis degrees. It is also shown that characteristics of Doppler ultrasound signals may be used to estimate the vessel's stenosis degree.     

Intraluminal fiber-optic Doppler imaging catheter for structural and functional optical coherence tomography

We describe a miniature fiber-optic Doppler imaging catheter for integrated functional and structural optical coherence tomography (OCT) imaging. The Doppler catheter can map blood flow within a vessel as well as image vessel wall structures. A prototype Doppler catheter has been developed and demonstrated for measuring the intraluminal velocity profile in a vessel phantom (conduit). A simple mathematical model is demonstrated to estimate the total flow rate. This estimation technique also enables the spatial range of flow measurements to be extended by approximately two times the normal OCT image-penetration depth. The Doppler OCT catheter could be a powerful device for cardiovascular imaging.     

Experimental study of the effects of pulsed Doppler sample volume size and position on the Doppler spectrum

With a pulsed Doppler system, the recorded Doppler spectrum is expected to vary depending upon the sample volume size relative to the diameter of the vessel, the position of the sample volume in the vessel and the velocity profile. In the in vitro experiments described in this paper, the velocity profile was kept constant by using steady parabolic flow in a flow model. As the Doppler sample volume size and position were changed, the maximum variations of quantitative measurements from the Doppler spectrum were determined. The maximum, mean and mode frequencies and spectral broadening index (SBI) were affected by the position of the sample volume but to a lesser degree by its length (1.5-5.0 mm) relative to the 9.5 mm beam path length across the tube. When the centre of the Doppler sample volume was moved within the central 25% of the tube, the maximum variations were as follows: maximum frequency 3-5%, mean frequency 8-9%, mode frequency 8-9% and SBI 16-18%, where the range indicates the effect of increasing the sample volume size. Based on these results obtained under steady flow conditions in vitro, it is concluded that quantification of pulsed Doppler spectra may be feasible if the sample volume is positioned within the central 25% of the vessel.     

Simulation of blood flow using extended Boltzmann kinetic approach

Lattice Boltzmann (LB) simulations are conducted to obtain the detailed hydrodynamics in a variety of blood vessel setups, including a prototype stented channel and four human coronary artery geometries based on the images obtained from real patients. For a model of stented flow involving an S-shape stent, a pulsatile flow rate is applied as the inlet boundary condition, and the time- and space-dependent flow field is computed. The LB simulation is found to reproduce the analytical solutions for the velocity profiles and wall shear stress distributions for the pulsatile channel flow. For the coronary arteries, the distributions of wall shear stress, which is important for clinical diagnostic purposes, are in good agreement with the conventional CFD predictions.     

Multiple contrast fast spin-echo approach to black-blood intracranial MRA: use of complementary and supplementary information.

Black-blood magnetic resonance angiography (black-blood MRA) could be considered an alternative to time-of-flight (TOF) MRA. In the cases of irregular flow conditions, it could be more advantageous than time-of-flight (TOF) MRA in providing vessel definition and delineation. Proton-density weighted (PDW) multi-slab three-dimensional fast spin-echo (3DFSE) sequences have been used to generate black-blood MRA. Unfortunately, multi-planar reformatted 3DFSE images often exhibit slab boundary artifacts (intensity variation in the slice direction) which create dark bands interfering with the identification of dark blood vessels. Furthermore, PDW measurements fail to darken slow flowing or re-circulating blood in some circumstances. In this work, a dual-contrast multi-slab 3DFSE acquisition is used to approach black-blood MRA. This sequence simultaneously provides proton-density weighted (PDW) and T(2)-weighted (T2W) images which can be further integrated together to produce black-blood angiograms gained by utilizing complementary contrast and supplementary vascular information. Additionally, a technique of suppressing slab boundary artifact has been incorporated into this sequence. This approach provides: i) good SNR measurement of anatomy for the PDW image and optimal black-blood angiograms from the T2W image; ii) scan time efficiency (dual-contrast image sets plus black-blood angiograms within one acquisition); and iii) suppressed slab boundary artifacts as well as minimized mis-registration error.     

基于超快超声多普勒的三维脑损伤成像方法研究

三维超声微血管成像可直观呈现血流信息,对于脑血管疾病诊断和治疗具有重要意义。本文旨在将超快超声成像技术、超快超声功率多普勒技术和机械扫描相结合,实现脑血管三维成像和脑缺血区域评价。通过工程实现,完成了可同步控制微型线性位移平台移动和超声阵列超快发射、高速采集与压缩存储的三维扫描数据采集序列与系统。利用GPU并行运算,高效实现了超声图像波束合成方法,对原始射频超声数据完成重建。进而,基于SVD杂波滤除技术,从重建三维超声数据中提取了脑部的动态小血管信号,并获得了各切面的功率多普勒成像和冠状面彩色多普勒超声小血管成像。最后,采用体素方法对三维脑血管进行重建。大鼠在体实验结果表明,该成像系统可用于三维脑血管网络在体成像,以及脑血管损伤区域定位与量化评价。本工作对脑病检测技术发展与诊断方法研究具有一定的借鉴意义。此外,相关检测系统和成像算法具有一定普适性,对其他富含微血流血管的组织检测也有一定的参考价值。     

Intracranial vascular feature changes in time of flight MR angiography in patients undergoing carotid revascularization surgery

PurposeTo characterize the intracranial vascular features extracted from time of flight (TOF) images and their changes from baseline to follow-up in patients undergoing carotid revascularization, using arterial spin labeling (ASL) cerebral blood flow (CBF) measurement as a reference.MethodsIn this retrospective study, brain TOF and ASL images of 99 subjects, acquired before, within 48 h, and/or 6 months after, carotid revascularization surgery were analyzed. TOF images were analyzed using a custom software (iCafe) to quantify intracranial vascular features, including total vessel length, total vessel volume, and number of branches. Mean whole-brain CBF was calculated from ASL images. ASL scans showing low ASL signal in the entire flow territory of an internal carotid artery (ICA), which may be caused by labeling failure, were excluded. Changes and correlations between time points were analyzed separately for TOF intracranial vascular features and ASL CBF.ResultsSimilar to ASL CBF, TOF vascular features (i.e. total vessel length, total vessel volume and number of branches) increased dramatically from baseline to post-surgery, then returned to a level slightly higher than the baseline in long-term follow-up (All P < 0.05). Correlation between time points was observed for all three TOF vascular features but not for ASL CBF.ConclusionIntracranial vascular features, including total vessel length, total vessel volume and number of branches, extracted from TOF images are useful in detecting brain blood flow changes induced by carotid revascularization surgery.