Frequent Links
Doppler echocardiography
Doppler echocardiography  

Intervention  
File:Doppler mitral valve.gif  
MeSH  D015150 
OPS301 code:  3052 
Doppler echocardiography is a procedure that uses ultrasound technology to examine the heart or blood vessels.^{[1]} An echocardiogram uses high frequency sound waves to create an image of the heart while the use of Doppler technology allows determination of the speed and direction of blood flow by utilizing the Doppler effect.
An echocardiogram can, within certain limits, produce accurate assessment of the direction of blood flow and the velocity of blood and cardiac tissue at any arbitrary point using the Doppler effect. One of the limitations is that the ultrasound beam should be as parallel to the blood flow as possible. Velocity measurements allow assessment of cardiac valve areas and function, any abnormal communications between the left and right side of the heart, any leaking of blood through the valves (valvular regurgitation), calculation of the cardiac output and calculation of E/A ratio^{[2]} (a measure of diastolic dysfunction). Contrastenhanced ultrasound using gasfilled microbubble contrast media can be used to improve velocity or other flowrelated medical measurements.
An advantage of Doppler echocardiography is that it can be used to measure blood flow within the heart without invasive procedures such as cardiac catheterization. For example, left ventricular filling pressure can be estimated by Doppler ultrasound, although only under certain conditions.^{[3]}
Although "Doppler" has become synonymous with "velocity measurement" in medical imaging, in many cases it is not the frequency shift (Doppler shift) of the received signal that is measured, but the phase shift (when the received signal arrives).
This procedure is frequently used to examine children's hearts for heart disease because there is no age or size requirement.
2D Doppler Imaging
Unlike 1D Doppler imaging, which can only provide onedimensional velocity and has dependency on the beam to flow angle,^{[4]} 2D velocity estimation using Doppler ultrasound is able to generate velocity vectors with axial and lateral velocity components. 2D velocity is useful even if complex flow conditions such as stenosis and bifurcation exist. There are two major methods of 2D velocity estimation using ultrasound: Speckle tracking and crossed beam Vector Doppler, which are based on measuring the time shifts and phase shifts respectively.^{[5]}
Vector Doppler
Vector Doppler is a natural extension of the traditional 1D Doppler imaging based on phase shift. The phase shift is found by taking the autocorrelation between echoes from two consecutive firings.^{[6]} The main idea of Vector Doppler is to divide the transducer into three apertures: one at the center as the transmit aperture and two on each side as the receive apertures. The phase shifts measured from left and right apertures are combined to give the axial and lateral velocity components. The positions and the relative angles between apertures need to be tuned according to the depth of the vessel and the lateral position of the region of interest.^{[5]}
Speckle tracking
Speckle tracking, which is a wellestablished method in video compression and other applications, can be used to estimate blood flow in ultrasound systems. The basic idea of speckle tracking is to find the best match of a certain speckle from one frame within a search region in subsequent frames.^{[5]} The decorrelation between frames is one of the major factors degrading its performance. The decorrelation is mainly caused by the different velocity of pixels within a speckle, as they do not move as a block. This is less severe when measuring the flow at the center, as the changing rate of the velocity is the lowest. The flow at the center usually has the largest velocity magnitude, called peak velocity. It is the most needed information in some cases, such as diagnosing stenosis.^{[7]} There are mainly three methods of finding the best match: SAD (Sum of absolute difference), SSD (Sum of squared difference) and Cross correlation. Assume <math>X_0 (i,j)</math> is a pixel in the kernel and <math>X_1 (i+\alpha,j+\beta)</math> is the mapped pixel shifted by <math>(\alpha,\beta)</math> in the search region.^{[8]}
SAD is calculated as: <math>D(\alpha,\beta)=\sum_{i=1} \sum_{j=1} X_0 (i,j)X_1 (i+\alpha,j+\beta)</math>
SSD is calculated as: <math>D(\alpha,\beta)=\sum_{i=1} \sum_{j=1} (X_0 (i,j)X_1 (i+\alpha,j+\beta))^{2}</math>
Normalized cross correlation coefficient is calculated as: <math>\rho(\alpha,\beta)=\frac{\sum_{i=1} \sum_{j=1}(X_0 (i,j)\bar{X_0})(X_1 (i+\alpha,j+\beta)\bar{X_1})}{\sqrt{(\sum_{i=1} \sum_{j=1} (X_0 (i,j)\bar{X_0})^{2})(\sum_{i=1} \sum_{j=1} (X_1 (i+\alpha,j+\beta)\bar{X_1})^{2})}}</math>
where <math>\bar{X_0}</math> and <math>\bar{X_1}</math> are the average values of <math>X_0 (i,j)</math> and <math>X_1 (i,j)</math> respectively. The <math>(\alpha,\beta)</math> pair that gives the lowest D for SAD and SSD, or the largest ρ for the cross correlation, is selected as the estimation of the movement. The velocity is then calculated as the movement divided by the time difference between the frames. Usually the median or average of multiple estimations is taken to give more accurate result.^{[8]}
Sub pixel accuracy
In ultrasound systems, lateral resolution is usually much lower than the axial resolution. The poor lateral resolution in the Bmode image also results in poor lateral resolution in flow estimation. Therefore, sub pixel resolution is needed to improve the accuracy of the estimation in the lateral dimension. In the meantime, we could reduce the sampling frequency along the axial dimension to save computations and memories if the sub pixel movement is estimated accurately enough. There are generally two kinds of methods to obtain the sub pixel accuracy: interpolation methods, such as parabolic fit, and phase based methods in which the peak lag is found when the phase of the analytic cross correlation function crosses zero.^{[9]}
Interpolation method(parabolic fit)
As shown in the right figure, parabolic fit can help find the real peak of the cross correlation function. The equation for parabolic fit in 1D is:^{[4]} <math>k_{int}=k_s\frac{(R_{12} (k_s+1)R_{12} (k_s1))}{2(R_{12} (k_s+1)2R_{12} (k_s )+R_{12} (k_s1))} </math>
where <math>R_{12}</math> is the cross correlation function and <math>k_s</math> is the originally found peak. <math>k_{int}</math> is then used to find the displacement of scatterers after interpolation. For the 2D scenario, this is done in both the axial and lateral dimensions. Some other techniques can be used to improve the accuracy and robustness of the interpolation method, including parabolic fit with bias compensation and matched filter interpolation.^{[10]}
Phase based method
The main idea of this method is to generate synthetic lateral phase and use it to find the phase that crosses zero at the peak lag.^{[9]}
The right figure illustrates the procedure of creating the synthetic lateral phase, as a first step. Basically, the lateral spectrum is split in two to generate two spectra with nonzero center frequencies. The cross correlation is done for both the up signal and down signal, creating <math>R_{up}</math> and <math>R_{down}</math> respectively.^{[9]} The lateral correlation function and axial correlation function are then calculated as follows: <math>R_{lateral}=R_{up}*R_{down}^{*}; R_{axial}=R_{up}*R_{down}</math>
where <math>R_{down}^{*}</math> is the complex conjugate of <math>R_{down}</math>.
It is obvious that they have the same magnitude, and the integer peak is found using traditional cross correlation methods. After the integer peak is located, a 3 by 3 region surrounding the peak is then extracted with its phase information. For both the lateral and axial dimensions, the zero crossings of a onedimensional correlation function at the other dimension’s lags are found, and a linear least squares fitted line is created accordingly. The intersection of the two lines gives the estimate of the 2D displacement.^{[9]}
Comparison between Vector Doppler and Speckle tracking
Both methods could be used for 2D Velocity Vector Imaging, but Speckle Tracking would be easier to extend to 3D. Also, in Vector Doppler, the depth and resolution of the region of interest are limited by the aperture size and the maximum angle between the transmit and receive apertures, while Speckle Tracking has the flexibility of alternating the size of the kernel and search region to adapt to different resolution requirement. However, vector Doppler is less computationally complex than speckle tracking.
Volumetric flow estimation
Velocity estimation from conventional Doppler requires knowledge of the beamtoflow angle (Inclination angle) to produce reasonable results for regular flows and does a poor job of estimating complex flow patterns, such as those due to stenosis and/or bifurcation. Volumetric flow estimation requires integrating velocity across the vessel crosssection, with assumptions about the vessel geometry, further complicating flow estimates. 2D Doppler data can be used to calculate the volumetric flow in certain integration plane.^{[11]} The integration plane is chosen to be perpendicular to the beam, and Doppler power (generated from power Doppler mode of Doppler ultrasound) can be used to differentiate between the components that are inside and outside the vessel. This method does not require prior knowledge of the Doppler angle, flow profile and vessel geometry.^{[11]}
Promise of 3D
Until recently, ultrasound images have been 2D views and have relied on highly trained specialists to properly orient the probe and select the position within the body to image with only few and complex visual cues. The complete measurement of 3D velocity vectors makes many post processing techniques possible. Not only the volumetric flow across any plane is measurable, but also other physical information such as stress and pressure can be calculated based on the 3D velocity field. However, it is quite challenging to measure the complex blood flow to give velocity vectors, due to the fast acquisition rate and the massive computations needed for it. Plane wave technique is thus promising as it can generate very high frame rate.^{[12]}
See also
 Medical ultrasonography section: Doppler sonography
 Echocardiography
 American Society of Echocardiography
 Christian Doppler
References
 ^ Doppler Echocardiography at the US National Library of Medicine Medical Subject Headings (MeSH)
 ^ [1] Abdul Latif Mohamed, Jun Yong, Jamil Masiyati, Lee Lim, Sze Chec Tee. The Prevalence Of Diastolic Dysfunction In Patients With Hypertension Referred For Echocardiographic Assessment of Left Ventricular Function. Malaysian Journal of Medical Sciences, Vol. 11, No. 1, January 2004, pp. 6674
 ^ Ommen, S. R.; Nishimura, R. A.; Appleton, C. P.; Miller, F. A.; Oh, J. K.; Redfield, M. M.; Tajik, A. J. (10 October 2000). "Clinical Utility of Doppler Echocardiography and Tissue Doppler Imaging in the Estimation of Left Ventricular Filling Pressures : A Comparative Simultaneous DopplerCatheterization Study". Circulation. Circulation. pp. 1788–1794. doi:10.1161/01.CIR.102.15.1788. Retrieved 12 July 2012.
 ^ ^{a} ^{b} J. A. Jensen, Estimation of Blood Velocities Using Ultrasound, A Signal Processing Approach, New York: Cambridge University Press, 1996.
 ^ ^{a} ^{b} ^{c} P. S. a. L. L. Abigail Swillens, "TwoDimensional Blood Velocity Estimation With Ultrasound: Speckle Tracking Versus CrossedBeam Vector Doppler Based on Flow Simulations in a Carotid Bifurcation Model," IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, pp. 327338, 2010.
 ^ R. S. C. Cobbold, Foundations of Biomedical Ultrasound, Oxford University Press, 2007.
 ^ G. Reutern, M. Goertler, N. Bornstein, M. Sette, D. Evans, A. Hetzel, M. Kaps, F. Perren, A. Razumovky, T. Shiogai, E. Titianova, P. Traubner, N. Venketasubramanian, L. Wong and M. Yasaka, "Grading Carotid Stenosis Using Ultrasonic Methods," Stroke, Journal of the American Heart Association, vol. 43, pp. 916921, 2012.
 ^ ^{a} ^{b} J. Luo and E. E. Konofagou, "A Fast Motion and Strain Estimation," in Ultrasound Symposium, 2010.
 ^ ^{a} ^{b} ^{c} ^{d} X. Chen, M. J. Zohdy, S. Y. Emelianov and M. O'Donnell, "Lateral Speckle Tracking Using Synthetic Lateral Phase," IEEE Transactions on Ultrasonics, Ferroelectrcs and Frequency Control, vol. 51, no. 5, pp. 540550, 2004.
 ^ X. Lai and H. Torp, "Interpolation Methods for TimeDelay Estimation Using CrossCorrelation Method for Blood Velocity Measurement," IEEE Transactions on Ultrasonics, Ferroelectrcs and Frequency Control, vol. 46, no. 2, pp. 277290, 1999.
 ^ ^{a} ^{b} M. Richards, O. Kripfgans, J. Rubin, A. Hall and J. Fowlkes, "Mean Volume Flow Estimation in Pulsatile Flow Conditions," Ultrasound in Med. & Biol., vol. 35, pp. 18801891, 2009.
 ^ J. Udesen, F. Gran, K. Hansen, J. Jensen, C. Thomsen and M. Nielsen, "High Frame Rate Blood Vector Velocity Imaging Using Plane Waves: Simulations and Preliminary Experiments," IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 55, no. 8, pp. 17291743, 2008.
External links
 Echocardiography Textbook by Bonita Anderson
 Echocardiography (Ultrasound of the heart)
 Doppler Examination  Introduction
 The Doppler Principle and the Study of Cardiac Flows
