- Research article
- Open Access
Volume Tracking: A new method for quantitative assessment and visualization of intracardiac blood flow from three-dimensional, time-resolved, three-component magnetic resonance velocity mapping
BMC Medical Imaging volume 11, Article number: 10 (2011)
Functional and morphological changes of the heart influence blood flow patterns. Therefore, flow patterns may carry diagnostic and prognostic information. Three-dimensional, time-resolved, three-directional phase contrast cardiovascular magnetic resonance (4D PC-CMR) can image flow patterns with unique detail, and using new flow visualization methods may lead to new insights. The aim of this study is to present and validate a novel visualization method with a quantitative potential for blood flow from 4D PC-CMR, called Volume Tracking, and investigate if Volume Tracking complements particle tracing, the most common visualization method used today.
Eight healthy volunteers and one patient with a large apical left ventricular aneurysm underwent 4D PC-CMR flow imaging of the whole heart. Volume Tracking and particle tracing visualizations were compared visually side-by-side in a visualization software package. To validate Volume Tracking, the number of particle traces that agreed with the Volume Tracking visualizations was counted and expressed as a percentage of total released particles in mid-diastole and end-diastole respectively. Two independent observers described blood flow patterns in the left ventricle using Volume Tracking visualizations.
Volume Tracking was feasible in all eight healthy volunteers and in the patient. Visually, Volume Tracking and particle tracing are complementary methods, showing different aspects of the flow. When validated against particle tracing, on average 90.5% and 87.8% of the particles agreed with the Volume Tracking surface in mid-diastole and end-diastole respectively. Inflow patterns in the left ventricle varied between the subjects, with excellent agreement between observers. The left ventricular inflow pattern in the patient differed from the healthy subjects.
Volume Tracking is a new visualization method for blood flow measured by 4D PC-CMR. Volume Tracking complements and provides incremental information compared to particle tracing that may lead to a better understanding of blood flow and may improve diagnosis and prognosis of cardiovascular diseases.
Blood flow patterns are closely linked to the function of valves, vessels and the myocardium, and functional and morphological changes in the heart influence the flow patterns [1–4]. Thus a better understanding of flow patterns may lead to improved diagnosis and prognosis of cardiovascular diseases.
Cardiovascular Magnetic Resonance (CMR) enables the measurement of the full three-dimensional and time-resolved velocity field in the heart and great vessels using four-dimensional Phase Contrast Cardiovascular Magnetic Resonance (4D PC-CMR) [5–7]. Due to the complexity of 4D flow fields, there exists no single natural representation that can show all details of the flow. To handle this abundance of information, flow visualization is often used for discovery, understanding and communication about flow . Visualization is a critical part of the discovery process, as it enables the researcher to interact with the complex data, and by using the powerful pattern recognition systems of the human visual system, new hypotheses about the flow can be formed and then be investigated scientifically. In the case of three-dimensional, time-resolved (4D) flow, visualization is not only beneficial, but a necessity, since it would be an extremely time-consuming and error-prone task to interpret the flow using two-dimensional slices of the velocity volume .
Several flow visualization methods have been developed, each showing different aspects of the flow [10, 11]. The most commonly used method for visualization of 4D PC-CMR data is particle tracing [9, 12], but streamlines, vector plots  and vortex cores  have also been used. Methods for quantitative analysis have also been presented, including retrospective flow quantification , relative pressure mapping , turbulence mapping , wall shear stress , separation of left ventricular flow into components  and quantification of helical flow .
Visualization plays an important role in generating hypotheses that guide research. The result and performance of flow visualizations has been shown to be heavily influenced by the choice of flow visualization method  and even by minor adjustments in the visualization methods themselves . This means that the impression of the flow changes with the visualization method, and therefore also the hypotheses generated and subsequently the studies performed and conclusions drawn. Since a single method may not show all aspects of the flow, research in alternative methods for blood flow visualization is motivated.
We have therefore developed a new visualization method, called Volume Tracking. Volume Tracking is similar to particle tracing in that it displays the path taken by the blood through the circulation, but follows and displays volumes and surfaces instead of points, changing the visual representation of the flow. Additionally, Volume Tracking adds a quantitative component to the visualizations. The purposes of this article are 1. To present and validate a newly developed flow visualization method, Volume Tracking, and 2. To investigate if Volume Tracking complements particle tracing and provides incremental information.
Study population and CMR protocol
Eight healthy volunteers were included in the study. All volunteers had normal blood pressure, normal ECG and no history of cardiovascular disease. One patient was also included to investigate if Volume Tracking was feasible in patients. This patient was hypertensive and had an apical aneurysm that developed after a major anterior infarction, resulting in an LV ejection fraction of 15%. A standard CMR examination consisting of Cine images in the 2-chamber, 3-chamber, 4-chamber and short-axis views was performed on all subjects. An experienced observer verified normal anatomy, myocardial function and valve function in the healthy volunteers. The study was approved by the Regional Ethical Review Board in Lund, Sweden. Written informed consent was obtained from all subjects. Subject characteristics are summarized in Table 1.
Three-dimensional, time-resolved, three-component velocity mapping (4D PC-CMR) was performed over the whole heart on all subjects using a three-dimensional Turbo Field Echo (TFE) Phase Contrast sequence  on a Philips Achieva 3T CMR scanner. The 4D flow and cine images were acquired in the same imaging session, so that the cine images could be used as spatial orientation when visualizing the 4D flow data. For the 4D flow measurements, the following parameters were used: Spatial resolution 3 × 3 × 3 mm, flip angle 8°, TE 3.7 ms, TR 6.2 ms, velocity encoding (VENC) 100 cm/s, and temporal resolution typically 50 ms. Parallel imaging (SENSE) was used with a factor of 2, and the turbo factor was 2. Respiratory gating to the end-expiratory phase was performed with a navigator beam for the diaphragm. Synchronization to the cardiac cycle was performed using retrospective ECG triggering, covering the whole R-R interval. All imaging was performed with the subjects at rest. Scan time for the flow measurement was typically 35 minutes. The flow measurement and standard CMR examination were performed in the same session, bringing the total scan time to about 50 minutes.
This section presents an overview of the Volume Tracking theory and implementation. The full mathematical formulation and technical details of Volume Tracking, necessary to reproduce and analyze the method, are given in Additional File 1: VT-Appendix.pdf.
Volume Tracking is based on the concept of a flow map , which is a function that maps the location of a particle at one time t 0 to its position at another time t 1, thereby answering the question "Where will this particle go?" or "Where did this particle come from?". Using the flow map, the evolution of any volume in the flow can be defined: If the volume at the starting time t 0 is V 0, and the volume at another time t is V(t), we define V(t) as the image of V 0 under the flow map as
To compute the flow map, an auxiliary function is introduced, mapping particles from their current position to their positions at time t 0, the starting time. Denoting the velocity field by , and decomposing ψ into its cartesian components (ψ x , ψ y , ψ z ), and decomposing = (x, y, z), the advection equations
then describe the motion of volumes in the flow. The initial conditions at t = t 0 are
These equations can be efficiently solved using a finite volume method, further described Additional File 1: VT-Appendix.pdf.
A starting volume V 0 is specified using a volume function f:
and the evolution of a volume can then be computed as
The solution to the Volume Tracking equations can be pre-computed, and the evolution of any volume is encoded in the flow map. Blood volumes of different shape, such as spheres and planes, can then be followed forward or backward in time interactively. Using the current implementation of Volume Tracking, the pre-computation step typically lasts 10-15 minutes per starting time. However, the pre-computation step can be performed offline without human interaction once the starting times have been chosen. Once the pre-computation has been completed, the flow map can be used to select volumes of several shapes and any position, and the movement and deformation of this volume can be displayed instantaneously.
The full mathematical formulation and technical details of Volume Tracking, necessary to reproduce and analyze the method, are given in Additional File 1: VT-Appendix.pdf.
Volume Tracking was implemented in MATLAB (The Mathworks, USA), using the freely available Finite Volume solver CLAWPACK [22–24] for the numerical calculations. User interaction, image display and animation were performed in the visualization software package Ensight (CEI, USA).
To correct for eddy currents and other phase background effects, a first-order polynomial was fitted to velocities in stationary tissue and subtracted from the acquired velocity field . Additionally, automatic phase unwrapping was performed to correct for phase aliasing effects . Concomitant gradient effects were compensated by the CMR scanner software. Correction of background effects and phase unwrapping of the velocity data was implemented as a plugin to Segment , a freely available software for cardiovascular image analysis based on MATLAB (The Mathworks, USA). The visualization software package Ensight 9.1 (CEI, USA) was used to display Volume Tracking surfaces, compute particle traces and to create images and animations. Particle traces were computed using Ensight's default fourth-order Runge-Kutta integration algorithm , which uses linear interpolation between timesteps and grid points. The implementation of Volume Tracking is described in Additional File 1: VT-Appendix.pdf.
Evaluation and validation of Volume Tracking
Particle tracing and Volume Tracking both display the path taken by blood through the cardiovascular system, although giving different visual impressions of the flow. Therefore, to evaluate the visualization aspects of Volume Tracking, Volume Tracking and particle trace visualizations were generated in the visualization software package Ensight (CEI, USA) and compared side-by-side. Side-by-side visualizations were generated for diastolic inflow in left ventricle and right ventricle in a healthy volunteer, and diastolic inflow in the left ventricle in the ischemic cardiomyopathy patient. In the visualizations of left ventricular flow, the Volume Tracking surface shows the boundary between inflowing blood and blood already in the ventricle from the previous heartbeat.
To validate Volume Tracking against particle tracing, combined visualizations were generated of LV inflow in all subjects, with Volume Tracking and particle tracing displayed superimposed in the same image (see Figure 1 and Additional File 2: VT-PT-Combination.mpg). Since both Volume Tracking and particle tracing compute and display the motion of blood through the heart, no particles should cross the Volume Tracking surface if the methods display the exact same blood motion. Therefore, a particle was defined as fulfilling the validation criterion if it was located on the same side of the Volume Tracking plane as it was released on. The agreement between Volume Tracking and particle tracing at mid-diastole and end-diastole was calculated as the number of particles fulfilling the validation criterion, divided by the number of particles released up to that timepoint, and expressed as a percentage. At least 1000 particles were released in each subject. To ensure that this validation methodology was not sensitive to the number of released particles, the number of particles was varied between 1000 and 8000 in subject 7.
Analysis of LV inflow patterns
To demonstrate one of the possible applications of Volume Tracking, Volume Tracking visualizations were used to describe blood flow patterns in the left ventricle (LV). Using a Volume Tracking plane, as shown in Figure 2, panels D-E, the boundary between blood flowing into the ventricle and blood already in the ventricle was determined for the whole diastolic period. To map the location of the inflow blood, the LV lumen was divided into segments according to Figure 3. The LV lumen was divided into Apical, Mid-ventricular and Basal parts along the long-axis (Figure 3, panel A). Each part was then divided into seven segments (Figure 3, panel B) adjusted from the AHA model . The intersection between the Volume Tracking surface and short-axis slices at the Apical, Mid-ventricular and Basal levels was then visualized (Figure 3, panel C). Two independent observers estimated whether each segment contained more or less than 25% blood. The location of the inflow blood was determined in this way at two times during diastole: at mid-diastole, defined as halfway between peak flow during the E-wave and peak flow during the A-wave, and at end-diastole, defined as the time just before systolic contraction. The E-wave and A-wave were measured in the mitral annulus in the 4D PC-CMR flow data and used for timing reference only. Agreement between the observers was measured using the kappa coefficient .
When using Volume Tracking, a starting volume of blood is first chosen. In Figure 2, Panel D the blood in the left atrium has been chosen, and the surface of this volume is displayed. As the blood flows, the volume will move and deform with the flow as if the volume's surface was infinitely flexible and stretchable, as seen in Figure 2, panel D-F, and Additional File 3: VT-Comparison.mpg. The Volume Tracking surface shows the boundary between inflowing blood and blood already in the ventricle from the previous heartbeat.
Although only one shape is used in this example, a wide variety of sizes and shapes can be used. Planes are used in Figures 1, 2 and 4, and Figure 5 shows a sphere as initial volume. A more detailed discussion of the possible shapes can be found in Additional File 1: VT-Appendix.pdf. The formulation of the method allows the volume selection to be performed interactively using a point-and-click user interface.
Evaluation and validation of Volume Tracking
In Figure 2 a comparison of visualizations of LV inflow using Volume Tracking and particle tracing is shown. Additional File 3: VT-Comparison.mpg shows an animated version of Figure 2, showing the effect of animation on Volume Tracking. Particle tracing shows an intricate pattern of swirling flows during diastolic filling, especially apparent in the animation (Additional File 3: VT-Comparison.mpg). Volume Tracking, in contrast, shows the boundary between the inflowing blood and the blood already in the ventricle. Additionally, Volume Tracking reveals an area close to the basal parts of the septum with no filling blood in subject 7 (Figure 2, lower right panel). This area was not observed when studying the particle trace visualization alone. Particle tracing shows a dynamic, complex scene, while Volume Tracking shows the combined movement of the whole inflow volume.
Figure 4 and Additional File 4: VT-Patient.mpg shows inflow into the left ventricle (LV) in a patient with a LV apical aneurysm, visualized using Volume Tracking and particle tracing. As in Figure 2, the visual impression of the LV inflow is different for the two methods. The inflow is distinctly different from the flow pattern in the healthy volunteer in Figure 2.
Figure 5 shows particle trace and Volume Tracking visualizations of inflow in the right ventricle (RV). Additional File 5: VT-RV.mpg shows an animated version of Figure 5. To demonstrate that Volume Tracking can be used with many different shapes, a sphere of blood in the right atrium (panel D) is used instead of a plane. The particle tracing visualization shows a vortex or helix motion of particles into the ventricle, and the particles coming to rest in the RV. Volume Tracking shows the initial volume deforming and moving into the right ventricle, but the vortex motion is not clearly visible. As the blood comes to rest in the RV during diastasis, it spreads out in a complex pattern. This complex pattern is not apparent in the particle tracing visualization. Additional File 5: VT-RV.mpg clearly shows the breakdown of ordered, vortical inflow to the complex deformation in late diastole.
Figure 1 shows the superposition of Volume Tracking and particle tracing visualizations of LV filling in a healthy volunteer. Additional File 2: VT-PT-Combination.mpg shows an animated version of Figure 1. Notice that very few particles pass through the plane. Table 2 shows the agreement between Volume Tracking and particle tracing, i.e. the fraction of particles that stay on the inflow side of the Volume Tracking surface. Mean agreement for all 9 subjects was 90.5% in mid-diastole and 87.8% in end-diastole. When varying the number of particles in subject 7 between 1000 and 8000, the standard deviation in the mid-diastolic agreement was 0.30% and in the end-diastolic agreement 0.17%.
Analysis of LV inflow patterns
Semi-quantitative analysis of the inflow patterns was feasible in all subjects. The results of the flow pattern analysis are given in Tables 3 and 4 for Observer 1, and in Additional File 6: Distribution of LV inflow blood.xls. Agreement between the observers was excellent with a kappa coefficient of κ = 0.91. In seven out of eight volunteers an in the patient, blood did not reach the apical level by mid-diastole. In end-diastole, blood has reached the apical level in the majority of healthy subjects (observer 1: 6 of 8, observer 2: 8 of 8). In the patient, the inflowing blood did not reach the apical level at all. There was a large variability inflow patterns between the healthy volunteers. In the patient with ischemic cardiomyopathy, the inflow blood is only present in segments 4, 5, 6 and 7 in the basal third of the ventricle. This is notably different from the normals, where blood is present in several segments in the basal and mid-ventricular thirds.
In this article, a new method for visualization and analysis of 4D PC-CMR data, called Volume Tracking, is presented for the first time. Visually, Volume Tracking provides a different perspective on the flow data when compared to particle tracing. Volume Tracking and particle tracing have a mean agreement of 90.5% in mid-diastole and 87.8% in end-diastole. In the analysis of LV inflow patterns, observer agreement was excellent. The flow patterns varied between the normals, and in the ischemic cardiomyopathy patient, the inflow was notably different from the normals.
Evaluation and validation of Volume Tracking
Three examples demonstrate that Volume Tracking provides incremental value to particle tracing: a) Volume Tracking shows the surface of the blood volume flowing into the left ventricle rather than the detailed dynamics; b) Volume Tracking emphasizes overall motion more than particle tracing; and c) Volume Tracking shows a complex flow pattern in RV filling not visible using particle tracing. This is in line with previous studies, where the visualization method influences the results of the visualization process [20, 21]. However, this is the first study investigating the influence of the choice of visualization method for 4D PC-CMR data. The significance of this result may be that when visualizing and analyzing 4D PC-CMR data, researchers can extract more information from a given dataset using Volume Tracking. This may lead to new insights and additional knowledge about blood flow and cardiac pumping.
In general, it may be beneficial to use several complementary visualization methods when studying 4D PC-CMR data, since it may reveal more aspects of the flow. One example of this is that the vortex in the right ventricular inflow was less apparent with Volume Tracking than with particle tracing. To detect this type of flow patterns, vortex detection methods  and quantitative indices of helix flow [19, 31] may be effective complements to Volume Tracking. Complementary to the large-scale flow dynamics shown by Volume Tracking, the flow dynamics in vessel wall boundary layers can be quantified using wall shear stress (WSS)  analysis. Simultaneous use of WSS analysis and Volume Tracking may provide a link between WSS and large-scale flow. Furthermore, the visual impression of the flow gained by Volume Tracking can be quantified by classifying particle traces according to origin and destination of flow [18, 32].
The high agreement between Volume Tracking and particle tracing shows that Volume Tracking displays the same blood motion as the particle traces, which in this study were computed using a numerical algorithm independent from Volume Tracking and known to be accurate . This suggests that the Volume Tracking numerical algorithm performs well in the presence of noise, acceleration and vortices in the 4D PC-CMR datasets used in this study. Additionally, the agreement between Volume Tracking and particle tracing shows that the differences in interpretation of the flow is not due to differences in blood transport calculations. The differences in interpretation of the flow may instead be explained by the differences in visual representation (surfaces instead of particles and lines). The low variation in agreement when changing the number of particles suggests that the agreement is not dependent on the number of particles.
In this study, only ventricular inflow has been studied. However, Volume Tracking is a general method, which may be used to study flow in all parts of the cardiovascular system. One example where Volume Tracking may also provide incremental value is the human aorta, where flow patterns may be influential in plaque formation and rupture.
Analysis of LV inflow patterns
The physiological significance of the present findings may be an increased understanding of the complex, and clinically intangible, hemodynamics that may help clarify normal cardiac pumping. The different filling patterns in the healthy volunteers suggest that normal hemodynamics may display additional degrees of freedom that could explain differences in normal physiology.
The pathophysiological significance of the present findings may be that the transition from normal to pathological LV hemodynamics which may take several years to develop can be better understood and discovered at an earlier sub-clinical stage. These gradual changes over time may also be coupled to early morphological and physiological changes. In line with previous studies , we agree that flow patterns may be a sensitive measure of LV function and may be able to detect early stages of dysfunction.
Errors in the velocity measurements accumulate in particle paths and Volume Tracking surfaces over time . This problem is minimized in this study by limiting the tracking to diastole.
Blood flow patterns may be influenced by factors which were not measured in the present experimental setting, such as changes in blood viscosity and hematocrit.
Long 4D PC-CMR acquisition times lead to averaging over several heartbeats, obscuring small flow details so that the resulting flow data describes large-scale flow features that are present in the majority of heartbeats. The 4D PC-CMR technique is also limited in spatial resolution, does not provide absolute pressure levels, and is not feasible during exercise. Noise in the 4D PC-CMR measurements and the complicated flow patterns in the human heart may also influence the Volume Tracking calcuations. The high agreement between particle tracing and Volume Tracking suggests that accuracy in the Volume Tracking algorithm is a minor problem for the present application.
Volume Tracking shows how a certain volume is deformed in the flow. This means that flow inside the volume is only seen indirectly. To overcome this limitation, Volume Tracking and particle tracing can be used simultaneously, as in Figure 1.
Ensight uses linear interpolation of the velocity fields in space and time to compute particle traces, which may introduce errors into the calculation. Since Ensight does not allow the user to choose interpolation method, the influence of the interpolation could not be checked. The authors are not aware of any studies of the effects of interpolation strategies in particle tracing, but the effects are probably small in relation to the smoothing, noise, low resolution and phase background effects of the 4D PC-CMR technique.
The evaluation of the visualization aspects of Volume Tracking was performed by visual comparison to particle tracing. An independent, quantitative evaluation for both methods would be desirable. However, the state of the art in evaluation methods for flow visualization has not yet been extended to treat time-resolved flow [20, 21] and may not adress the hypothesis-generating aspects of visualization.
In this study Volume Tracking, a new method for visualization of blood flow from 4D PC-CMR, is presented. Compared to particle tracing, Volume Tracking can provide incremental information that may lead to a better understanding of blood flow and may improve diagnosis and prognosis of cardiovascular diseases. The diastolic inflow pattern in the patient is notably different from the volunteers, suggesting that blood flow patterns may be used as a measure of LV function and as an early indicator of disease.
Beppu S, Izumi S, Miyatake K, Nagata S, Park Y, Sakakibara H, Nimura Y: Abnormal blood pathways in left ventricular cavity in acute myocardial infarction. Experimental observations with special reference to regional wall motion abnormality and hemostasis. Circulation. 1988, 78 (1): 157-
Reiter G, Reiter U, Kovacs G, Kainz B, Schmidt K, Maier R, Olschewski H, Rienmueller R: Magnetic Resonance-Derived 3-Dimensional Blood Flow Patterns in the Main Pulmonary Artery as a Marker of Pulmonary Hypertension and a Measure of Elevated Mean Pulmonary Arterial Pressure. Circulation: Cardiovascular Imaging. 2008, 1 (1): 23-10.1161/CIRCIMAGING.108.780247.
Thanigaraj S, Chugh R, Schechtman K, Lee L, Wade R, Pérez J: Defining left ventricular segmental and global function by echocardiographic intraventricular contrast flow patterns. The American journal of cardiology. 2000, 85 (1): 65-68. 10.1016/S0002-9149(99)00608-6.
Garcia M, Smedira N, Greenberg N, Main M, Firstenberg M, Odabashian J, Thomas J: Color M-mode Doppler flow propagation velocity is a preload insensitive index of left ventricular relaxation: animal and human validation. Journal of the American College of Cardiology. 2000, 35 (1): 201-208. 10.1016/S0735-1097(99)00503-3.
Wigström L, Sjövist L, Wranne B: Temporally resolved 3D phase-contrast imaging. Magn Reson Med. 1996, 36 (5): 800-3. 10.1002/mrm.1910360521.
Dyverfeldt P, Kvitting J, Sigfridsson A, Engvall J, Bolger A, Ebbers T: Assessment of fluctuating velocities in disturbed cardiovascular blood flow: In vivo feasibility of generalized phase-contrast MRI. Journal of Magnetic Resonance Imaging. 2008, 28 (3): 655-663. 10.1002/jmri.21475.
Markl M, Kilner P, Ebbers T: Comprehensive 4D velocity mapping of the heart and great vessels by cardiovascular magnetic resonance. Journal of Cardiovascular Magnetic Resonance. 2011, 13: 7-10.1186/1532-429X-13-7.
DeFanti T, Brown M, McCormick B: Visualization: Expanding Scientific and Engineering Research Opportunities. Computer. 1989, 22 (8): 12-25. 10.1109/2.35195.
Napel S, Lee D, Frayne R, Rutt B: Visualizing three-dimensional flow with simulated streamlines and three-dimensional phase-contrast MR imaging. J Magn Reson Imaging. 1992, 2 (2): 143-153. 10.1002/jmri.1880020206.
McLoughlin T, Laramee R, Peikert R, Post F, Chen M: Over two decades of integration-based, geometric flow visualization. Computer Graphics Forum. 2010, 29 (6): 1807-1829. 10.1111/j.1467-8659.2010.01650.x.
Post FH, Vrolijk B, Hauser H, Laramee RS, Doleisch H: Feature Extraction and Visualisation of Flow Fields. Proceedings of the Eurographics Conference 2002. Edited by: Fellner D, Scopigno R. 2002, Saarbrücken, Germany: Eurographics Association
Wigström L, Ebbers T, Fyrenius A, Karlsson M, Engvall J, Wranne B, Bolger A: Particle Trace Visualization of Intracardiac Flow Using Time-Resolved 3D Phase Contrast MRI. Magnetic Resonance in Medicine. 1999, 41 (4): 793-799. 10.1002/(SICI)1522-2594(199904)41:4<793::AID-MRM19>3.0.CO;2-2.
Buonocore M: Visualizing blood flow patterns using streamlines, arrows, and particle paths. Magnetic resonance in medicine. 1998, 40 (2): 210-226. 10.1002/mrm.1910400207.
Heiberg E, Ebbers T, Wigstrom L, Karlsson M: Three-dimensional flow characterization using vector pattern matching. Visualization and Computer Graphics, IEEE Transactions on. 2003, 9 (3): 313-319. 10.1109/TVCG.2003.1207439.
Westenberg J, Roes S, Ajmone Marsan N, Binnendijk N, Doornbos J, Bax J, Reiber J, de Roos A, van der Geest R: Mitral Valve and Tricuspid Valve Blood Flow: Accurate Quantification with 3D Velocity-encoded MR Imaging with Retrospective Valve Tracking. Radiology. 2008, 249 (3): 792-10.1148/radiol.2492080146.
Ebbers T, Wigström L, Bolger A, Wranne B, Karlsson M: Noninvasive Measurement of Time-Varying Three-Dimensional Relative Pressure Fields Within the Human Heart. Journal of Biomechanical Engineering. 2002, 124 (3): 288-10.1115/1.1468866.
Stalder A, Russe M, Frydrychowicz A, Bock J, Hennig J, Markl M: Quantitative 2D and 3D phase contrast MRI: optimized analysis of blood flow and vessel wall parameters. Magnetic Resonance in Medicine. 2008, 60 (5): 1218-1231. 10.1002/mrm.21778.
Eriksson J, Carlhäll C, Dyverfeldt P, Engvall J, Bolger A, Ebbers T: Semi-automatic quantification of 4 D left ventricular blood flow. Journal of Cardiovascular Magnetic Resonance. 2010, 12: 9-10.1186/1532-429X-12-9.
Morbiducci U, Ponzini R, Rizzo G, Cadioli M, Esposito A, De Cobelli F, Del Maschio A, Montevecchi F, Redaelli A: In vivo quantification of helical blood flow in human aorta by time-resolved three-dimensional cine phase contrast magnetic resonance imaging. Annals of biomedical engineering. 2009, 37 (3): 516-531. 10.1007/s10439-008-9609-6.
Laidlaw D, Kirby R, Jackson C, Davidson J, Miller T, Da Silva M, Warren W, Tarr M: Comparing 2D vector field visualization methods: A user study. IEEE Transactions on Visualization and Computer Graphics. 2005, 11 (1): 59-70. 10.1109/TVCG.2005.4.
Forsberg A, Chen J, Laidlaw D: Comparing 3D Vector Field Visualization Methods: A User Study. IEEE Transactions on Visualization and Computer Graphics. 2009, 15 (6): 1219-1226.
LeVeque R, Berger M: CLAWPACK Software 4.3. [Accessed 2009-10-23], [http://www.clawpack.org/]
LeVeque R: High-Resolution Conservative Algorithms for Advection in Incompressible Flow. SIAM Journal on Numerical Analysis. 1996, 33 (2): 627-665. 10.1137/0733033.
Langseth J, LeVeque R: A Wave Propagation Method for Three-Dimensional Hyperbolic Conservation Laws. Journal of Computational Physics. 2000, 165 (1): 126-166. 10.1006/jcph.2000.6606.
Ebbers T, Haraldsson H, Dyverfeldt P, Sigfridsson A, Warntjes M, Wigström L: Higher order weighted least-squares phase offset correction for improved accuracy in phase-contrast MRI. Proceedings of the International Society for Magnetic Resonance in Medicine. 2008, 16: 1367-
Xiang Q: Temporal phase unwrapping for CINE velocity imaging. J Magn Reson Imaging. 1995, 5 (5): 529-34. 10.1002/jmri.1880050509.
Heiberg E, Sjögren J, Ugander M, Carlsson M, Engblom H, Arheden H: Design and validation of Segment-available software for cardiovascular image analysis. BMC Medical Imaging. 2010, 10: 1-10.1186/1471-2342-10-1.
Darmofal D, Haimes R: An Analysis of 3D Particle Path Integration Algorithms. Journal of Computational Physics. 1996, 123 (1): 182-195. 10.1006/jcph.1996.0015.
Cerqueira M, Weissman N, Dilsizian V, Jacobs A, Kaul S, Laskey W, Pennell D, Rumberger J, Ryan T, Verani M: Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart: a statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Circulation. 2002, 105 (4): 539-542. 10.1161/hc0402.102975.
Cohen J: A coefficient of agreement for nominal scales. Educational and psychological measurement. 1960, 20 (1): 37-10.1177/001316446002000104.
Morbiducci U, Ponzini R, Rizzo G, Cadioli M, Esposito A, Montevecchi F, Redaelli A: Mechanistic insight into the physiological relevance of helical blood flow in the human aorta: an in vivo study. Biomechanics and Modeling in Mechanobiology.
Bolger A, Heiberg E, Karlsson M, Wigstrom L, Engvall J, Sigfridsson A, Ebbers T, Kvitting J, Carlhall C, Wranne B: Transit of Blood Flow Through the Human Left Ventricle Mapped by Cardiovascular Magnetic Resonance. Journal of Cardiovascular Magnetic Resonance. 2007, 9 (5): 741-747. 10.1080/10976640701544530.
Hong G, Pedrizzetti G, Tonti G, Li P, Wei Z, Kim J, Baweja A, Liu S, Chung N, Houle H, et al: Characterization and quantification of vortex flow in the human left ventricle by contrast echocardiography using vector particle image velocimetry. JACC: Cardiovascular Imaging. 2008, 1 (6): 705-717. 10.1016/j.jcmg.2008.06.008.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/11/10/prepub
The authors would like to thank the Cardiovascular MR Research Group at Linköping University, Linköping, Sweden for help and support in our work on intracardiac blood flow.
Ann-Helen Arvidsson, Christel Carlander and Martin Ugander at the Cardiac MR Group, Department of Clinical Physiology, Lund University Hospital, Lund, Sweden assisted in data collection.
Freddy Ståhlberg and Anders Nilsson at the Department of Medical Radiation Physics, Lund University Hospital, Lund, Sweden and Karin Markenroth-Bloch at Philips Medical Systems, Best, the Netherlands provided expertise in MRI physics and data acquisition.
This study was supported by Swedish Research Council grants VR 621-2005-3129, VR 621-2008-2949 and VR K2009-65X-14599-07-3, National Visualization Program and Knowledge Foundation grant 2009-0080, the Medical Faculty at Lund University, Sweden, the Region of Scania, Sweden and the Swedish Heart-Lung Foundation.
The authors declare that they have no competing interests.
JT developed the details of the method and implemented it in software, wrote the manuscript, created all figures and animations, performed analysis of LV inflow patterns and merged feedback from the other authors. MC collected data, supplied expertise in cardiac physiology and revised the manuscript for intellectual content. GS contributed to the design and development of the numerical method. HA supplied expertise in cardiac physiology and revised the manuscript for intellectual content. EH conceived the method, performed analysis of LV inflow patterns, contributed to the method design and revised the manuscript for intellectual content. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Appendix: Theory and implementation. Appendix: Document describing the mathematical theory and technical details of the implementation of Volume Tracking. (PDF 429 KB)
Additional file 2: . Composite particle trace and Volume Tracking visualization of LV filling in volunteer 7, also shown in Figure 1. Anatomical four-chamber Cine images are displayed in the background for orientation. See Figure 1 for a description of the anatomy. The image is semi-transparent to show flow behind the four-chamber plane. Time is counted from the start of ventricular systole. At the beginning of the movie a Volume Tracking plane is visible at the level of the mitral valve, and a collection of particles are visible just above the plane. Particles are released every 20 milliseconds in the same location. The Volume Tracking plane can be seen as an infinitely flexible and stretchable sheet, deformed effortlessly with the flow. Between t = 0.337 s and t = 0.5 s, the Volume Tracking plane is deformed as blood flows into the ventricle. The particles show an ordered inflow. The Volume Tracking plane shows that the front of inflowing blood moves quickly towards the apex. After t = 0.5 s the particle traces show a complex arrangement of vortices, which slows down almost to a standstill around t = 0.8 s. The Volume Tracking surface keeps on deforming slowly even after t = 0.8s, when the particles seem to have stopped moving. Notice the correspondence between the two methods. Specifically, very few particles cross through the Volume Tracking plane. This is due to the fact that all particles are released in the inflowing blood and that the Volume Tracking plane separates the inflowing blood from the blood already in the ventricle. (MPG 2 MB)
Additional file 3: . Visualization of LV filling in volunteer 7, also shown in Figure 2. Left: particle tracing. Right: Volume Tracking. The left and right panels differ only in the visualization method used. Anatomical four-chamber Cine images are displayed in the background for orientation. See Figure 2 for a description of the anatomy. The four-chamber plane image is semi-transparent to show flow behind the plane. The color scale indicates velocities between 0 and 1.0 meters per second. Time is counted from the start of ventricular systole. The video starts just before the opening of the mitral valve. In the left panel, particles are released every 20 milliseconds. In the right panel, a plane is placed near the mitral valve. The plane is infinitely flexible and stretchable and deforms effortlessly with the flow, showing the boundary between inflow blood and blood already in the ventricle. Between t = 0.337 s and t = 0.5 s, particle tracing shows an ordered inflow, and Volume Tracking shows the inflow blood progressing quickly into the ventricle. After t = 0.5 s, particle tracing shows a complex arrangement of vortices decelerating almost to standstill around t = 0.8 s. In contrast, Volume Tracking clearly shows a clockwise rotation and further deformation in the apical region after t = 0.8s. Additionally, Volume Tracking shows that inflow blood is located closer to the lateral wall of the LV in this subject. A region near the basal parts of the septum contains no filling blood (see Figure 2, white arrow), something which is not clearly visible using particle trace visualization. (MPG 2 MB)
Animation: Particle trace and Volume Tracking visualizations of LV inflow in patient with apical aneurysm
Additional file 4: . Particle tracing (left) and Volume Tracking (right) visualizations of LV inflow in a patient diagnosed with a large LV apical aneurysm. Anatomical four-chamber Cine images are displayed in the background for orientation. See Figure 4 for a description of the anatomy. The four-chamber image is transparent to show flow behind the four-chamber plane. The color scale in the lower left corner indicates velocities between zero and 1.0 meters per second (m/s). Time is counted from the start of ventricular systole. At the start of the movie, just before the filling of the LV begins, the starting location of the particles is visible just above the mitral valve (left). A Volume Tracking plane is visible at the level of the mitral valve (right). As the movie plays, the particles flow into the ventricle. The plane deforms to show the boundary between inflowing blood and blood already in the ventricle. (MPEG 2 MB)
Additional file 5: . Particle tracing (left) and Volume Tracking (right) visualizations of RV diastolic inflow. Anatomical four-chamber Cine images are displayed in the background for orientation. See Figure 5 for a description of the anatomy. The four-chamber image is transparent to show flow behind the four-chamber plane, and the scene has been rotated in comparison to Figures 1, 2 and 4. Due to the rotation, the RV is in the foreground, with the LV behind. The color scale in the lower left corner indicates velocities between zero and 1.0 meters per second (m/s). Time is counted from the start of ventricular systole. The animation starts at the beginning of RV filling. After this the animation can be divided into two parts. First an organized inflow into the right ventricle is shown between t = 0.337 s and t = 0.63 s. Particle tracing (left) shows a distinct swirling vortex of flow into the right ventricle. Volume Tracking does not show the vortex clearly. During the remaining time, from t = 0.63 s to t = 1.00 s, the flow decelerates during diastasis. Particle tracing shows the blood coming to rest in an uncomplicated manner. However, Volume Tracking shows the blood volume becoming stretched and twisted, finally occupying most of the right ventricle in a complex structure. The volume is spread out in three dimensions, outside the image plane, but is still completely within the RA and RV. (MPEG 2 MB)
About this article
Cite this article
Töger, J., Carlsson, M., Söderlind, G. et al. Volume Tracking: A new method for quantitative assessment and visualization of intracardiac blood flow from three-dimensional, time-resolved, three-component magnetic resonance velocity mapping. BMC Med Imaging 11, 10 (2011). https://doi.org/10.1186/1471-2342-11-10
- Cardiovascular Magnetic Resonance
- Wall Shear Stress
- Right Ventricle
- Particle Trace
- Visualization Method