Accuracy of MRI wall shear stress estimation
Journal of Cardiovascular Magnetic Resonance volume 14, Article number: W6 (2012)
Wall shear stress (WSS) is a potential biomarker for vascular disease. The aim of this work was to investigate the accuracy of WSS estimation using MRI. The influence of spatial resolution, wall segmentation and voxel location were investigated over a range of WSS values by using numerical simulations.
Three methods for WSS estimation were studied. These methods are based on 1) linear extrapolation (LE) of MRI velocity data, 2) MRI velocity data in combination with estimation of location of vessel wall, and 3) Fourier velocity encoding (FVE).
Numerical velocity fields representing axisymmetric 2D velocity profiles were generated for WSS values ranging from 1-20 N/m2. Based on the numerical velocity fields, phase-contrast MRI data voxels were simulated as follows: A jinc-function was used to model the 2D point spread function (PSF), and this PSF was used to obtain each voxel’s intravoxel velocity distribution. The phase-contrast MRI signal of each voxel was simulated by taking the Fourier transform of this distribution. To account for the fact that voxels cannot be positioned exactly at the wall in an MR-experiment, all simulations were carried out for ten different voxel positions uniformly distributed over one voxel length. In the LE method the spatial velocity derivative was estimated as the velocity difference between the two adjacent near-wall voxels divided by the distance between them. In the wall-based method, WSS was estimated by dividing the linear interpolated velocity at one voxel distance from the wall by the distance to the wall. Errors in segmentations of wall position were accounted for by modeling them as normally distributed with a standard deviation of 1/4 voxel size. In the FVE-based method, the WSS was obtained by first estimating the intravoxel velocity profile via a simulated FVE measurement and then computing the spatial velocity derivative near the wall. Note that the FVE-method uses larger voxels.
All methods were sensitive to spatial resolution, especially for high WSS (Fig. 1). The velocity-based methods generally underestimated the WSS, and were unable to resolve the highest WSS values (Fig. 1 a,c). The FVE-method was most sensitive to the voxel position (Fig. 1c). The wall-based method was sensitive to errors in segmentation (Fig 1 b). Linear regression results for all three methods are shown in Table 1.
MRI WSS estimation is hampered by important limitations that require attention in studies where it is applied. WSS obtained from MR velocity data can be substantially different from WSS obtained by FVE. Although WSS was underestimated and influenced by parameter settings and segmentation errors, distinguishing areas of low and elevated WSS may be feasible.
Swedish Research Council, Swedish Heart-Lung Foundation.
About this article
Cite this article
Petersson, S., Dyverfeldt, P. & Ebbers, T. Accuracy of MRI wall shear stress estimation. J Cardiovasc Magn Reson 14 (Suppl 1), W6 (2012). https://doi.org/10.1186/1532-429X-14-S1-W6