- Workshop presentation
- Open Access
Accuracy of MRI wall shear stress estimation
© Petersson et al; licensee BioMed Central Ltd. 2012
- Published: 1 February 2012
- Wall Shear Stress
- Point Spread Function
- High Wall Shear Stress
- Voxel Position
- Shear Stress Estimation
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.
Linear regression results for the entire interval (1-20 N/m2).
Voxel Size [mm]
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.
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.