Vectorial wall shear stress calculations in vessel structures using 4D PC-MRI
© Potters et al; licensee BioMed Central Ltd. 2012
Published: 1 February 2012
We propose a fully automated method for calculating vectorial wall shear stress (WSS) in-vivo based on 4D PC-MRI data.
Wall shear stress (WSS) is the tangential force of flowing blood on the vessel wall. WSS directly influences remodeling of the vessel wall.
Velocity data were corrected for aliasing and phase offsets and subsequently filtered using a median filter. Inward unit normal vectors were determined on the wall, after which a coordinate transformation was performed for each point at the wall such that the normal vector coincided with the z-axis of the transformed coordinate system. Velocities at fixed points along the normal were calculated using natural 3D interpolation in the original data. Any perpendicular velocity components were removed as only tangential velocity components contribute to the viscous forces at the wall. Smoothing splines were then fitted to the x- and y- velocity components along the inward unit normals. The x- and y-derivatives at the vessel wall were derived analytically and multiplied with the viscosity, which resulted in the WSS. Finally, all WSS vectors were transformed back to the original coordinate system.
This method was validated using a synthetic dataset of a rigid straight tube (diameter 6mm) with parabolic flow, in which the theoretical WSS could be derived analytically (Poiseuille). Effects of resolution, segmentation errors and noise were assessed using this phantom data. Secondly the algorithm was tested in in-vivo PC-MRI data. In vivo PC-MRI of the common carotid artery was performed on a 3T MRI system (Philips Healthcare, Best, The Netherlands) using a dedicated 8 channel carotid coil, 5 heart phases, FOV 80x80x12 mm, isotropic voxel size 0.4 mm (non-interpolated), sense factor 2, Venc of 30 (ap), 30 (rl) and 70 (fh) cm/s, scantime 10 minutes.
This work presents a new method to determine WSS in-vivo independent of assumptions on the flow profile. Results suggest that a resolution of at least 0.6 mm should be used for a 6mm diameter vessel to calculate the WSS within 10% of theoretical values. High-resolution images are needed to avoid segmentation errors. Further validation of the algorithm in-vivo is planned.
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.