This article has Open Peer Review reports available.
Can echocardiographic particle image velocimetry correctly detect motion patterns as they occur in blood inside heart chambers? A validation study using moving phantoms
- Christian Prinz†1, 4,
- Reka Faludi†1, 5,
- Andrew Walker2,
- Mihaela Amzulescu1,
- Hang Gao1,
- Tokuhisa Uejima3,
- Alan G Fraser3 and
- Jens-Uwe Voigt1Email author
© Prinz et al.; licensee BioMed Central Ltd. 2012
Received: 3 April 2012
Accepted: 6 June 2012
Published: 6 June 2012
To validate Echo Particle Image Velocimetry (PIV)
High fidelity string and rotating phantoms moving with different speed patterns were imaged with different high-end ultrasound systems at varying insonation angles and frame rates. Images were analyzed for velocity and direction and for complex motion patterns of blood flow with dedicated software. Post-processing was done with MATLAB-based tools (Dflow, JUV, University Leuven).
Velocity estimation was accurate up to a velocity of 42 cm/s (r = 0.99, p < 0.001, mean difference 0.4 ± 2 cm/s). Maximally detectable velocity, however, was strongly dependent on frame rate and insonation angle and reached 42 cm/s under optimal conditions. At higher velocities estimates became random. Direction estimates did depend less on velocity and were accurate in 80-90%. In-plane motion patterns were correctly identified with three ultrasound systems.
Echo-PIV appears feasible. Velocity estimates are accurate, but the maximal detectable velocity depends strongly on acquisition parameters. Direction estimation works sufficiently, even at higher velocities. Echo-PIV appears to be a promising technical approach to investigate flow patterns by echocardiography.
Spectral- and color Doppler are powerful echocardiographic methods for imaging and quantifying blood flow velocities. However, they measure only the velocity component along the direction of the ultrasound beam and therefore cannot provide information on the direction or pattern of blood flow.
Recently introduced echocardiographic tracking algorithms for the assessment of myocardial motion stimulated the development of software which is capable of tracking image features in the blood pool in order to estimate blood motion in any direction within the image plane. This concept of Particle Image Velocimetry (PIV) , applied to contrast enhanced echocardiographic images [2, 3], may allow not only to display and quantify blood flow velocity and direction, but may also provide new insights into typical cardiac flow patterns, such as vortices. The obtained information could then be used even to estimate energy dissipation in the flow field [4, 5].
Previous studies have tested the software in comparison to other imaging techniques  or in virtual models  using flow patterns mimicking cardiac flow. So far, no detailed analysis of the influence of imaging geometry in combination with different acquisition settings has been performed.
In this study we used moving phantoms to investigate to what extent a newly developed echocardiographic PIV software (Omega Flow, Siemens, Mountain View, CA, USA*a) is able to correctly describe flow direction and velocity and to recognize flow patterns within the image.
The rotating phantom was used to test both, the flow recognition at different angles as well as the vortex recognition algorithm. A flat cylinder of agar of 3 cm thickness and 10 cm in diameter was cast on a disc of acrylic glass. Holes were punched in the agar with a fine needle in a way that an evenly distributed pattern of distinct speckles was obtained when the phantom was imaged. The disc with the agar phantom was directly fitted to the shaft of the DC-motor of the same testing device  (Figure 1B). The tangential velocity of a certain point in the agar phantom could be calculated by multiplying its distance from the motor axis with 2π and the rotational speed of the motor.
Setup and Imaging Protocol
Phantoms were submerged in a rubber damped tank which was filled with degassed tap water. The ultrasound transducers were placed in a tripod and positioned just below the water surface. A regular attenuator, as used to avoid near field artefact in soft tissue imaging, was placed under the probes (Figure 1A and 1B).
The string phantom could be turned to obtain different insonation angles. Care was taken to center the ultrasound transducer over the tilt axis of the string phantom which was then imaged at a depth of approximately 7.5 cm. Similarly, the rotating phantom was positioned in the center of the image with its rotational axis at a depth of approximately 9.5 cm.
Images were collected with high-end ultrasound systems (Siemens Acuson SC 2000, Siemens Medical Solutions, Mountain View, USA, probe: 4V1c-S, 1.0-4.0 MHz; GE Vivid E9, Vingmed Ultrasound, Horten Norway, probe: M3S, 1.5-4.0 MHz; Philips iE33, Philips Ultrasound, Bothell, WA, USA, probe: S5-1, 1.0-5.0 MHz). Image loops of at least three seconds were digitally stored at three different gain settings (−10, 0, +10 dB) to allow the selection of the optimal image for subsequent off-line analysis. The optimal image was selected, if already in part a pure visual tracking of specific speckle-patterns was possible.
In order to investigate the influence of frame rate, three different acquisition settings resulting in comparable frame rate ranges (38–44 fps, 56–62 fps and 74–81 fps) were used in all ultrasound machines.
For determining the accuracy of velocity and angle estimates, the rotating phantom was driven at different continuous speeds resulting in circumferential velocities between 5 and 400 cm/s.
In order to investigate the limits of tracking, both the string and the rotating phantom were driven using a ramp pattern with increasing string speed up to a varying plateau between 10 and 100 cm/s. The ramp was repeated 60 times per minute.
Particle Image Velocimetry (PIV)
In a second step, the circumscribed area of interest is analyzed by the feature tracking algorithm of the program which provides a matrix of instantaneous local flow velocity vector estimates based on the tracking of the moving scatterers (Figure 3B and E). Generally, the highest possible matrix density setting was used (approx. one vector estimate per 8*8 image pixels). For evaluating the vortex recognition in images from different machines, medium and low matrix density settings were tested as well.
The velocity matrix is then used to calculate the local in-plane vorticity. The software provides an estimate of the position and size of the dominant clockwise and counterclockwise vortices within the circumscribed region. For the latter, the threshold for including detected vortical flow in the main vortex area could be set as 25%, 50%, 75% or 100% of the maximum heartbeat-averaged steady-streaming field (Figure 3C). We defined the biggest continuous area of vortical flow as steady-streaming field.
For string phantom data, the sample volume of the dFlow program was placed in the midline of the image sector directly on the string. In the rotating phantom clips, the sample volume was positioned in a way that different insonation angles could be investigated. The distance of the sample volume to the phantom center of rotation was measured in the image in order to scale the known velocity profile of the phantom to the true velocity in sample position (Figure 2 A and B).
Proof of principle in human hearts
For a proof of principle, we collected contrast enhanced echo data from routine patients of the clinical echo lab in Leuven who had good apical echo quality and normal cardiac structure and function on the routine echocardiogram. The study had been approved by the ethical committee of our institution and all patients had given their informed consent prior to the examination. 20 echocardiographic image loops of good quality but varying frame rate (25 – 179 fps) showing the three- or four chamber view of human hearts after a bolus injection of 0.2 ml SonoVue were analyzed with the PIV software as described above. Flow structures were displayed and cut-off velocities were determined by comparing the tracking derived flow velocity component towards the probe at the tip of the mitral or tricuspid valve leaflets to pulsed wave Doppler traces from the same position during ventricular filling (Figure 4).
Determination of the cut-off velocity
Determination of the accuracy of the velocity estimate
In the string phantom, velocity estimates from the center of the image were compared to the known velocity using the method of Bland and Altman . In the rotating phantom, error histograms of the entire flow field were calculated. In ramp profiles, data below and above the cut-off velocity were considered separately.
Determination of the accuracy of the direction estimate
In the string phantom, a direction estimate from the center of the image was analyzed. A deviation from the true angle of less than ± 10° was arbitrarily accepted as accurate and the relative frequency of accurate measurements was determined in relation to the true velocity of the phantom. In the rotating phantom, the deviation of the angle estimate from the truth was analyzed using error histograms for the entire flow field. In ramp profiles, data below and above the cut-off velocity were considered separately (see example in Figure 5 C and D).
In general, continuous variables are presented in the following as mean ± standard deviation (SD). Linear regression was used to investigate the relation between two parametric variables. Continuous measurements were compared by Spearman correlation and a Bland-Altman-Analysis. A two-tailed p value of < 0.05 was considered significant. Statistical analyses were performed with MATLAB (dito) and the PASW™ software (SPSS Inc., Chicago, Illinois).
The echocardiographic PIV software was able to follow the moving speckle pattern produced by the phantoms. Only acquisitions from the string phantom at 90 degree insonation angle were not analyzable for technical reasons, since the speckles produced by the surgical silk merged to a continuous white line. Estimates of motion velocity and direction were related to the true motion of the phantoms as reported below. The rotating phantom was properly recognized as a vortical structure.
Accuracy of PIV velocity estimates was found to depend strongly on acquisition settings.
Using the string phantom, we could further show a significant correlation between cut-off-velocity and insonation angle with (r = −0.88, p < 0.001) (Figure 6B). The highest cut-off velocity was reached when the motion direction was close to the beam direction and lowest, when close to perpendicular to it.
Direction estimates proved to be less sensitive to changing acquisition settings. We found no significant relation between the frequency of accurate direction estimates and frame rate. However, we noticed a relation to the insonation angle. The error of the direction estimate was lowest when the motion was close to the beam direction (Figure 6C).
Estimates became already accurate at very low velocities (above 4–8 cm/s) and remained stable up to the cut-off-velocity. Above the cut-off-velocity, estimates were less, but still sufficiently accurate (see example in Figure 5C and 5D). Under optimal conditions (low insonation angle in string phantoms), direction estimates deviated less than ±10° from the true direction in 86.7 ± 1.5% of all frames.
The PIV software was able to identify the rotating agar phantom as vortex (see example in Figure 3C). The vortex center was accurately identified with a minor offset of < 5.5 mm. Interestingly, the size of the indicated vortex region varied with the rotational speed of the phantom. The identified vortex had at maximum an approximate diameter of 2/3 of the region of the phantom in which the tangential velocities were still below the cut-off velocity for tracking. In general, the accuracy of the detection of the vortex center (p = 0.5) and the estimation of vortex size (p = 0.4) did not depend on the density settings of the flow vector field.
The maximally detectable velocity (cut-off velocity) in the different high-end ultrasound machines was comparable (44 cm/s, 42 cm/s, and 42 cm/s, n.s.). We found further no significant differences for direction estimation in the different machines (80–90% of all frames within the ±10° range of the true direction for all machines).
The error for the vortex center detection was comparable (4.7 ± 3.5 mm, 5.5 ± 3 mm and 4.0 ± 3.1 mm). Acquisition and post-processing settings had no relevant impact on the accuracy of the vortex pattern recognition in image loops from the different machines.
In concordance with the phantom experiments, the cut-off velocity proved to be linearly dependent on frame rate (r = 0.94, p < 0.01) (Figure 4C).
Flow can be described as laminar, vortical or turbulent. Typically, laminar flow is observed in straight blood vessels, vortex formation occurs in heart chambers, and turbulence may arise from stenotic valves. Vortex and recirculation may also occur in vessel bifurcations and turbulences in vessel stenoses [4, 12–17]. A vortex is a structure able to store kinetic energy while turning. Turbulent flow, characterized by chaotic appearance of vortical structures on different scales, however, leads to rapid dissipation of kinetic energy. Flow patterns may therefore be relevant for an energy-efficient cardiac function and their non-invasive detection by echocardiography may be of clinical interest.
In this study, we used an in-vitro model of string and rotational phantoms to mimic laminar and vortical flow patterns. Our results indicate that PIV-echocardiography is able to estimate flow velocities and direction of moving speckle patterns in these in-vitro models. Accuracy of velocity estimation was strongly dependent on acquisition settings, but turned out to be sufficiently accurate for lower velocities below a certain cut-off velocity. The estimation of flow direction was less sensitive to changing acquisition settings and sufficiently accurate. The mimicked vortical flow pattern was recognized, even in varying settings with different ultrasound machines.
Both the existence of a cut-off velocity and its frame rate dependence, can be explained with the underlying tracking algorithm: a speckle pattern can only be recognised and found back again in the following frame if its displacement is less than the maximum search distance of the algorithm. If the flow velocity is too high or the frame rate too low, speckle patterns move further than this and the tracking fails. In general terms, the cutoff velocity is about the length of the search window (in pixels) * pixelsize * framerate. In a normal echocardiographic study of a heart chamber (pixelsize = 0.25 mm and frame rate 50frames/sec) using a 32x32 window, this means a cutoff velocity of 0.4 m/s. This theoretical assumption is in good agreement with our measurements with a maximal measured cut-off velocity of 0.44 m/s.
This low cut-off velocity is a current limitation of the technique. It will be regularly exceeded already under physiologic conditions in the inflow region of the ventricle. For the remaining areas of the ventricle or the atria , however, it will be fully sufficient. Besides that, pattern recognition will be less affected, since our data demonstrate that direction estimates are still reasonable above the cut-off velocity.
The angle dependency of the cut-off velocities as seen in the practically one-dimensional string phantom can be best explained by the changing image resolution in motion direction when the phantom is tilted: a good resolution in-line with the beam results in higher cut-off velocities while the lower image resolution across beams delivers worse results.
Surprisingly, flow direction estimates remained reasonably accurate even above the cut-off velocity (Figure 5C and 5D). We hypothesize that this phenomenon could be explained if the correlation of speckle patterns between frames is higher in the direction of flow than it is perpendicular to flow.
The displayed vortex size was dependent on the rotational speed of the phantom, but always smaller than the phantom size. While the former is explained by the loss of tracking in the outer regions of the phantom which exceed the cut-off velocity, the latter effect is due to the particularities of the rigid phantom. In contrast to a fluid vortex in the heart which has centrifugally decreasing angular velocities and therefore its highest vorticity in the center, the rigid phantom shows constant angular velocities in all parts which results in highest vorticity values at the outer edge. Under such circumstances, the implemented algorithm - which looks for a compact vorticity region above a certain threshold percentage of the heart beat averaged maximum vorticity – results in a radius of the displayed vortex (rv) proportional to the phantom radius (R): (personal communication with G. Pedrizzetti, University of Trieste, Italy).
We observed no relevant difference between the tracking results from image loops acquired with the different high-end ultrasound machines. This may not fully reflect clinical reality since our phantoms, providing close-to-ideal imaging conditions, did not challenge the imaging capabilities of the different systems. Our results show clearly, however, that care must be taken to preserve original image resolution and frame rate during DICOM-conversion.
Our findings indicate that requirements for image acquisition as revealed by our phantom experiments, apply fully to the clinical settings. Therefore, care must be taken to have ideal imaging conditions, in particular, good spatial resolution and frame rates as high as possible. Since such conditions can not always be found in clinical patients, a limited feasibility of the approach in a routine setting must be expected.
In this study we used phantoms for the validation of PIV echocardiography which provided a defined two-dimensional pattern of speckle motion. In clinical practice, out-of plane motion of contrast speckles may lead to inaccuracies in the assessment of complex vortex structures in human heart chambers [19, 20]. Interpretation of 2D tracking results must always consider, that at best the in-plane component of a 3D flow structure can be assessed. Three-dimensional echocardiography may help to overcome this limitation, but volume rates and spatial resolution of current machines are not sufficient for flow tracking in humans.
Our phantoms provided close-to ideal image quality without artefacts and homogeneous distribution of speckles. Less favorable scanning conditions in cardiac patients and inhomogeneities in contrast distribution may lead to tracking problems and subsequent misinterpretations of the flow field.
In our study, speckle patterns were generated by solid bodies, while in the clinical application they are generated by a non-Newtonian fluid. Together with the spatial character of fluid vortices, this may limit applicability of the model to the in vivo situation.
Given the underlying methodology, it must be assumed that tracking results do also depend on other imaging parameters, such as speckle distribution and density . Therefore the sensitivity of the algorithm to the number of speckles per unit area was not evaluated. The current setup was not suited to investigate this relation. Furthermore, it is unclear how the software determines the periphery of vortex structures, since velocities are gradually declining when increasing the distance from the center. Further studies are needed to address those questions.
Particularly in diastole, kinetic energy of the blood entering the ventricle should be stored and the blood flow should be re-directed for ejection with a minimum of energy loss. A vortex structure appears ideally suited for this and may help to maximize ventricular efficiency. Therefore, echocardiographic PIV can provide new insights into diastolic function and dysfunction of the heart chambers. Further, the method might help to better understand the effect of regional myocardial dysfunction, e.g. after myocardial infarction, or help to optimize surgical interventions, such as valve replacements.
The new method of echocardiographic particle image velocimetry (PIV) appears feasible on high quality, high frame rate DICOM images from different ultrasound machines. Accuracy of tracking results depends strongly on acquisition settings, image quality and absolute flow velocities. While velocity estimates are only accurate below a certain cut-off velocity, motion direction estimates appear more robust, even at higher velocities. This may allow a correct interpretation of flow patterns (relative vorticity) even if absolute flow velocity estimates become inaccurate. The demand on high quality, high frame rate image data appears as a current limitation to the use of the method in the clinical setting.
aAt the time point of study. Currently, the software is marketed under the name “HyperFlow”, AMID, Trieste, Italy.
This study was supported by the Research Foundation Flanders (FWO), research grant G068508N.
We thank Dr. M. Ressner, Dept of Biomedical Engineering, Linköping University, Linköping, Sweden, for his help with the preparation of the agar phantoms.
CP receives a personal grant from the German Heart Foundation (“Deutsche Herzstiftung”) (grant: S/11/10; Frankfurt am Main, Germany) and a research grant from the Ruhr-University of Bochum (FoRUM grant: AZ 713-11; Bochum, Germany).
- Adrian RJ: Particle-imaging technique for experimental fluid mechanics. Annu Rev Fluid Mech. 1991, 23: 261-304. 10.1146/annurev.fl.23.010191.001401.View ArticleGoogle Scholar
- Kim HB, Hertzberg JR, Shandas R: Development and validation of echo PIV. Experiments in Fluids. 2004, 36: 455-62. 10.1007/s00348-003-0743-5.View ArticleGoogle Scholar
- Zheng H, Mukdadi O, Hertzberg J, Shandas R: Advantages in using multi-frequency driving ultrasound for optimizing echo particle image velocimetry techniques. Biomed Sci Instrum. 2004, 40: 371-6.PubMedGoogle Scholar
- Hong GR, Pedrizzetti G, Tonti G, Li P, Kim JK, Wei Z, Kim JK, Baweja A, Liu S, Chung N, Houle H, Narula J, Vannan MA: Characterization and quantification of vortex flow in the human left ventricle by contrast echocardiography using vector particle image velocimetry. J Am Coll Cardiol Img. 2008, 1: 705-17.View ArticleGoogle Scholar
- Faludi R, Szulik M, D’hooge J, Herijgers P, Rademakers F, Pedrizzetti G: Voigt JU. Left ventricular flow patterns in healthy subjects and patients with prosthetic mitral valves. An in-vivo-study using Echocardiographic Particle Image Velocimetry. J Thorac Cardiovasc Surg. 2010, 139: 1501-10.PubMedGoogle Scholar
- Kheradvar A, Houle H, Pedrizzetti G, Tonti G, Belcik T, Ashraf M, Linder JR, Gharib M, Sahn D: Echocardiographic particle image velocimetry: a novel technique for quantification of left ventricular blood vorticity pattern. J Am Soc Echocardiogr. 2010, 23: 86-94. 10.1016/j.echo.2009.09.007View ArticlePubMedGoogle Scholar
- Gao H, Claus P, Amzulescu MS, Stankovic I, D'hooge J, Voigt JU: How to optimize intracardiac blood flow tracking by echocardiographic particle image velocimetry?. Eur J Echocardiogr, Exploring the influence of data acquisition using computer-generated data sets, 2011, [Epub ahead of print].Google Scholar
- Walker AR, Phillips DJ, Powers JE: Evaluating Doppler devices using a moving string test target. J Clin Ultrasound. 1982, 10: 25-30. 10.1002/jcu.1870100106View ArticlePubMedGoogle Scholar
- Walker A, Olsson E, Wranne B, Ringqvist I, Ask P: Accuracy of spectral Doppler flow and tissue velocity measurements in ultrasound systems. Ultrasound Med Biol. 2004, 30: 127-32. 10.1016/j.ultrasmedbio.2003.08.020View ArticlePubMedGoogle Scholar
- Smith SW, Insana MF, Lopez H: New contrast-detail phantoms for improved precision in lesion detection measurements. Ultrasound Med Biol. 1989, 15: 383-93. 10.1016/0301-5629(89)90050-1View ArticlePubMedGoogle Scholar
- Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986, 8476: 307-10.View ArticleGoogle Scholar
- Wigström L, Ebbers T, Fyrenius A, Karlsson M, Engvall J, Wranne B, Bolger AF: Particle trace visualization of intracardiac flow using time-resolved 3D phase contrast MRI. Magn Reson Med. 1999, 41: 793-9. 10.1002/(SICI)1522-2594(199904)41:4<793::AID-MRM19>3.0.CO;2-2View ArticlePubMedGoogle Scholar
- Kim WY, Walker PG, Pedersen EM, Poulsen JK, Oyre S, Houlind K, Yoganathan AP: Left ventricular blood flow patterns in normal subjects: a quatitative analysis by three-dimensional magnetic resonance velocity mapping. J Am Coll Cardiol. 1995, 26: 224-38. 10.1016/0735-1097(95)00141-LView ArticlePubMedGoogle Scholar
- Yang GZ, Mohiaddin RH, Kilner PJ, Firmin DN: Vortical flow feature recognition: a topological study of in vivo flow patterns using MR velocity mapping. J Comput Assist Tomogr. 1998, 22: 577-86. 10.1097/00004728-199807000-00013View ArticlePubMedGoogle Scholar
- Kvitting JP, Ebbers T, Wigström L, Engvall J, Olin CL, Bolger AF: Flow patterns in the aortic root and the aorta studied with time-resolved, 3-dimensional, phase contrast magnetic resonance imaging: implications for aortic valve-sparing surgery. J Thorac Cardiovasc Surg. 2004, 127: 1602-7. 10.1016/j.jtcvs.2003.10.042View ArticlePubMedGoogle Scholar
- Kilner PJ, Yang GZ, Wilkes AJ, Mohiaddin RH, Firmin DN, Yacoub MH: Asymmetric redirection of flow through the heart. Nature. 2000, 404: 759-61. 10.1038/35008075View ArticlePubMedGoogle Scholar
- Mathieu J: Scott J. 2000, An introduction to turbulent flow, Cambridge University Press.Google Scholar
- Fyrenius A, Wigström L, Ebbers T, Karlsson M, Engvall J, Bolger AF: Three dimensional flow in the human left atrium. Heart. 2001, 86: 448-55. 10.1136/heart.86.4.448View ArticlePubMedPubMed CentralGoogle Scholar
- Bolzon G, Zovatto L, Pedrizetti G: Birth of three-dimensionality in a pulsed jet through a circular orifice. J Fluid Mech. 2003, 493: 209-18.View ArticleGoogle Scholar
- Domenichini F, Pedrizetti G, Baccani B: Three-dimensional filling flow into a model left ventricle. J Fluid Mech. 2005, 539: 179-98. 10.1017/S0022112005005550.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.