- Workshop presentation
- Open Access
Comparison of divergence-free algorithms for 3D MRI with three-directional velocity encoding
Journal of Cardiovascular Magnetic Resonance volume 14, Article number: W64 (2012)
Two divergence-free algorithms (finite difference and radial basis functions) are evaluated for improving the flow visualization quality of PC-VIPR 4D flow data. Both methods improve streamline length, with the finite difference method performing slightly better.
The velocity-to-noise ratio in 4D MR Flow acquisitions can suffer, particularly for low velocity segments, because the velocity encoding settings are adjusted to the highest velocities in the volume. Blood is essentially incompressible so vascular velocity fields should be divergence free. Therefore, any remaining divergence in an acquired 3D or 4D MR flow dataset can be assumed to originate from image noise, which can be reduced by imposing divergence free conditions. Several methods have been proposed to reconstruct a divergence free velocity vector field, yet to our knowledge no comparison studies exist to explore their properties in flow MRI data. This study compares two methods in their effectiveness at improving flow visualizations, the finite difference method (FDM) and the radial basis function (RBF) based approach, which differ in their handling of noise and boundary conditions.
MR flow data were acquired in a straight tube flow phantom and the internal carotid arteries of a volunteer on a clinical 3T scanner with a 3D radially undersampled acquistion (PC-VIPR). Time resolved PC-VIPR images were acquired as follows: scan time: 6 min acquired isotropic spatial resolution: 0.68 mm, TR/TE = 8.3/2.9 ms, temporal resolution: 39.5 ms. The FDM algorithm is solved with a fast Fourier solver for the Poisson equation. The RBF algorithm includes a normalized convolution, and is solved using iterative least squares. The quality of the resulting flow fields was quantified by average streamline lengths and streamline counts that extended from the emitter plane to a distal portion in the flow phantom or vessel. Data fidelity was measured by comparing flow values pre and post processing in the emitter and analysis planes.
Percentage of streamlines traveling 10cm in the flow phantom: original dataset = 9.2%, FDM = 59.4%, RBF = 40.4%. Average length of streamlines in the flow phantom: original dataset = 6.29cm, FDM = 10.17cm, RBF = 8.10cm. Percentage of streamlines traveling approx. 6 cm along the internal carotid artery: original dataset = 7.4%, FDM = 26.2%, RBF = 12.5%. Average length of streamlines in volunteer: original dataset = 4.41cm, FDM = 5.67cm, RBF = 5.44cm. Variation in flow values from uncorrected data: FDM = 0.1±5.9%, RBF = by 0.6±0.4%. Time to compute: FDM = ~2 min, RBF = ~60 min.
Both methods, FDM and RBF based, show a significant increase in the streamline length and the number of streamlines that reach a distal portion of the vessel. In our comparisons, the FDM algorithm performed better and was computed significantly faster than the iterative RBF method, whereas the latter introduced fewer variations in flow values. We believe that divergence free methods could become an important component of 4D-MRI flow reconstructions to improve data consistency and visualization.
NIH grant R01HL072260.
Song SM: JMRI. 1993, 587-96.
Busch J: ISMRM. 2011, #1201-
Johnson KM: MRM. 2008, 1329-36.
About this article
Cite this article
Loecher, M., Kecskemeti, S., Turski, P. et al. Comparison of divergence-free algorithms for 3D MRI with three-directional velocity encoding. J Cardiovasc Magn Reson 14 (Suppl 1), W64 (2012). https://doi.org/10.1186/1532-429X-14-S1-W64
- Internal Carotid Artery
- Radial Basis Function
- Original Dataset
- Finite Difference Method
- Velocity Encode